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SUMMARY 


* 

} 


This report describes the development of a digital autopilot program for a 
transport aircraft and the evaluation of that system’s performance on a transport 
aircraft simulation. The digital autopilot includes three axis attitude stabili- 
zation, automatic throttle control and flight path guidance functions with emphasis 
on the mode progression from descent into the terminal area through automatic land- 
ing. The study effort involved a sequence of tasks starting with the definition of 
detailed system block diagrams and control laws followed by a flow charting and 
programming phase and concluding with performance verification using the transport 
aircraft simulation. The autopilot control laws were programmed in FORTRAN IV in 
order to isolate the design process from requirements peculiar to an individual 
computer. These control laws were grouped into the following categories: 

• Longitudinal Stabilization (Pitch) 

• Lateral - Directional Stabilization ' 

• Autothrottle Control 

• Vertical Guidance — Non— Landing 

• Vertical Guidance - Landing 

• Lateral Guidance 

Stability, response time and accuracy performance criteria for these various 
functions were identified and the simulations verified that these criteria can be 
met using the specified digital autopilot program. 
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SECTION I 


INTRODUCTION 


This report summarizes work performed by Sperry Flight Systems Division and 
NASA Ames Research Center on a joint study to develop a simulation program for 
an advanced, three-axis, digital autopilot. The effort involved a sequence of 
tasks starting with the definition of detailed system block diagrams and control 
laws followed by a flow charting and programming phase and concluding with 
performance verification using a complete digital simulation of a representative 
transport aircraft. The digital autopilot includes three axis attitude 
stabilization, automatic throttle control, and flight path guidance functions, 
with special emphasis on the mode progression from descent into the terminal 
area through automatic landing. 

The first phase of the study was performed by Sperry. It was a preliminary 
design activity in which all control laws and detailed system block diagrams were 
generated in Laplace Transform format. All filters and compensators associated 
with realistic requirements of airborne equipment were defined. Stability con- 
siderations for all control modes were identified in terms of generalized root 
loci which scoped the stability problems for the jet transport class of vehicles. 
Nonlinear constraints within the aircraft’s flight controls such as actuator 
velocity and travel limits were identified for proper vehicle simulations. 
Appropriate non-linear controls such as maneuvering contraints and control 
authority limits were incorporated in the autopilot control laws in accordance 
with the program objectives of including practical operational considerations. 
Finally, quantitative and qualitative performance criteria and test procedures 
were specified to permit evaluation of the autopilot designs. 

The second phase of the study involved the translation of control equations 
and system block diagrams into the mathematical form and flow charts needed to 
generate the digital autopilot program. This work was done by NASA Ames Research 
Center with Sperry participating in program review and debugging. It is noted 
that neither the scope of the problem or the programming procedures that were 
followed were fully representative of an actual airborne digital autopilot design 
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activity. The process of reducing Laplace Transform specified control laws and 
filters to flow charts and difference equations is certainly applicable to the 
real design procedure. In this study the difference equations were derived by 
converting the Laplace Transform equations to equivalent Z transforms as dis- 
cussed in reference 5. Subroutines were used to compute difference equation co- 
efficients for different sampling rates. In actual practice, when higher 
frequency compensators are programmed (such as those associated with elastic mode 
stabilization) it may be desirable to arrive at the required difference equation 
via the bilinear W transformation to the Z transform in order to improve the 
discrete representation of the continuous filter. (Reference 6). 

A more significant departure from real design practice for an airborne system 
was the use of a high level language to code the guidance and control computations, 
In this study the autopilot control laws were programmed in FORTRAN IV in order 
to isolate the design process from requirements peculiar to an individual com- 
puter. The FORTRAN IV programs permitted the problem to run in nonreal time on 
the IBM 360/367 or in real time on the EAI 8400. While a higher level programming 
language is certainly a convenient approach to the problem, in actual design 
practice, compiler efficiencies do not seem able to compete with direct machine 
language or assembler coding in terms of computer time and memory consumption. 

The most significant difference between the computer programs developed in 
this study and those that would be used in a practical airborne system relates 
to the fact that in this study the problem was restricted to guidance and con- 
trol law computation. A practical airborne system would have to devote a major 
part of its computer program to data formatting, packing and unpacking, I/O 
operations, test, and monitoring routines. Experience in designs (such as the 
SST digital autopilot) shows that these functions would consume 50 to 90 percent 
of the airborne program. 

The organization of this report is based on the autopilot control modes 
rather than on the chronological phases of the study. Six major groups of 
control functions are covered. They are: 

• Longitudinal Stabilization (pitch) 

• Lateral-Directional Stabilization 

• Autothrottle Control ' 

• Vertical Guidance - Nonlanding 
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• Vertical Guidance - Landing 

• Lateral Guidance 

The first two stabilization functions are the autopilot inner loops. The 
guidance functions may be viewed as commands applied to these inner loops. 

Trends in modern control theory tend to neglect this concept of multiple loop 
closures and instead, treat the entire problem in terms of a single state and 
control vector. In this study we have taken the approach of successive loop 
closures because it leads to a better perspective of the design problem. The 
criteria associated with these different loops are considerably different. The 
inner or stabilization loops are concerned primarily with stability and the 
interaction with the actuation system dynamics. In treating these loops we 
must recognize that aircraft rigid body equations are only approximations and 
due consideration must be given to the possibility of elastic mode coupling 
even when the elastic mode data is not available. The guidance modes are, in 
general, uncoupled from the higher frequency modes associated with attitude 
stabilization. 

The criteria for guidance modes are essentially those related to accuracy. 
Guidance (or steering) system design involves such factors as selection and 
I blending of the proper state measurements and the use of appropriate compen- 
sating techniques for winds and aircraft asymmetries. Stability considerations 
are important but different from those associated with the inner loops . While 
the damping ratios of some higher frequency modes associated with the inner 
loops may be tolerated at values of 0.3 or lower, the guidance mode damping 
ratios as high as 0.6 may be objectionable. Relatively low damping ratios 
associated with guidance modes can be made acceptable if special switching logic 
is used to prevent excitation of those modes. 

In the case of the auto throttle modes, the criteria are accuracy and minimi- 
zation of throttle activity. The constraints necessary to prevent excessive 
throttle activity necessitate the use of open loop or predictive controls to 
achieve the desired performance. 
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In this report the individual loops are considered separately and then 
in combination with the other modes. Thus longitudinal and lateral stabili- 
zation functions are treated first without the guidance modes engaged. The 
control laws are developed and stability factors identified. The digital pro- 
grams are described and then performance is verified with the digital simula- 
tions. When this process is repeated for the guidance modes, the attitude 
stabilization and autothrottle modes are already included with their optimized 
parameters. The combination of all control modes defines a complete autopilot 
control law computation program. A brief discussion of mode controller and 
mode interlock requirements that provide for the proper integration and progress- 
ion of the modes is given in the last section. The judicious application of the 
control laws and digital programs described in this report can, therefore, serve 
as a guide to the design of a digital autopilot. The parameters must, of 
course, be optimized for each vehicle and the actual airborne programs must be 
coded efficiently for the specific machine being used. The control functions 
defined herein and the techniques used, however, should be generally applicable 
to any transport aircraft. 
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SECTION II 


LONGITUDINAL STABILIZATION 


A. DESCRIPTION OF CONTROL LAWS 
1 . General 

The longitudinal stabilization functions of the transport autopilot are 
those associated with pitch damping, pitch attitude control, and the pitch 
steering or guidance laws. The general form of these control functions for a 
transport employing powered (hydraulically boosted) elevator control is illustra- 
ted in Figure 2-1 . Pitch control is achieved through an elevator/hydraulic servo 
system and pitch trim is accomplished through a movable horizontal stabilizer. 

Some jet transports operational today use aerodynamically-boosted flight controls. 
The particular aircraft used for the simulations in this study actually employs 
a servo tab system for controlling elevator deflection. Figure 2-1 suggests that 
the elevator servo dynamics are expressible as a linear transfer function with 
appropriate acceleration, velocity and position limits. For hydraulically- 
boosted controls a third order model is usually an adequate representation for 
rigid body simulation and stability studies. A similar third order model is also 
a reasonable representation of the aerodynamically-boosted system. In the case 
of the hydraulic servo a first order lag represents the power boost stage and a 
second order lag represents the secondary actuator (autopilot servo) that strokes 
the power boost control valve. In the case of the aerodynamic boost system, a 
first order lag represents the autopilot actuator while the second order lag 
represents the elevator-to-elevator tab dynamics although the elevator to tab 
dynamics are generally variable with flight condition. The resultant third order 
dynamic representation is, therefore, compatible with both types of boost systems. 
Since the more recent trends in transport aircraft has been toward the use of all 
hydraulically-powered control surfaces, the actuator model used in this study were 
assumed to represent hydraulic systems. 


6 



ATTITUDE 

SYNCHRO- 


MANUAL COMMANDS 
(VARIOUS REFERENCES) 


MISCELLANEOUS • 


L t - RATE LIMIT 


L 2 - ACCELERATION LIMIT 
L 3 - POSITION LIMIT 



BODY AXIS f _ 

1 

• 




PITCH RATE ) . M 6 fC0M 

* 5 . | E ! 

+ l 1 

c 

1 




BU 


v +1 


STABILIZER 

TRIM 

SERVO 



HORIZONTAL STABILIZER 


Figure 2-1 

Basic Pitch Stabilization Block Diagram 
Typical Transport 









2. Servo Response 


Elevator response to autopilot command 




( 2 - 1 ) 


with the following constraints 


\ax = Li 


5e max = L2 
5e max ’ S 

Note that 5 results from force authority limiting if an autopilot parallel 

J^MAX 

servo is used (servo moves the column) . The 5^ would result from authority 

limit stops if an autopilot series servo is used (servo motion is not reflected 
at the stick) . The prevalent commercial transport practice has been to use 
parallel servos although recent trends indicate that series servos are becoming 
more acceptable. When the limit L 3 is associated with a parallel servo, that 
limit is generally a function of flight condition if any aerodynamic load is re- 
flected at the servo or if the feel spring is adjusted as a function of dynamic 

pressure . 

3. Pitch Stabilization 

The basic control law for pitch stabilization is 


V = K 1 + * error] 


(2-2) 


where <j = body axis pitch rate 


and 


(2-4) 


0 = = (0 + B )* 

ERROR' Esc 


4 


th 


L = 9 maximum constraint 

4 E 

6 - SYNCHRONIZED PITCH ATTITUDE 


where 


= (0 - e ± ) 

9 = pitch attitude existing at instant of autopilot engagement. 


(2-5) 


e = A 1 0 G 1 (s) + A 6 G, 
c 1 c 1 1 ^ c 2 


,(s) + ... A n e c G n (s) +« 
n 


(2-6) 


where 


6 = existing value of prior pitch command at time of mode 

c 

° . engagement. 

A to A mode logic = 1.0 for mode ON; 

1 n =0 for OFF** 

» to 9 - various pitch connnands associated with different inodes 

c - % 


1 


G^s) to G n (s) 


= transfer function (including gain) of various pitch 
guidance modes 


.The polarity”conventions that define the polarity of 9, are not standardized. 

_ / vot-Vior t-han subtracts attitude and attitude command, 

ftTei.rr °tha 2 t1h S eTtt r i a tr c^ndtr defined as the difference between 
actual state and reference state as follows. ? 

9 c - f(h - \e F ) 

. fl — (0 — 9 then B would , for the 

If equation 2-4 had been written as ^ (9 V* X c 

above illustration be defined as 0 c = f ^ h REF “ * 

**In nractice the engagement and disengagement of pitch command modes when finite 
error signals exist, must not cause- transient disturbances to the aircraft. 

A function referred to as an easy engage, engage smoother, of fader, s b e 
employed in analog type autopilots to achieve this requirement. The dig a 
autopilot provides this capability by virtue of the attitude command accumula- 
tor function provided by the 9 c term in equation 2-6 and maneuver constraints 

to be discussed in the sections° on guidance law descriptions. 
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The following typical pitch commands are associated with the various 
pitch guidance modes: (Note that the subscripts used are arbitrary and do not 

correspond to the notations used in the description of the guidance laws developed 
in subsequent sections of this report.) 


(1 - cos $ ) 

„ « C 

8 = banking maneuver lift compensation = cQg ^ * 


where 


C a f(Q ) (impact pressure or airspeed) 


6 = roll command 

c 

8 = flap position lift compensation = f(S p ) 

c 2 

0 = throttle compensation = f(8 ) 

C 3 . 

6 = altitude control corrective command = f ( h error » b, M » or v x^ 

C 4 . 

8 = glide path capture steering command = f(beam error, h) 

C 5 ... 

s 

8 - glide path control corrective command = f(beam error, h, 

c e x 

radio alt' 

0 - flare out control command = f(h radiQ , \ nertia i’ 

7 b. , GtC ^ 
oaro 

8 = vertical speed control 

C 8 

0 = airspeed or Mach control 

C 9 


(2-7) 


( 2 - 8 ) 

(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 

( 2 - 12 ) 

(2-13) 


Automatic 
dance with Figure 


trim associated with pitch stabilization is provided in accor- 
2-1 and the following control laws: 


8 =? + B. deg/sec for e > (e + h) for turn ON (2-14) 

u 1 I 1 

8 = + B deg/sec for e > e for turn OFF (2—15) 

Si 

8g “ o for - e.|<e>+e 1 - (2-16) 


_J0 



6g - - B 1 deg/sec for (- e 1 - h) < e 1 for turn ON 


(2-17) 


6 = - B 1 deg/sec for - e <e for turn OFF (2-18) 

(where h is the switch hystersis) 

Note that the threshold € t as shown in Figure 2-1, is based on elevator deflec- 
tion. In autopilot practice where parallel servos are employed, the threshold 
is often based on elevator servo force or displacement of the feel spring. The 
use of an on-off autotrim system of this type is dictated by practical considera- 
tions regarding the stabilizer drive mechanizations. Aircraft that provide pitch 
trim by transferring the control moment from an elevator deflection to horizontal 
tail deflection employ fixed speed trim motor drives. While a proportional 
autotrim mode would have advantages from the standpoint of performance and 
stability, the on-off system is specified for compatibility with the actuating 
mechanisms that are used in practice. 

4. Stability Considerations 

The stability factors associated with the aircraft pitch stabilization 
functions can be derived from the basic block diagram shown in Figure 2-2. The 
body axis rate signal q is assumed to be equal to 0 , the pitch Euler angle rate 
for the purpose of this analysis. Note that the washout (time constant T J is 
added to the pitch rate feedback term in order to eliminate the steady state 
azimuth rate coupling that is sensed by the body axis pitch rate sensor during 
a constant altitude banking maneuver. In Figure 2-2, the aircraft dynamics are 
those obtained from three-degree-of-freedom perturbation equations (forward 
speed, pitching moment, and normal force) . The stabilization of such a plant 
is easily obtained in theory with a simple attitude rate plus attitude displace- 
ment feedback. The phugoid poles are driven into the real axis with a relatively 
low gain attitude feedback. The sum of attitude rate plus displacement creates 
a zero on the real axis that draws the short period poles. This is illustrated 
in Figure 2-3, where servo dynamics and filter effects are neglected. 

This ideal situation does not usually exist in practice. Even the 
phugoid mode which is easily stabilized with simple pitch attitude feedback leads 
to difficulties when pitch attitude is obtained from a vertical gyroscope that 
employs long term gravity slaving to maintain verticality. In that case, a new 
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Figure 2-2 

Block Diagram For Pitch Stability Analysis 
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Figure 2-3 

Idealized Pitch Stabilization Root 
Locus with Servo and Filter Dynamics Neglected 


pair of zeros appears in the open loop transfer function and, as the gravity 
slaving is increased in gain (faster erection to the pendulum) , these zeros 
approach the phugoid poles . The dipoles formed in this manner prevent the 
damping of the phugoid poles. However, this phenomenon tends to be signifi- 
cant only at high cruise speeds where the phugoid period is very long. When 
a tight guidance loop such as glide path control or constant altitude control 
is used, the phugoid mode is eliminated as a source of any problem. Also, the 
use of airspeed control loops through the throttles prevents the excitation of 
a phugoid mode. The point of the gravity erecting gyro has been mentioned here 
because that phenomenon does represent one of the sources of difficulty in some 
automatic approach and landing systems. Erecting gyros, responding to speed 
transients, or correcting for previous speed transients are one of the contribu- 
tors to flight path control errors. Even very small changes in the pitch angle 
reference can cause a few feet of error on final approach. This is a time when 
every foot of error seriously strains the total error budget. . 

The pitch stabilization root locus departs from the ideal form of 
Figure 2-3 when the various filters and servo dynamics are included. Now ex- 
cessive gains will begin to excite high frequency modes created from the upward 
movement of the servo and aircraft short period poles (Figure 2-4). It is 
theoretically possible to design lead compensators that would appear to produce 
a more desirable situation, but this is generally not feasible or desirable be- 
cause of the following factors. Elastic modes (not included in this analysis or 
incorporated in the simulations associated with this program) should be gam 
stabilized with about 10 to 12 db of margin. High frequency compensators aimed 
at correcting servo phase characteristics would result in coupling with elastic 
inodes. Also, servo rate and acceleration limits and nonlinearities are usually 
incompatible with high frequency, high gain compensators. The approach, 
therefore, is to use roll-off filters (such as r 2 of equation 2-3) to ensure 
gain stabilization of elastic modes and to accept the penalty of attitude 
stabilization loops that are not as tight as one could achieve with higher 
gains and idealized vehicle models. 

Note that the stability analysis shown in Figures 2-2, 2-3, and 2-4 does 
not include the automatic stabilizer trim loop. This loop has the effect of an 
integration term. Its linearized representation would add a pole at the origin 
with a zero near the origin. Its effect is usually to draw the short period poles 
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toward the j<0 axis at low gains. The destabilizing effect becomes less sign 
cant at higher gains. The low value of stabilizer trim rate pins a reasons e 
Wide threshold helps eliminate any oscillatory problems as a result of e 
automatic trim loop. 

5. Authority Limit Considerations 

The acceleration, rate and displacement limits (L,, V “d L y associ- 
ated with equation 2 - 1 ) are usually variable with flight condition loading of the 
surfaces. In fully powered systems,, actuator flow limits in both the power 
1st stages and autopilot actuators dictate 1,, the velocity limit. The acce ora- 
tion limits are usually imposed by the static force limits (pressure * effective 
piston area) where the actuator is sized to handle the maximum anticipated su - 
face loads. Control valve dynamics (pressure and flow characteristics in response 
to valve displacements) also enter into this limit. Finally, the displacement 
limit is imposed by physical stops or force limits on the autopilot servos n 
aerodynamic boost systems (control tabs), the force limits of the — ™ 
usually impose the displacement limit L 3 . Figure 2-1 shows servo physical it. 

L , and L S . and an electronic co-and limit L,' • In practice, it is easy to 

vary L • so that the electronic authority limit always corresponds to a reasonabl 
aircraft acceleration limit. A more difficult problem is to achieve a variation 
of L with speed. When electromechanical parallel servos are employed, the torque 
limits on these servos may he adjusted by controlling the current limits as a 
function of aircraft flight condition. Also, the reflection of changing aero y- 
namic loads on the servo or a Q spring feel system serves to change the value o 
L with flight condition. A reasonable approach to simulating this requiremen 
is to compute a continuous value of electronic limit !.,• and assume a second but 
larger physical constraint L 3 for safety purposes. If we select an increments 
"g" limit for commercial transports of + 0.7g's and -O.lg’s, the passengers will 
certainly not object to these g levels. The surface deflection limits may then 
be estimated as follows . 

Let « and 8 represent incremental values from the existing equilibrium 
conditions. Then, to achieve an incremental normal acceleration N z> 


c_ 0 £QS + c 

L a L 5, 


6 QS = - N 
E x g z 


(2-19) 
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Figure 2-4 

Pitch Stabilization Root Locus with 
Servo and Filter Dynamics 





and 



(2-20) 


or 


lflt 


a = 


E 


m. 


a 


Substituting equation 2-20 into equation 2-19 and solving for gives 


S E‘- 


N 

z 


1 

o, B ° 

1 

“ 



W ^ 

J* 

1 

+ C L 
' \ 

QS 


? 2 

for (+)0.7g = 22.5 ft/sec and ('-)0.3g = 9.66 ft/sec contraints 


( 2 - 21 ) 




Note that these limits are functions of dynamic pressure Q. Note also, that 
equation 2—21 represents an approximation to the actual surface deflection that 
will give the specified acceleration limit. The approximation is only as good 
as the accuracy of the linearized coefficients. This accuracy is not adequate 
for defining the specific limit values to be incorporated in an autopilot design. 
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Equation 2-21 scopes the range of limits that must be considered and thereby 
aids in scaling the computation and control designs. The actual limits designed 
into the autopilot control logic should be determined by simulation studies where 
the non-linear aerodynamic effects can be considered more accurately. 

Reasonable limits for , the velocity limit, are 10 to 30 degrees per 
second. Acceleration limits should allow the peak velocity to be attained in 
0.25 seconds. 

6. Longitudinal Stabilization Response Criteria 

There are no specifications or criteria for pitch autopilot response that 
have received wide acceptance. Military handling quality specifications 
(MIL-STD-8785, for example) are applicable to basic airframe response but they 
do not deal with phenomena that are specific to an autopilot. (For example, how 
should the aircraft's pitch response to a pitch command rather than stick command 
be specified.) In this section a set of somewhat arbitrary response criteria 
are defined. They are based on experience obtained with many autopilot designs. 
The suggested criteria are considered achievable and experience has shown that 
when this level of performance is attained, pilot evaluation is good and overall 
guidance accuracy objectives can generally be met. 

/ 

For a step input pitch command, the vehicle attitude response should 
reach 90 percent of the commanded value within about 1.2 seconds. Overshoots 
associated with an oscillatory response should be governed by an equivalent second 
order system damping ratio criterion that permits the damping ratio to have a 
minimum value of about 0.5 (approximately 30 percent overshoot). The pitch re- 
sponse will tend to have a long taxi associated with convergence to its final 
value. If tight airspeed control is maintained, this tail effect is usually 
negligible. A 'reasonable specification on this tail effect is that the response 
converges to within 95 percent of the commanded value wxthin 6 seconds and must 
always remain within 90 percent of the final value within 2.5 seconds. An ideal 
response would be the transient that has the shape of a 0.8 damped second order 
system with frequency of about 4 to 6 radians per second. 
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7. Control Law Parameter Summary 


The range of control law and control system parameters for pitch stabili- 
zation are fairly well bounded for all jet transports. Table 2-1 summarizes 
these parameters in terms of typical values and the possible range of variation 
for all jet transports. 
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TABLE 2-1 


CONTROL LAW AND CONTROL SYSTEM PARAMETER SUMMARY 
PITCH STABILIZATION 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 

W 1 

10 Rad/Sec 

15 Rad/Sec 

20 Rad/Sec 

May represent 1st 
order lag of 
hydraulic power 
servo. 


15 Rad/Sec 

20 Rad/Sec 

30 Rad/Sec 

Represents auto- 
pilot servo or 
tab to elevator 
dynamics . Lower 
values typical of 
parallel servo. 
Higher values for 
series servo. 

L 1 

10 Deg/Sec 

20 Deg/Sec 

30 Deg/Sec 


L 2 

2 

AO Deg/Sec 

2 

60 Deg/Sec 

2 

100 Deg/Sec 


L 3 

±10 Deg 

±15 Deg 

±25 Deg 


L,' 

15 percent 

Per 

25 percent 

f(Q, c ) 

3 

Below 

Nominal 

Equations 
2-22 and 2-23 

Above 

Nominal 

a 

K 1 

2.0 

3.0 

6.0 

Low Q 

1.4 

2.0 

A .0 

Med Q 


0.6 

1.0 

2.0 

High Q 


0.2 

0.5 

1.0 

Ratio of rate to 
displacement . 
These values can 
be increased at 
Low Q and de- 
creased at Hi Q. 

T 1 

2 Sec 

A Sec 

6 Sec 

Washout time con- 
stant - can be 
programmed as 
function of Q 
(increase for 
Low Q) . 
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TABLE 2-1 (cont) 


CONTROL LAW AND CONTROL SYSTEM PARAMETER SUMMARY 
PITCH STABILIZATION 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 

t 2 

0.02 Sec 

0.05 Sec 

0.1 Sec 

Roll-off filter. 

±B 1 

0.03 

Deg/Sec 

0.06 

Deg/Sec 

0.1 

Deg/Sec 

for auto- 

STAB 

trim (lower for 
reduced K^) . 

«1 

0.1 Deg 

0.25 Deg 

0.5 Deg 

Trim threshold 
in equivalent 
degrees of 
autopilot 
command . 

h 

0.05 Deg 

0.08 

0.15 

Hysteresis of 
trim threshold. 

T 3 

(See 

Figure 2-1) 

0.25 Sec 

0.75 Sec 

1.5 Sec 

Trim threshold 

detector filter 
on signal e . 

L 4 

(See 

Figure 2-1) 

±3 Deg 

±5 Deg 

±10 Deg 

Attitude error 
limit. 
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B. DIGITAL PROGRAM - PITCH STABILIZATION 

r'' 

1. Control Law Conversion 

In order to facilitate the programming of the control equations defined 
In the previous section, a notation compatible with the programming language 
must be defined. The mathematical symbols given thus far are convenient for 
analyzing systems and communicating concepts but for machine programming a more 
cumbersome set of symbols must be used. A set of FORTRAN permissible notations 
corresponding to the mathematical symbols used previously was defined. The 
FORTRAN notation has a mnemonic relationship to the original symbol when possible 
A tabulation of the FORTRAN designations and their definitions is given in 
Table 2-2. A general description of the FORTRAN subroutines used for the combin- 
ation of control modes that make up the autopilot is given in Appendix A. The 
equivalent pitch stabilization block diagram using the FORTRAN designations is 
shown in Figure 2-5. Note that this block diagram shows the specific dynamic 
pressure gain control function that optimized performance for the specific 
test vehicle used in the simulations. 

2 . Program Flow Chart 

The Initial condition computations which are performed in the SASIC sub- 
routine for the pitch stabilization mode is given in the following summary and 
the flow chart is shown in Figure 2-6. 


TABLE 2-2 


PITCH STABILIZATION NAMELIST 


Name 

FORTRAN 

Name 

Definition 

q 

QB 

Body axis pitch rate. 

e 

THET 

Pitch angle. 

« E 

DELTE 

Elevator angle; positive-nose down. 

« s 

XIH 

Horizontal stabilizer angle. 

0 

0 

THECOM 

Pitch angle command; positive-nose down. 

5 ec 

DELTEP 

QBW 

Elevator command. 

Filtered body axis pitch rate. 

«sc 

\ 

DHSD 

Horizontal stabilizer rate command signal. 

e 1 

e 2 

BKPT1 

Breakpoint in threshold detector. 

BKPT2 

Breakpoint in threshold detector. 

• 

5 

XMAXIH 

Horizontal stabilizer drive rate. 

6 

IETAB 

Logic switch for horizontal stabilizer drive; 
-1 drive up, 1 drive down, 0 no drive. 

L 3 ( +) 

EL3UP 

Upper elevator command position limit; function 
of dynamic pressure 

l 3 (-) 

EL3DN 

Lower elevator command position limit; function 
of dynamic pressure. 

L 1 

PL1 

Elevator command rate limit. 

Q 

QBAR 

Dynamic pressure. 

Forward pitch loop gain. PSKR = PSK • 57 3 rad/deg 

K 1 

PSK, PSKR 

Kr 

QFBK 

Pitch rate feedback gain. 

T 1 

r 2 

TAUP1 

Pitch rate filter time constant. 

TAUP2 

Pitch rate filter time constant. 

T 3 

TAUP3 

ITESTP 

Automatic trim filter. 

Logic variable for pitch loop synchronization. 
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TABLE 2-2 (cont) 

PITCH STABILIZATION NAMELIST 


Name 

FORTRAN 

Name 

Definition 

\ 

CLALPH 

Lift coefficient curve slope. 

C 

mg 

e 

CMDE 

Elevator moment stability derivative. 

c % 

CMALPH 

Angle of attack moment stability derivative 

% 

e 

CLDE 

Elevator lift stability derivative. 

W 

WAIT 

Aircraft weight. 

S 

AREA 

Aircraft wing area. 

— 

DT2 

Subroutine sample time interval. 

— 

R2D 

Radian to Degree Conversion. 


D2R 

Degree to Radian Conversion . 




AUTOMATIC TRIM CONTROL BLOCK DIAGRAM 


DELTE 



IETAB 


Pitch Stabilization Control 
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Block Diagram Fortran Notation 
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STABILIZER DOWN 

HORIZONTAL STABILIZER 
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(a) Difference equation coefficients for pitch rate feedback filter. 

CQ1 = EXP (-DT2/TAUP1) + EXP (-DT2/TAUP2) 

CQ2 = EXP(-DT2/TAUP1 - DT2/TAUP2) 

CQ3 = (QFBK * R2D/TAUP2)/ (1/TAUP2 - 1/TAUP1) * 

(EXP (-DT2/TAUP 1 ) - EXP (-DT2/TAUP2) ) 

(b) Difference equation coefficients for horizontal stabilizer command 

filter. 

CHI = EXP(-DT2/TAUP3) 

CH2 = 1 - CHI 

(c) Elevator command rate limit 

ERATE = PL1 * D2R * DT2 

(d) Threshold detector limits (converted for input in radians) 

BKPT1 = BKP1 * D2R 
BKPT2 = BKP2 * D2R 

(e) Feedback gain converted for input in radians 

PSKR = PSK * D2R 

(f) Position limits for elevator commands. (Must be multiplied by 
dynamic pressure to get actual value.) 

WSCL - WAIT/ ( (CLAPH/ CMALPH * CMDE + CEDE) * AREA) 

EL3UP = -0.7 * WSCL 
EL3DN = 0.3 * WSCL 

(g) Initialize elevator deflection 

DELTEP = 0 

The chosen pitch control loop parameters for best performance are 
TAUP1 = 4.0 BKP1 = 0.3 

TAUP2 = 0.05 BKP2 = 0.45 

PL1 = 20.0 XIHDOT = 0.06°/sec 


TAUP3 =0.75 


PSKR =5.0 
QFBK = 0.6 




DHSD - CHI • DHSD1 + CH 2 • DELTE 
0HS01 - DHSD 


THRESHOLD DETECTOR f 


HORIZONTAL 

STABILIZER 

COMMAND 

FILTER 


- -1 / \ = 1 
■ IETAB — 


AUTOMATIC 

TRIM 

CONTROL 


PHSD > 0 


PHSD > 0 


PHSD>0 


r DHSD 

. < - bkpti 


F / DHSD 
\ > BKPT2 


DHDS 
< — BKPT2 , 


DHSD 
> BKPTI 


IETAB = 0 


_JLL 

IETAB = 1 


IETAB *= - 1 


IETAB = 0 


BODY AXIS 
PITCH RATE 
WASHOUT 
FILTER 


QBW 

o CQ1 * QBW1 - CQ2 * QBW2 + CQ3 • 


(QB1 — QB2) • (1.096-0.00144 • QBAR) 

QBW2 

= QBW1 

QBW1 

= QBW 

QB2 

= QB1 

QB1 

= QB 


0 c IN DEG 
QBW IN DEG/ 
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DEG 


DELP = PSKR • (THECOM + THET + QBW) • 
(1.17 - 0.00258 • QBAR) 


Figure 2-6a 

Pitch Stabilization Flow Chart 











Figure 2-6b 

Pitch Stabilization Flow Chart (cont) 






C. SIMULATION TEST RESULTS 


The pitch stabilization response was determined using the reference jet 
transport simulator described in Appendix B. The pitch command response at 
the approach condition (V = 141 knots, 50 degree flaps deployed, autothrottle 
engaged and h - 1500 feet is illustrated in Figure 2-7. The response time 
criteria specified in Section IIA-6 are shown shaded. It is apparent that the 
response time objective is met but at the expense of a deterioration in the 
high frequency mode damping. This high frequency mode response can be improved 
theoretically with a lead-lag compensator designed to expand control bandwidth. 
Removing the pitch rate feedback high frequency roll-off filter would also help. 
However, these techniques cannot be considered practical without taking the 
vehicle’s elastic mode characteristics into account. Also, precise definitions 
of the actuator dynamics and sensor dynamics must be included in any such com- 
pensator design. Consequently, the most practical method of improving the 
damping of the high frequency mode is a loop gain reduction of about 20 percent. 
This would compromise the time response. The results shown in Figure 2-7, 
however, are not unacceptable. With a parallel autopilot servo there might be 
some pilot objection to excessive control activity in turbulence. With a series 
autopilot servo, pilot comment on such a system might be an opinion that the 
autopilot is "too tight" but would otherwise be quite acceptable. 

A more interesting phenomenon is observed in the responses illustrated in 
Figure 2-7. The response is characteristic of a statically unstable vehicle but 
the reference jet transport has adequate static margin. The pitch attitude 
overshoots the 5.0 degree reference value after the up-elevator command is 
inserted. The elevator must reverse polarity to hold the new attitude; conse- 
quently the overshoot error must occur. This phenomenon is a characteristic 
that had actually been encountered in the design of the autopilot and autothrottle 
system used in the real vehicle. The problem is the result of the pitch moment 
coupling from the autothrottle system. As pitch attitude increased, the aircraft 
started to decelerate. In order to maintain speed, the auto throttle system 
commanded a thrust increase. The location of the engines below the aircraft c.g. 
resulted in a nose-up pitching moment. This moment must be countered by down 
elevator; hence the sustained overshoot. The overshoot would eventually be 
minimized by the autotrim system. The correction for this problem would have 
been the incorporation of the throttle compensation pitch command per equation 2-9 
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Pitch Stabilization Test Approach Condition 
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Figure 2-8 

Pitch Response of Aircraft to a Step 
Change in Theta Command of w 5° , V = 500 fps 


This compensation was not used for this test. It is noted that in the actual 
autopilot design for the reference jet transport, a compensation term having the 
form of equation 2-9 was included after this phenomenon was encountered in 
flight tests. 

The response at a higher speed condition prior to entering final approach 
is illustrated in Figure 2-8. At this condition, the speed is 296 knots, clean 
(zero flaps) and autothrottle is not engaged. The pitch response meets the 
time and damping criteria. Note that the elevator position limit is reached and 
consequently the initial response is slowed somewhat by saturation. Also, the 
absence of the auto throttle pitch moment coupling prevented the type of overshoot 
seen in Figure 2-7. 



SECTION III 


LATERAL STABILIZATION 


A. DESCRIPTION OF CONTROL LAWS 
1 . General 

The lateral stabilization functions of the transport autopilot are those 
associated with yaw (or dutch roll) damping, turn coordination, and roll attitude 
control. The dutch roll damping is provided by a yaw damper that drives the rud- 
der in response to various sensor inputs. For many flight conditions, a roll 
damper is also an effective dutch roll damper. The use of a separate roll damper, 
however, has significance only for manual control. Roll damping is implicit in 
the roll attitude control function and is not treated as a separate stabilization 
loop in the subsequent description of lateral stabilization control laws. 

Turn coordination is considered as part of the yaw stabilization function. 

In addition to providing for damping of dutch roll oscillations, the yaw (or rud- 
der) channel is also used to turn the aircraft into the relative wind for coor- 
dinated (zero sideslip) turns. (Note that zero sideslip and zero lateral 
accelerometer or pendulum angle are not identical. The differences are illustrated 

in Appendix C. 

Roll attitude stabilization is achieved through the actuation of rolling 
moment producing surfaces, generally ailerons. Many of the jet transports use 
combinations of ailerons and spoilers. Supersonic vehicles use differential tail 
controls or elevons where roll and pitch controls are accomplished with common 
surfaces. The reference vehicle on this study uses inboard and outboard ailerons 
that are controlled by aerodynamic tabs. For aerodynamically-boosted control, 
the control surface deflection is determined by the hinge moment equations. Rolling 
moment characteristics associated with ailerons and spoilers tend to be quite non- 
linear and very dependent upon other variables such as angle of attack, flap deflec- 
tion, and sideslip. These nonlinear effects could easily be masked by high gain, 
roll rate, and roll displacement feedback loops, but the problem is usually com- 
pounded by mechanical nonlinearities in the control system (backlash, cable stretch, 
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friction, etc). In large transports, these mechanical nonlinearities often 
determine the quality of an autopilot's roll control performance. In the present 
study, an idealized control actuation system is assumed. It is represented by 
realistic, linearized dynamic transfer functions. This linear representation of 
the control system will permit performance that is, in general, better than that 
which is obtainable in the real environment of control nonlinearities. The roll 
performance results obtained in this study should therefore be viewed as an upper 
bound or goal that may be difficult to achieve with state-of-the-art mechanical 
and hydromechanical controls . 

V 

2. Yaw Stabilization and Turn Coordination 


a. Feedback Variables 

The yaw stabilization and turn coordination control laws are those 
associated with the 5 R (rudder) feedback on Figure 3-1. These feedbacks include 
the following information: 

* r » body axis yaw rate (yaw rate gyro) 

m a lateral acceleration (eg, mounted accelerometer) 

• <p = roll angle 

*8 = commanded aileron deflection 

AC 

• 8 =» rudder deflection 

R 

# => skid command (from lateral steering system) 

These variables may be used in different combinations ranging from the simple, 
single variable (r) yaw rate dampers to the use of all of the above variables in 
a computed P (sideslip rate) damper . 
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b. Yaw Rate Damper 


The typical jet transport dutch roll mode is readily damped with yaw 
rate feedback as illustrated in Figure 3-2a. The [ r/5 R ] (s) aircraft dynamics 

^^ire of the form: 




(s) 


( r R s + ') (V + 1 ) (— 2 

V°D 



(3-D 


where 


x = is the roll convergence time constant 
R 


s 


is the spiral mode time constant (usually slightly negative) 

= is the dutch roll frequency 

t «= is the dutch roll damping ratio 
D 

The location of the numerator quadratic zeros determine the effectiveness of the 
yaw damper in damping the dutch roll. When approaches ^ to form a dipole, 

^ he yaw damper becomes ineffective. This often occurs at high angles of attack 
^for some transports and is also a common phenomenon in vehicles designed for 

hypersonic flight. The optimum yaw damper gain is usually the one yielding maxi- 
mum damping, as illustrated on Figure 3-2a. Higher gains cause the dutch roll 
period to be stretched with the response becoming more oscillatory. 

It is noted that the simple yaw rate-to-rudder feedback described in 
Figure 3-2a also provides good spiral mode stabilization. It will be demonstrated, 
however, that this spiral mode stabilization effect is not generally attained when 
some of the practical problems begin to be considered. The first problem (not 
revealed in the root locus diagrams) relates to achieving good turn coordination. 
In a steady state, coordinated turn, the steady body yaw rate would deflect the 
rudder to cause the aircraft to develop sideslip in a direction that drives the 
aircraft out of the turn. That is, if the yaw damper control equation is 


+ 6_ = k.r 


(3-2) 
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and the steady yaw rate (for a constant altitude, coordinated turn) is 

r = y sin 0 cos 9 

the. steady rudder deflection would be 

+ S R = ^ sin 0 cos 9 


(3-3) 


(3-4) 


This is, of course, unacceptable. The standard practice for coping with this 
problem has been to incorporate a washout function in the yaw rate feedback con- 
trol law. Thus, in terms of the control law functions identified on Figure 3-1, 
the simple yaw rate damper with washout to eliminate steady turning rates would 

have the following control law 

+ « R ar G t (s) G 4 (s) ys) (3-5) 


where 


G.j (s) 


G 4 (s) 


H r(3) 


lag filter = 

V 

washout « k 2 

rudder servo transfer function 


(3-6) 

(3-7) 



(3-8) 


Neglecting the rudder servo dynamics and the lag filter (for elastic 
mode and noise attenuation), the root loci for the yaw damper with washout (j J 
are illustrated on Figures 3-2b and 3-2c. Note that for fast washout time con- 
stants, the attainable damping improvement decreases. Fast washouts are desirable 
for sideslip minimization following turn entry and exit, but they are undesirable 
because of the loss in attainable damping. One method often used to help minimize 
sideslip, both in turns or in response to gust disturbances, is lateral accelera- 
tion feedback. 

c. Lateral Acceleration Feedback 

The transfer function of lateral acceleration' for rudder deflection 

inputs is 


x -r -j 





(3-9) 



(v + ’) (v + 0 ( r c s + 0 (v + 0 

(v + ') ( t r s + ’) ft ' + zf + ') 


where typically 


T «T 
A S 


t ~ r 
B R 


T and r are approximately equal and of opposite polarity 

v U 


Of the two possible polarities for feedback, one provides sideslip 
minimization, but usually at the expense of dutch roll damping. The other polar- 
ity of A feedback can damp the dutch roll, but it will cause sideslip divergence 

y 

by driving the spiral mode pole deeper into the right-half plane. This is illu- 
strated in Figures 3-3a, and 3- 3b, which show the root loci for both polarities 
of feedback. Figure 3-3a, the correct control for sideslip minimization, causes 
a deterioration in dutch roll damping, but the yaw rate feedback associated with 
the complete control law should have adequate gain to compensate for the adverse 
shift in the dutch roll poles. The proper polarity for sideslip control is: 


Right rudder command for aircraft acceleration to the left. 
(Acceleration to the left is caused by -40.) 


d. The P Damper 

In large jet transport aircraft, the problem of achieving adequate 
dutch roll damping without sideslip deterioration is aggravated by the size of 
the aircraft. The larger values of yaw and roll inertias result in longer period 
dutch roll oscillations. Washout time constants on the yaw rate feedback should 
be sufficiently high so that the proper phase of the yaw rate feedback exists at 
the dutch roll frequency. For example, a 12.56-second (w = 0.5 rad/sec) dutch 
roll period may not be adequately damped by a yaw rate feedback signal with a 
2.0-second washout. The phase of the yaw rate signal would be the phase of 
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The desirable phase at the dutch roll frequency is zero degrees; but in the above 
case, it would be +45 degrees. 


One method of coping with this problem is to account for the steady 
yaw rate during a bank angle maneuver in some manner other than with a washout 
filter. An intuitive solution would be the establishment of a dynamic yaw rate 
reference other than zero. That is, let the yaw rate reference be zero only when 
the aircraft is in nonbanking flight. The yaw damper feedback signal would there- 
fore be r error wh ere 


r error r r REF . 


(3-10) 


r REF. 


^ sin 0 


(3-11) 


Such an approach is often referred to as the P damper because subtracting the 
steady turn rate from the damping signal approximates a damping signal proportional 
to sideslip rate. In principle, therefore, the damper would be responsive to 
motions with respect to the relative wind rather than to an inertial reference 
frame. Even if perfect compensation is not achieved, a longer time constant (G^(s) 
in Figure 3-1) could be used without causing the adverse yaw sideslip effect en- 
countered with the conventional yaw rate feedback plus long washout time constant. 


The P damper is not commonly used because of its computation complex- 
ity. It should be noted that yaw dampers are generally simple devices that do not 
contain the computation sophistication of an autopilot. There is also a problem 
regarding access to a bank angle signal and to a velocity signal since yaw dampers 
must usually be autonomous of other aircraft subsystems. The perfect P damper 
equations can be derived from the small perturbation side force equation (in body 
axis coordinates) 

S F=mV/3+r-ap--§-0 • (3-12) 

y o V 


or 

P =» — — ^ -r + ap+^0 (3-13) 

mV o V 

where 

F is the sum -of aerodynamic and thrust forces. 

y 
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Since a lateral accelerometer is responsive to the sum total of aerodynamics and 
thrust forces on the aircraft (see Appendix C) , we can substitute for F y using 
the following relationship 



mV 



(3-14) 


where 

I 

A is the effective acceleration sensed by a lateral accelerometer. 

y 

Substituting equation 3-14 into equation 3-13 and solving for (-0), 
the correct polarity for a damping feedback term yields 

HS - r - f t - ^ - V 0-15) 

Note that the first two terms represent the simple or approximate P damper defined 
by equations 3-10 and 3-11. Adding the lateral accelerometer with its sensitivity 
inversely proportional to V helps converge the result to perfect P . The a^p term 
can be significant at landing angles of attack although roll rates. p, for transports 
in landing approaches are usually kept low. Feeding back perfect P (even if it can 
be computed) has one slight disadvantage. (Note that the polarity of the A y term 
is opposite that required for sideslip minimization, as discussed in the previous 
paragraph . ) 
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The root locus for a typical P damper is shown on Figure 3-4. It is 
seen that unlike the typical yaw rate damper, the spiral mode divergence may be 
increased slightly. Damping capability for the dutch roll mode, however, is 

excellent . 


e. Other Sideslip Control Techniques 


Major sources of sideslip in turning maneuvers are the yawing moments 
and side forces associated with turn entry and exit. Yawing moments due to ailer- 
ons, spoilers, and roll rate cause sideslip. The Q^p term in equation 3-13 can 
result in large sideslip transients if the roll rates are high. In jet transports, 
these dynamic turn coordination effects are minimized by restricting roll rate 
maneuvers to relatively low values. If maneuvering roll rates above about 10 
degrees/second are allowed and if large deflections of ailerons (or spoilers) are 
permitted, then a corrective, feedforward compensation proportional to aileron 
displacement should be fed into the rudder channel. 

As shown on Figure 3-1, the compensation term is: 


8 



G 6 (s) 


(3-17) 


where 


g 6 <3 > - k 6 


r t 3 s i 


1 

T^S + 1 


r s + 1 
_ 4 


(3-18) 


The required compensation gain can be computed approximately as follows: 


%A S A + "SR \ - 0 


(3-19) 


8 


R. 


n sa 5 a 

N c 


“ k 6 5 AC 


(neglecting washout) 


or 


k 6 



(3-20) 


(3-21) 
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f. Rudder Centering 


A typical yaw damper may develop offset errors as a result of_.the 
k accumulation of electronic and electrohydraulic unbalances. To insure a zero 
steady state rudder deflection for undisturbed flight, the position feedback 
associated with the rudder servo mechanization may incorporate an integration 
term. This will provide the effect of a washout on all rudder commands. It is 
represented on Figure 3-1 as G 5 (s). This long term washout should be at least 
ten times larger than the yaw damper washout time constant (r . In Figure 3-1 , 
the rudder servo dynamics H^s) are represented as having unity static gain. 

Hence, the effective long term washout time constant will be the reciprocal of 
the integrator loop gain. Thus, 

^5 

g 5 (s) * “ (3-22) 

and 

T « —1 (3-23) 

x k^ 

where T is the effective long term washout time constant on rudder commands. 

, Since this function is added to cope with specific equipment mechanization problems 
it need not be included in the present study; hence k 5 should be assumed to be 

zero. 

3. Roll Stabilization 

With an effective yaw damper operating through the rudders, roll attitude 
stabilization is usually a very simple task if a linear control system can be mech- 
anized. Reference to the roll- to-aileron transfer function illustrates this point: 



for the undamped vehicle. This transfer function retains the same form (in the 
lower frequency regions) when a yaw damper loop is closed through the rudder. 



However, the yaw damper loop closure moves both the poles and zeros of equation 
3-24 into regions that considerably simplify the achievement of a tight roll 
attitude system. With yaw damper operative, the roll transfer function becomes 



Where the notation ( )' represents the new value due to yaw damper loop closure. 
The root loci for a roll attitude control law of the form 

5 A = k 7 + G 7 (s)p ] V s) (3 “ 26) 

are shown on Figure 3-5. G 7 (s) is a cutoff filter that may be approximated by a 

first order lag and H (s) may be represented by a third order model identical to 

A 

that of equation 3-8. On Figure 3-5, the lag in G 7 (s) is neglected and it is 
assumed that 

G 7 (s) = a 1 (3-27) 

and 

p = s0 (3-28) 


Note that the effect of the yaw damper loop closure is to move 
and 0 into well damped regions where they form a dipole that is essentially un- 
coupled from the roll stabilization dynamics. Also, the spiral root has been 
moved from the right-half plane to near the origin. However, even if the spiral 
root remained in the right-half plane, the closure of a roll loop rapidly moves 
the spiral pole into a stable region. 





Figure 3-5 

Roll Stabilization Root Loci with Yaw Damper Loop Closed 


4. Lateral Stabilization Control Law Summary 


The recommended control laws are*: 
Rudder 


+ = (r - r ) 

RC c 


T 2 s 


T 2 s + 1 


r ^s + 1 


k 2 +A y 


1 

r 5 s + 1 


+ 5 


AC 


T 3 S 


r 3 s + 1 


t.s + 1 
4 


where 


r c “ V Sln 0 c ( 573 ) deg/sec 


Aileron 


(3-29) 


(3-30) 


* Ok + &1 r p] (3-31) 

^ AC k 7 [j* error t 7 s + 1 J 

Typical values for lateral stabilization parameters are summarized in Table 3-1. 


*The polarity of the surface commands depend upon the polarity convention used 
for the vehicle equations. There is no universally accepted convention for 
these polarities. The rudder and aileron control laws given here assume that 
positive rudder is trailing edge to the left and positive aileron is left wing 
aileron trailing edge down (right rolling moment) . 
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TABLE 3-1 
PARAMETER SUMMARY 



0.4 
0.05 

5 deg/scc 
3 deg/sec 

♦Values are given for mechanical Halt 


0.8 
0.1 

10 deg/see 
5 deg/sec 

Electrical liolt should be set 
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5. Lateral Stabilization Performance Critetia 


a. Procedure for Establishing Yaw Damper 

• Establish free-aircraft simulation and obtain transient 
response using initial condition of 2.0 degrees in sideslip 

• Close yaw rate feedback loop without washout 

• Increase yaw rate gain until damping ratio of greater than 
0.5 is obtained, or until damping begins to deteriorate with 
long-period oscillation detectable, or until signs of high 
frequency mode appear on. 5 R response. The S R response should 
never show oscillations with damping below 0.20. 

• Add yaw rate washout time constant, starting with value of 
about 4.0 seconds and decreasing until dutch roll damping ratio 
shows deterioration to below about 0.5. 

b. Procedure for Establishing Roll Stabilization Parameters 

• Use nominal control law - § A = k (0 g + 0.5 p) . 

• Apply 5-degree step roll command (with roll rate constraint L^) . 

• Raise k to nominal or higher value until roll overshoot is held 
to less than 10 percent. Roll angle should settle at final 
value within 3.0 seconds. If any evidence of high frequency 
oscillations appear on the § A response, the gain k should be 
lowered. High frequency § A motions with damping ratios below 
0.2 are objectionable. 

t 

c. Turn Coordination 

• Apply 20-degree bank command, 0 c » through a 1.0-second lag 
filter plus a command constraint that restricts 0 fi to 10 deg/sec. 

• Observe peak excursions of 0 and A^. The lateral acceleration 
A should not exceed O.lg (3.2 ft/sec 2 ) during this maneuver. 

y 



• Observe the effect of the yaw rate washout time constant on P 
and A during the bank maneuver. If reduction in washout time 
constant improves turn coordination but causes a dutch roll 
damping deterioration, then the P damper may be desirable or 
the inclusion of an loop should be considered. 

• If transient turn coordination is poor because of yaw moments 
due to aileron, then aileron compensation loop should be added. 

The criterion should be O.lg peak lateral accelerometer reading 
for the most severe bank maneuver. Steady state g’s in the^tum 
for 20-degree bank should always be below 0.5 to 1.0 ft/sec . 

B. DIGITAL PROGRAM - LATERAL STABILIZATION 
* v 1. Control Law Conversion 

The representation of the lateral stabilization block diagram and control 
laws in FORTRAN notation is summarized in Figure 3-6 with a tabulation of the 
FORTRAN namelist given in Table 3-2. Note that, as in the case of the pitch sta- 
bilization control laws, a proportional gain program as a function of dynamic 
pressure Q is incorporated. 




Figure 3-6 

Lateral Stabilization Block Diagram 
with Fortran Notation 


/ 













TABLE 3-2 


VARIABLE NAMELIST 
SUBROUTINE LATSTB 


Variable 

FORTRAN Name 

Definition 

P 

PB 

Body axis roll rate 

0 

PHI 

Euler roll angle 

«c 

PHICOM 

Roll angle command 

r 

RB 

Body axis yaw rate 

A 

y 

AYACC 

Lateral body-mounted accelerometer output 

V 

VT 

Aircraft ground speed 

S R 

DELTR 

Rudder angle 

5 rc 

DELTRP 

Rudder angle command 

5 ac 

DELTAP 

Aileron angle command 

a 1 

LAI (ULA1 ) * 

Roll rate gain 

K 7 

ULK7R (ULK7) 

Roll stabilization gain 

K 2 

LK2 

Yaw damper gain 

*3 

LK3 

Accelerometer feedback gain 

K 4 

LK4 

• 

Roll angle feedback gain; P damper 

K 5 

LK5 

Rudder centering integrator gain 

K 6 

LK6 

Aileron feedback gain 

S' 

L2P (PIP) 

Aileron command rate limit 

Y 

L2P (P2P) 

Aileron command position limit 

S' 

L3P (P3P) 

Rudder command rate limit 

S’ 

UP (P4P) 

Rudder command position limit 

T 1 

TAUL1 

Yaw damper filter time constant 

T 2 

TAUL2 

Yaw damper filter time constant 

T 3 

TADL3 

Aileron compensation filter time constant 

T 4 

TAUL4 

Aileron compensation filter time constant 


[ 53 . 


*Fortran variables in parenthesis are unsealed values. 








TABLE 3-2 (cont) 

VARIABLE NAMELIST 
SUBROUTINE LATSTB 


Variable 

FORTRAN . Name 

. . .Definition 

r 5 

TAUL5 

Accelerometer filter time constant 

T 7 

TAUL7 

Roll rate filter time constant 


DT2 

Subroutine sample time 

- 

R2D 

Radians-to-degrees conversion 

- 

D2R 

Degrees— to— radians conversion 

- 

ITESTL 

Logic variable for lateral stabilization loop 
synchronization 

Q 

QBAR 

Dynamic Pressure 


2. Program Flow Chart 

The initial condition computations which are performed in the SASIC sub 
routine (see Appendix A) for the lateral stabilization mode is given in the 
following summary. The flow charts are shown in Figure 3-7. 

Lateral Stabilization Subroutine LATSTB IC Ca lculation 

a. Set Outputs 

SKIDCM = 0 
DELTRP =0 
DELTAP = 0 

b. Scale Limits 

LIP = PIP * D2R * DT2 
L2P - P2P * D2R 
L3P = P3P * D2R * DT2 
L4P = P4P * D2R 

c. Scale Gains 
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ULK7R = ULK7 * D2R 
LK3 = ULK3 * D2R 
LK4 - 32.2 * LK2 * NOP 


/ 





d. Difference Equation Coefficients 


• Roll Rate Filter 

CPI = EXP (-DT2/TAUL7) 

DPI *» LAI * (1.0 - CPI) * R2D 

• Aileron Compensation Filter 

CEP1 « EXP(-DT2/TADL3) + EXP(-DT2/TAUL4) 

CEP2 = EXP[-DT2 * (1.0/TAUL3 + 1.0/TAUL4)] 

DEP1 = LK6/TAUL4 * [ EXP(-DT2/TAUL3) -EXP(-DT2/TAUL4)] / 
(1.0/TADL3 - 1.0/TAUL4) 

• Rudder Centering Integrator 
LK5 ULK5 * DT2 

• Yaw Rate Filter 

CRB1 = EXPC-DT2/TAUL1) + EXP (-DT2/TAUL2) 

CRB2 - EXP[ -DT2 * (1.0/TAUL1 + 1.0/TAUL2)] 

DRB1 = LK2/TAUL2 *[ EXP (-DT2/TAUL1) - EXP (-DT2/TAUL2)] / 
(1.0/TADL2 - 1.0/TADL1) 


• Lateral Accelerometer Filter 

CA1 - EXP (-DT2 / TAUL5 ) 

DAI - LK3* (1.0 - CA1) 


The chosen Lateral Stabilization parameters for best performance are: 

3.0 


-3.0 


TAUL1 

= 0.1 

PIP = 

25.5 

ULK7 « 

TAUL2 

= 2.4 

P2P = 

14.9 

LK2 = 

TAUL3 

TADL4 

= 2.0 
■ 0.2 

P3P - 
P4P = 

17.0 

15.3 

ULK3 ' 


TAUL5 

TAUL7 


0.15 

0.1 


LK6 

LAI 


-0.15 

0.5 
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EPOUT - CEP1 * EPOUT1 - CEP2 . EPOUT2 + DEP1 • (EPIN1 - EPIN2) 

EPOUT2 » EPOUT 1 

EPOUT 1 - EPOUT 

EPIN2 - EPIN1 

EPIN1 “ DELTAP 


[Aileron compensation filter!] 

[turn COORDINATION J 


" RBOUT - CRB1 • RBOUT1 - CRB2 • RBOUT2 + DRB1 • (RBI - RB2) 
(0.6531 - 21 .25133/Q BARKj: 

RBOUT 2 - RBOUT1 
RBOUT 1 - RBOUT 
RB2 - RBI 
PHICMR - PHICOM ♦ D2R 

RBI i RB — (32.2/VTI « SIN (PHICMR) — 


[YAW RATE FILTER: YAW DAMPER] 


AYOUT = CA1 • AYOUT + DAI • AYIN • (-0.15624 + 70.8371/QBARK) 
AY ACC - AY - 32.2 • SPHI 
AYIN » AY ACC 


SDOUT - SDOUT + LK5 • SDIN 
SDIN - SKIDCM - DELTR 



[lateral accelerometer! 
[turn coordination J 


[rudder centering] 
[integrator J 



E jdder command] 
IMMER J 


Figure 3-7b 

Lateral Stabilization Flow Chart (cont) 








Figure 3-7c 

Lateral Stabilization Flow Chart (cont) 
































C. SIMULATION TEST RESULTS 


1 . Yaw Damper 

Yaw damper performance achieved with the recommended control law parameter 
is illustrated in Figure 3-8. The transient responses are at the 141-knot and 
295-knot, sea level flight conditions. They illustrate sideslip in response to a 
step gust corresponding to 2.0 degrees of sideslip. It is apparent from these 
responses that dutch roll damping (by the yaw damper) is excellent at the higher 
speed condition but only marginally acceptable at the low speed condition. The 
relative ineffectiveness of the yaw rate into rudder loop as a dutch roll damper 
at the landing approach flight condition is a common aircraft phenomenon. It is 
a consequence of the close proximity of the zeros to the ^ poles in equation 
3-1. Fortunately, this does not cause any autopilot problems because the dutch 
roll is readily damped at this flight condition by a roll rate feedback into the 
ailerons. This feedback is inherent in the roll stabilization system. 

2. Roll Stabilization 

The absence of any dutch roll oscillations and hence the effectiveness of 
the roll stabilization loop as a dutch roll damper at low speeds is apparent from 
the roll stabilization test results illustrated in Figure 3-9. Roll angle re- 
sponses to filtered 5.0-degree bank angle commands are shown for the 141-knot and 
296-knot flight conditions. Damping is excellent and static accuracy is perfect. 

3. Turn Coordination 

The turn coordination criteria defined in Section III A-5 were easily 
achieved as illustrated in Figures 3-10 (a) and 3-10(b) for the low and high speed 
conditions respectively. Those criteria can probably be viewed as worst case 
numbers. The level of performance achieved in Figures 3-10 (a) and 3-10 (b) would 
probably be rated good to excellent by a flight test pilot who would evaluate 
turn coordination performance on the basis of his ball-bank inclinometer. One 
"ball" on this instrument is equal to a 4.7-degree deflection of the apparent 
gravity angle. The apparent gravity (or pendulum) angle, o t is 


(—Sr) 

\cos w 


—2. cos 0 
g 





a = 



SIDESLIP (DEG) 0- SIDESLIP (DEG) 


YAW DAMPER PERFORMANCE: V- 141 KNOTS 



YAW DAMPER PERFORMANCE: V - 296 KNOTS 



SECONDS 

YAW DAMPER PARAMETERS: 

TAUL1 = 0.1 SEC 
TAUL2 » 2.5 SEC 

LK2 - -2.0 (THE YAW DAMPER GAIN IS REDUCED 
AUTOMATICALLY AS THE DYNAMIC 
PRESSURE INCREASES) 

P3P - 17.0 DEG/SEC 

P4P - 15.3 DEG 


Figure 3-8 

Yaw Damper Performance 




0 ROLL ANGLE (DEG) 0 ROLL ANGLE (DEG) 







A pilot would consider turn coordination good if his inclinometer shows less than 
1/8 ball. At 0 “ 30 degrees, 1/8 ball corresponds to 0.0118 g's or 0.38 ft/seC . 

At the low speed condition the transient acceleration reaches a peak of 0.039 g s 
but always reduces to less than 0.01 g's steady state. The rudder control system's 
turn coordination includes the lateral acceleration (gain controlled with Q) and 
the aileron-to-rudder crossfeed. Any additional sophistication to improve turn 
coordination performance does not appear to be warranted. 
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TURN COORDINATION TEST: V- 141 KNOTS 



\ 




t (SEC) 



10 20 30 40 

t (SEC) 


TURN COORDINATION PARAMETERS: 

TAUL5 = 0.15 SEC 
TAUL3“ 4.0 SEC 
TAUL4 » 0.1 SEC 
LK6 “0.1 

ULK3 - 3.0 DEG/FT/SEC 2 

YAW DAMPER AND ROLL STABILIZATION PARAMETERS GIVEN IN PREVIOUS TESTS. 


Figure 3- 10a 

Turn Coordination Performance 
(141 Knot Condition) 
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SECTION IV 


AUTOTHROTTLE SYSTEM 


A. DESCRIPTION OF CONTROL LAWS 


1 . General 

The general autothrottle system block diagram is shown in Figure 4-1. 

The basic inputs are calibrated airspeed (v.) in knots, a commanded airspeed 
reference /v \, compensated inertial acceleration (x), pitch angle (9), and 

\ c ref7 ^ 

elevator position (fig). Compensated forward acceleration refers to the sub- 
traction of the component of gravity due to pitch angle from the body-mounted 
fore-aft accelerometer. This resolution of the required acceleration signal 
from an accelerometer is shown in Figure 4-2. The cosine of the angle-of-attack 
compensation is not required for the usual range of angle of attack. 

2. Control Laws 

~ t , 

The complete throttle control equation may be expressed as 


V 

C ERROR 





+ 5 


E 




(4-1) 


where 5 

P 


is a predictive feedforward command based on computation of 


future throttle requirements 


T 

C— “ 

*T 


- G 4 (s) 


pL + i1 

h J 


Throttle Servo Dynamics 


(4-2) 


with following constraints 


6 = L_ 

t max * 


(4-3) 


W L 3 


... 66 ... 


(4-4) 



WEIGHT ► 

AIRSPEED 

REFERENCE 


Sf 


COMPUTATION 


FLAP 




POSITION 


COMPENSATED 

FORWARD 

ACCELERATION 


V„ 

e REF 


AIRSPEED 

ERROR 

COMPENSATOR 




RATE AND 

DISPLACEMENT 

LIMITS 






THROTTLE 

SERVO 

DYNAMICS 


Figure 4-1 

Autothrottle System Block Diagram 







ACCELEROMETER READS 
. —3.22 FT/SEC2 (DECELERATION) 
BECAUSE OF PITCH COUPLING 


Figure 4-2 

Pitch Coupling In Fore-Aft Body 
Mounted Accelerometer 
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It is noted that most mechanizations of autothrottle servo systems use 
velocity servos rather than position servos (equation 4-2). This is a practical 
consequence of the requirement that the throttle servo have a floating reference 
position established by the manual throttle input. The velocity servo mechaniza- 
tion usually eliminates the need for specifying an integral compensation term. 

For the purposes of this study, however, the mathematical form of the control laws 
as throttle increments will be valid and all peculiarities of mechanization will 
be omitted. 


The combination of G 1A (s) and G 1B (s) includes a complementary filter that 
creates a wide bandwidth speed error signal using inertial data for short term 
and pressure data for long. A detailed block diagram of this function is shown 
in Figure 4-3. As seen in this figure 


and 


g ia ( s) 


1 

T S + 1 
C 



“v 


G 1B^ “ T 


K T 
v c 

S + 1 


+ K. 


(4-5) 


(4-6) 


where provides for a lead compensator, if required. 
The compensated signal x^ is defined as 


X ACCEL- g Sl " 6 " *C 


(4-7) 


and 


•• • 

x • a V cos Of + g sin 9 
X ACCEL 6 


(4-8) 


where V is true inertial acceleration along the velocity axis. A problem is 

often encountered because of imperfect cancellation of the pitch angle coupling 

term. The higher the complementary filter time constant T the higher the 

accelerometer loop gain. Any erroneous null signal including the imperfect 

measurement, scaling and subtraction of g sin 0 can cause large errors in the 

airspeed control loop accuracy. Hence, in practice, a washout prefilter is often 

added to the x filter. This filter produces a new signal x ' of the form 
c 


i 



X T 

c w 


T S 
W 


+ 1 


(4-9) 
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where r is about 100 seconds. Since this is a problem encountered in actual 
w 

aircraft installations but need not appear in simulator studies (unless one 
wishes to simulate instrument installation accuracies, etc.), there is no reason 
^to include this filter in a study of control laws. 

The elevator compensation control, G^(s) is 
G 3 (s) - 5 e 

The washout is needed to remove steady state elevator deflections from the 
control signal. Time constants in the order of about 30 seconds are needed. 

The pitch compensation G 2 (s) has a similar form as G^Cs). If T g above is 
sufficiently large, the elevator compensator actually provides the same speed 
change anticipation that would be obtained from the pitch signal. Both are shown 
for generality, but it is quite possible that only one would suffice. The 9 or 
6^ control inputs act as feedforward predictors that allow the throttle loop to 
begin correcting for an anticipated error. That is, the correction occurs before 
the error develops. In aircraft where the engine thrust change produces a pitch- 
ing moment, the possiblity of coupling the pitch-to-throttle compensation and 
^throttle-to-pitch compensation into a throttle-pitch instability exists. The 
filters on the compensation signals must be set for the proper phasing and 
attenuation that avoids this type of instability. 


V 

T ? S + 1 


V + 1 


(4-10) 


An important computational requirement for an autothrottle control law is 

the definition of the calibrated airspeed reference command, V . Even if the 

REF 

speed reference is selected on a display device by the pilot, the effective command 
to the autothrottle system must be processed for inclusion of a rate constraint 
before being applied to the autothrottle control law. In general, a command rate 
limit, often referred to as the retard rate limit, should be about 1 knot per 
second. If the airspeed reference program is automated to follow the correct 
speed requirements as flaps are deployed, the actual value of V c can be determined 
using the following procedure on a simulator: 


• Hold throttles fixed. 

• Deploy flaps at standard flap rates. 

• Constrain flight path to desired trajectory (hold aircraft on glide 
path) . 
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• Measure actual V 


• Use this value of V as V , and the final value of calibrated airspeed 


as the new V 


"REF 


"REF 


• The V program will therefore be 
c 


'REF 


V , V 
C o + C REF 


(4-11) 


until the desired final value is reached, 
final value. 


The V 


is held at that 


"REF 


The effect of this type of V programming is to allow the aircraft to 

C REF 

decelerate in response to flap deployment without necessitating any throttle 
adjustments. In effect, the airspeed error is maintained zero by predicting the 
airspeed transient and using that transient as the reference airspeed. 

The reference airspeed may be adjusted as a function of weight. A 
method that allows the correct reference speed to be determined for any aircraft 
weight uses an angle-of-attack outer loop on the as follows 


'REF 


"REF. 


+ k «/(“ -“e ep) 


dt 


( 4 - 12 ) 


where a 


REF 


’NOMINAL 

is the desired or optimum angle of attack for the approach condition. 


The oc value must be adjusted as a function of flap position. A practical 
REF 

problem associated with this method involves the accurate measurement of a for 
various flap conditions. 

3. Throttle-Thrust Scaling and Authority Limits 

For simulation purposes, the suggested relationship between the throttle 
limits and thrust is illustrated in Figure 4-4. Note that this figure shows that 
full throttle quadrant can produce greater than 100 percent rated thrust. The 
42-degree throttle limit is selected to provide 100 percent rated thrust under the 
most optimistic conditions favoring increased thrust (cold day , high barometric 
pressure). For nominal conditions, it is suggested that the 42-degree limit cor- 
respond to 90 percent full rated thrust. Hie 13-degree lower limit would then 


72 




Figure 4-4 

Autothrottle Servo Thrust Authority 
Limits (Nominal Day) 




correspond to about 10 percent full rated thrust. Note that for about. the first 
ten degrees of throttle quadrant rotation, the thrust output is about zero; the 
engine output drives only the accessory loads . 


4-1. 


Autothrottle Control Law Summary 
In summary, the autothrottle control law is 

■**. = [ v w + V‘][vHt 1+ rb 


+k A X c + k 2 5 E 


r_v_i r i 

L T 7 S + U L T 8 S + V 


+ 5 


TP 


(4-13) 


where ^ 2 ' 9 may be used in place of ^2 ^E' 

The range of typical control law parameters is summarized on Table 


TABLE 4-1 


AUTOTHROTTLE PARAMETER SUMMARY 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 


4 Rad/Sec 

6 Rad/Sec 

12 Rad/Sec 

High bandwidth 
not important 

L 2 

4 Deg/Sec 

8 Deg/Sec 

12 Deg/Sec 

Throttle rate 
limit - 
(degrees of 
throttle 
quadrant 
rotation) 

L 3 


13 Degrees 
above idle 
stop to 
42 Degrees 


Throttle limits 

T 

C 

2 Sec 

4 Sec 

8 Sec 

Complementary 

filter 

k 

V 

1 .0 Deg 5 t 
per knot 

3 Deg 5 t 

6 Deg 5 t 

Air speed 
error gain 
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TABLE 4-1 (cont) 
AUTOTHROTTLE PARAMETER SUMMARY 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 

k I 

0.025 Deg/Sec 
8^ per knot 

0.1 Deg/Sec 
per knot 

0.25 Deg/Sec 
8^ per knot 

Airspeed 
integral gain 

k A 

2 Deg 8 X per 
knot/sec 

5 Deg 5 X per 
knot/sec 

10 Deg 5 t per 
knot/sec 

Acceleration 

compensation 

gain 

k 2 

1.0 Deg 5^ 
per deg § E 

2 Deg 5 t 
per deg $ E 

5 Deg 8 t 
per deg $ E 

Elevator com- 
pensation 

L 1 

r w 

2.0 knot 

5 Knots 

10 Knots 

Error limit 

30 

100 

150 

If this 
washout 
is used, it 
should be 
disabled for 
large errors 

T 7 

15 

30 

50 

Elevator 

washout 

T 8 

20 

4 

6.0 

Compensation 

filter 

• 

V 

C REF 

0.75 Knots/Sec 

1.0 Knots/Sec 

2.0 Knots/Sec 

Slew rate of 
calibrated 
airspeed 
reference 
(may also be 
computed for 
perfect 
compensation) 
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5 . Autothrottle Control Response Criteria 


a. Airspeed Transients 

Because speed control uses a frequency-weighted combination of pressure 
references and inertial references, the dynamic response to speed changes depends 
upon the actual inertial and wind velocity. In general, the control law is de- 
signed to respond to inertial changes at high frequency and impact pressure 
changes in the lower frequency region. A step increase in calibrated airspeed 
because of a step input of forward wind will cause an aircraft inertial velocity 
deceleration. The inertial terms of the control law will then tend to increase 
throttles while the pressure terms will tend to retard them. Ideally, a wind 
pulse of about 0.5 second duration should cause no throttle response. A step in- 
crease in headwind should produce a response compatible with the final determina- 
tion of the required equilibrium speed. If pitch attitude is constrained, a step 
increase in headwind will cause a positive normal acceleration with a resulting 
increase in flight path angle. With no throttle adjustment the aircraft will ex- 
perience a deceleration until the original airspeed equilibrium is restored. An 
autothrottle system will retard throttles initially, permitting a faster restora- 
tion of the original airspeed. The throttles will then return to their original 
value. If flight path angle with respect to a ground reference is constrained, 
then changes in the equilibrium throttle position will be required. Consider 
descent on a 3 -degree glide path, for example. A step increase in headwind will 
cause a deceleration with respect to the ground track. In order to maintain the 
aircraft's position on the glide path an increased angle of attack is needed. If 
the aircraft is on the stable side of the thrust required curve, it can continue 
on the original glide path but at a lower airspeed. If airspeed is to be main- 
tained, the throttles must be advanced. The autothrottle system should provide 
that advance with a minimum of initial retard activity trying to cope with the 
initial airspeed increase. The complementary filter blend of inertial accelera- 
tion and airspeed helps achieve this desired response. 

The dynamic coupling between the longitudinal stabilization modes and 
the throttle control loops is always significant. Ultimately, the best response 
criteria are those which result in the tightest flight path control. However, 
for initial throttle loop adjustment, the following transient response criteria 

are suggested. 
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Input 

#0.5 second step wind pulse 

• Step 5 knot speed change 
(pitch attitude constrained 
via m avtnmm gain pitch 
autopilot) 


• Apply 5.0 knot speed reference 
change in the form of a 1.0 
knot /second ramp command 


Response Criteria 

• Minimum throttle activity 

• Min imum forward throttle 
transient . 

• Airspeed error reduced. by 90 
percent in 8 seconds . 

• Maximum 10 percent speed 
overshoot (0.5 knots). 

• Airspeed error should never 
exceed 1.5 knots. 

• Speed overshoot should not 
exceed 0.5 knots. 

• Final value within 0.25 knots 
should be attained within 8 
seconds of ramp command com- 
pletion. Error should be 
within 0.5 knots after 4 
seconds of ramp completion. 


B. DIGITAL PROGRAM - AUTOTHROTTLE CONTROL 
1. Control Law Conversion 

The representation of the autothrottle block diagrams and control laws in 
FORTRAN notation is summarized in Figure 4-5 and the FORTRAN namelist is given in 
Table 4-2. Note that in the system configuration shown in Figure 4-5, the pitch 
compensation loop rather than the elevator compensation loop is used. As 
discussed previously these two loops may usually be used interchangeably (with 
appropriate gains) . Also noted on Figure 4-5 is the absence of a feedforward 
throttle predictive input. The predictive input could be used advantageously for 
improving command response but it will have no bearing on the disturbance 

response . 




'■'J 

00 . 



TPCL 


Figure 4-5 

Autothrottle Block Diagram with 
Fortran Notation 







TABLE 4-2 


AUTOTHROTTLE NAMELIST 


Name 

FORTRAN Name 

Definition 

• • 

X ACCEL 

AXB 

Body axis longitudinal accelerometer output. 

V 

VRW 

Aircraft airspeed. 

V ' 

cREF 

VCREFP 

Airspeed reference command. 

• 

V 

cREF 

RL4 (UTH4) 

Slew rate of airspeed reference command. 

S 

RL1 (UTH1) 

Airspeed error magnitude limit. 

V 

TRATE (UTHR) 

Thrust command rate limit. 

S' 

TPMIN 

Thrust command position limit. 

V 

TP MAX 

Thrust command position limit. 

II 

X 

AXC 

Compensated body axis longitudinal accelero- 

c 


meter output . 

e 

c 

THECOM 

Pitch angle command. 

- 

STHET 

Sine of pitch angle command. 

- 

TPCTR 

Trim thrust value in percent. 

- 

TP CL 

Position and rate limited total thrust command 
in percent. 

- 

TPC1 ,2,3,4 

Individual engine thrust commands in 
percent. 

- 

DTMAX 

Percent thrust change per degree throttle 
change . 

- 

ITESTT 

Logic variable for synchronizing autothrottle 
at engagement. 

K 

a 

KA 

Longitudinal acceleration feedback gain. 

K 

V 

KV 

Airspeed error feedback gain. 

K‘K t 
v I 

KINT (KV-KI) 

Integral error feedback gain. 

T 

TAUC 

Time constant for complementary airspeed 

C 


filter. 

K 2 

DP (UDP) 

Feed-forward pitch compensation gain. 
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TABLE 4-2 (cont) 


AUTOTHROTTLE NAMELIST 


Name 

FORTRAN Name 

Definition 

r 7 

TAUTH7 

Pitch compensation filter time constant. 

T 8 

TAUTH8 

Pitch compensation filter time constant. 


DT3 

Subroutine sample time interval. 


2. Program Flow Chart 

The initial condition computations which are performed in the SASIC 
subroutine (See Appendix A) for the autothrottle mode is given in the following 
summary. The flow charts are shown in Figure 4-6. 

Autothrottle I.C. Calculations 

a. Set Engine Thrust Commands to Trim Value 

TPC1 = TPCTR 
TPC2 = TPCTR 
TPC3 = TPCTR 
TPC4 = TPCTR 

b. Convert Gains and Limits to Proper Units 


RL1 

= UTH1 *1 ,688 

RL4 

= UTH4 *1,688 *DT3 

DP 

= UDP/DTMAX 

KA 

= UKA/DTMAX/1 .688 

KV 

= UKV/DTMAX/1 .688 

TRATE 

= UTHR/DTMAX *DT3 


c. Airspeed Error Integrator Gain 

KINT = KI*KV*DT3 

d. Pitch Compensation Filter Difference Equation Coefficients 

Cl = EXP (-DT3/TAUTH7) + EXP (-DT3/TAUTH8) 

C2 = EXP(-DT3*(1 .0/TAUTH7 + 1.0/TAUTH8)) 

D1 = DP/TAUTH8/(1 .0/TAUTH8 - 1.0/TAUTH7) 

*(EXP (-DT3/TAUTH7) - EXP (-DT3/TAUTH8) ) 




/ 

e. Complementary Airspeed Filter Difference Equation Coefficients 

C = EXP (-DT3/TAUC) . 

D = KV*(1 .0-C) 

The autothrottle parameter test values selected for best performance are: 


UDP 

2.5 

Deg Throttle/Deg Theta 

UKA 

4.08 

Deg Throttle/Knot/Sec 

UKV 

6.0 

Deg Throttle/Knot 

KI 

0.05 

Dimensionless 

UTH1 

5.0 

Knots 

UTH4 

1.0 

Knots/Sec 

DTMAX 

36.2 

Deg Throttle/% Max Thrust 

UTHR 

8.0 

Deg/Sec 

TAUTH7 

30.0 

Seconds 

TAUTH8 

2.0 

Seconds 

TAUC 

4.0 

Seconds 


C. SIMULATION TEST RESULTS - AUTOTHROTTLE 

Figures 4-7 and 4-8 summarize the autothrottle transient responses for 
approach and cruise speeds achieved with the control law parameters listed in 
paragraph B.2 above. At both speeds, the 5-knot step headwind response shows a 
small overshoot of about 10 percent maximum. This is within the specified 
criteria but it could have been improved with a higher gain (degrees throttle per 
knot) and switching logic on the integration function. The higher autothrottle 
loop gain would tend to move the response characteristics into the range where 
excessive throttle activity occurs. Logic switching of the integrator involves 
holding the integral gain at zero until an error plus error rate criterion is 
satisfied. With such a technique the integrator would not have started until 
the speed returned to about 240 ft/second in Figure 4-7 or 502 ft/second in 
Figure 4-8. 

The 0.5-second pulse response illustrates the effect of the complementary 
filter in minimizing throttle activity although it obviously has not done the 
job completely. The effectiveness can be demonstrated with a simple calcula- 
tion. The error is 5 knots for 0.5 second. Without the complementary filter 
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[engage] 


I 


X5 

= 0 

X6 

= 0 

DTT1 

= 0 

DTT2 

= 0 

THEC2 

= THECOM 

THEC1 

= THEC2 

VCREFP 

= VRW 

TPCL 

= (TPC1+PTC2+TPC3+TPC4)/4.0 

X8 

= TPCL 

TPC 

«= TPCL 

VCREF 

= VCREFP 

ITESTT 

= 1 




Figure 4-6a 

Auto throttle Flow Chart 
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Figure 4- 6b 

Autothrottle Flow Chart (cont) 















the throttle retard rate would have been 8 degrees per second (neglecting throttle 
servo dynamics). With a thrust gradient of 2.76 percent full thrust command per 
degree of throttle, the commanded thrust change in 0.5 second would have been 11 
percent. As shown in the response, the blend of inertial and air data results in 
a significant attenuation of this throttle activity. A high time constant on the 
complementary filter would result in greater throttle motion attenuation but at 
the expense of excessive stretching of the recovery time to a step wind change. 

The airspeed command responses shown on the third set of traces in Figures 
4-7 and 4-8 do not meet the desired response time criteria but this can easily 
be remedied with the use of a feedforward compensator. The responses obtained 
in the simulation tests depended entirely on the closed loop system. Thus, an 
airspeed error had to exist before throttle movement is initialed. Also, the 
acceleration loops, both in the complementary filter and in the damping loop (KA) , 
oppose the build-up of the corrective throttle. The use of a predictive feed- 
forward compensator would eliminate this problem by providing the throttle change 
needed to yield the 1 .0 knot/second change in airspeed reference plus the A 
throttle needed to bias the acceleration loop (4.08 degrees of throttle per 
knot/second). It is also noted that the use of the switching logic on the 
integral loop, as discussed previously, would help eliminate the small overshoot 
and the long tail in the convergence to the reference airspeed. 
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Figure 4-8 

Autothrottle Tests (Cruise Speed) 
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SECTION V 


VERTICAL FLIGHT PATH GUIDANCE LAWS (NONLANDING) 


A. DESCRIPTION OF CONTROL LAWS 

1 . General 

The vertical guidance functions covered in this section relate to the 
flight path steering commands associated with climb and descent, and the acqui- 
sition and holding of constant altitude. Not covered herein are the final 
approach glide path descents or cruise modes based on airspeed and Mach control. 

Vertical guidance laws are represented as pitch steering commands into 
the autopilot pitch stabilization inner loop. Figure 5-1 shows how the series 
of pitch command outputs of specific modes are transmitted through a synchro- 
nizing function and summed with incremental pitch attitude to create the pitch 
error signal of the pitch stabilization loop. The synchronization function 
causes the last value of pitch command to be held as the initial value for a new 
control mode during mode transitions. (Note that the polarity of the summation 
of pitch command and pitch attitude is positive. See the discussion in Section 
II for the explanation of this polarity.) 

In general, the performance of the vertical guidance modes is dependent 
upon good forward speed control. This may be achieved through the use of the 
automatic throttle control system. 

2. Vertical Speed Control 

a. Control Laws 

Vertical speed control is accomplished with two submodes: Vertical 

Speed Hold and Vertical Speed Command. Both submodes are represented by the 
block diagram on Figure 5-2. The basic vertical speed control law is: 



(5-1) 




0 «* SYNCHRONIZED PITCH ATTITUDE ... 

PITCH ATTITUDE INCREMENT FROM 
VALUE EXISTING AT TIME OF 
AUTOPILOT ENGAGEMENT 


Figure 5-1 

Pitch Guidance Function Engage Synchronization 
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PREDICTIVE 

PITCH 

COMMAND 


Figure 5-2 

Vertical Speed Control Block Diagram 






where 


h 


error 



h 

c 


compensated vertical speed 


(5-2) 


^i T 1 +h B (5-3) 

“ T^S + 1 

h = reference vertical speed 
ref 

U . = inertial vertical acceleration 

h = barometric vertical velocity 
B 

Q = predictive pitch command . . . based on mode sequencing and 
c • 

^ recognition of changes in ^ re £ 

The gain should be made a function of velocity with the following 
characteristic: 


200 

^ ^ * nominal ^ nominal 


for V > V or V > 200 ft/sec 
o 




nominal 


(5-4) 


(5-5) 


for V < 200 ft/sec 

Note that in many practical mechanizations, continuous gain programs on are 
not used. The gain is switched at two or three discrete velocity (or Mach 
number) points to achieve an approximate gain program that permits reasonable 
stability and performance. 
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b. Vertical Speed Reference 


> 


h 

c 


) prior 


For the hold mode, h ref is synchronized to the existing h (using 
to vertical speed hold engage. Thus, 



at t 


t 

o 


(5-6) 


where 

t = time of vertical speed hold engage 
o 

For the vertical speed command mode, the vertical speed reference 
is established as a numerical input to the autopilot via a Mode Control Panel. 

In that case , 




(5-7) 


where [ L] represents a limiting function based on an acceleration constraint and 

h' is the numerical input value of desired vertical speed. The limit L pro- 
ref 

vides the following function: 


ref 


i» re f + [ Slgn ( ^ref “ ^rep] |^max| (t V 

O u 


for .|*ief| > |^ref | ( 5 - g > 

2 

where |h I represents the maneuvering acceleration constraint in feet/sec , 

(t - t ) represents the time from the entry of the new h reference, and h rgf 
o • o 

is the initial value of the h reference. 


For h = h* the ramp function described by equation (5 8) is terminated i 
and the h reference is clamped at the desired value until a new reference is 
entered at the Mode Control Panel. When programmed in a digital autopilot, 
equation (5—8) is implemented as : 

X, ’ X.-, + [ 8lgn C" £ ' 

where At is the computation cycle time. 

2 

A typical value for h mgv = 0.07g*s = 2.25 ft/sec . 


max 


At 


(5-9) 
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c. Predictive Pitch Commands 

To minimize the h error when the h reference is changed, a feed- 
forward (or predictive) pitch command 6 is added to the h control law. A 

p 

good approximation of the required angle-of-attack increment that will yield 

the h associated with the ramp change in h can be computed from the fol- 
max rei 

lowing equations: 

C_ QS Ax + C. QS AS = 7 h (5-10) 

a \ 


C Aot + C 


m. 


a 


ID r 




(5-11) 


Substituting (11) into (10) yields 



( 5 - 12 ) 


Using approximate values for C , C , C , and C can yield a reasonable 

(Ot ^Sj. a 

estimate of the aerodynamic coefficient part of equation (5-12). The main 
variables that will change over the different flight conditions will be weight, 
W, and dynamic pressure, Q, although there is also a variation in the aerody- 
namic terms. The predictive pitch term must include the Aa of equation (5-12), 
plus the flight path angle change associated with the acceleration maneuver. 

For Kef - Vf| > 2 £t/seC 


or 



(5-1 3a) 




dt 


(5-1 3b) 
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(5-13c) 




h 

max 

g 




where 



and the control input is 



(5-14) 


(5-15) 


For 


*ref - ‘‘ref I ^ 2 ft/sec 


0 ' 

c 


a 


The significance of bounding the range of h^^ errors for which the 
predictive feedforward correction is made is the recognition that equations (5-13) 
and (5-14) are approximations. Hence, as the maneuver approaches terminal condi- 
tions, the predictive or open-loop command is removed and the closed-loop control 
takes over the entire task of error correction. 


3. Altitude Preselect and Capture 
a. General 

Preselect altitudes are entered at a Mode Control Panel. When the 
altitude preselect mode is engaged, the autopilot is armed to automatically capture 
the reference altitude. The vertical steering mode that directs the aircraft 
toward the preselected reference altitude may then be selected by the pilot. The 
aircraft is not automatically maneuvered to a descent or climb toward the target. 
The pilot must select the manner in which the preselected altitude will be ap- 
proached. With the altitude preselect mode engaged, the autopilot will automati- 
cally sense the approach to the preselected altitude and initiate an altitude 
| capture maneuver. 
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The computational requirements for the altitude preselect mode involve 
four programs. They are: 

• Compute Capture Initiate Altitude 

• Capture Phase (A) Control Law Computation 
(Vertical Speed Command) 

• Capture Phase (B) Control Law Computations 
(Exponential Flare) 

• Capture Phase (C) Altitude Hold Computation 


b. Altitude Capture Computations and Control Laws 
Phase A - Initial Capture Maneuver 

In order to acquire the reference altitude without exceeding an accel- 
eration limit, the capture maneuver must be started at a distance from the desired 
altitude that is proportional to the square of the vertical speed. For a fixed 
vertical acceleration maneuver representing the maximum acceleration constraint, 

ii , the maneuver must be initiated at an altitude displacement Ah from the refer- 
max 

ence given by: 


| Ah = 


(5-16) 


2h 


max 


If we make provision for the response time required to achieve the 

maximum ii and recognize that the final flare into the reference altitude involves 

accelerations below ii , then Ah should be biased on the high side by assuming an 

max 

actual value of h about 10 percent below the specified maximum. A reasonable 
max 2 

value of h is 0.07g's = 2.25 ft/sec . A 10-percent reduction in this value for 
max 2 

the Ah calculation would result in 0.063g's = 1.84 ft/sec . Thus, a 6000-ft/min 

(100 ft/sec) climb toward a preselected altitude will require capture initiation at 


Ah = 


2 

100 

2(1.84) 


= 2720 feet 


2 

if we seek a 0.07g = 2.25 ft/sec maximum maneuver constraint. 







The first phase of the capture maneuver, designated Phase (A), occurs 
when the altitude error, h g , is defined by 


Ah 

> 

h 



e 


> 


k^h 


(5-17) 


where is a control law gain defined below. Note that this phase is required 
only when very large initial vertical speeds are involved. The pitch command con' 
trol equation for this phase is: 



(5-18) 


This is seen to be the same equation as the vertical speed control law. Note also 

that the h signal is synthesized exactly as in equations (5-2) and (5-8) . The 

error 

Phase A altitude capture mode is actually a vertical speed command mode defined by: 

• Equation (5-18) - Basic Control Law 

• Equations (5-16), (5-17) - Initiation and Duration of Phase A 

• Equation (5-8) - Programming of ^ re f (with desired value = 0) 

• Equations (5-13), (5-14), (5-15) - Computation of Predictive 
Pitch Command 

Figure 5-3 is the block diagram illustrating the altitude capture sequence. 

Phase B - Capture Flare Maneuver 

The transition to Phase B should occur when h^ on Figure 5-3 comes 
out of limiting. At that time, the integral gain a 1 is reduced to zero and the « 
component of the predictive pitch command is decayed to zero through the time 
constant r ^ The transition criteria for Phase B are illustrated in Figure 5-4. 
This is a phase plane of h and h. During Phase A, we follow the constant h tra- 
jectory. If we started the ii maneuver at the precise altitude and maintained the 
precise value of h , the acceleration trajectory would intersect the origin. 

The intersection of the acceleration trajectory with the h + k^ switching boundary 
corresponds to the point at which h^ comes out of limiting (for the ideal) accel- 
eration trajectory. When the state of the aircraft intersects (below the switching 
boundary), it corresponds to the case where the h £ + k^ control law required lower 
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accelerations to reach the final state than E^. The Phase B control law is 
.initiated when the switching boundary is reached. The control law is: 


e 


'2B 



K 


6 ' 

c 

E 

T^s + 1 



(5-19) 


where 6 * is stepped to zero at the time of transition to Phase B while 0 

c p ^ a 

Is retained in 6* , During Phase B, the instantaneous value of h ref continues 
c 

to decay toward zero at the maximum acceleration constraint rate. If altitude 
capture is initiated from relatively low vertical speeds, then Phase A will never 
occur and the initial capture phase will be Phase B. 


An additional criterion must be added to allow switching from the 
Phase A to Phase B control mode. If the acceleration maneuver is too large or 
started too soon, then h = 0 will be reached before the terminal altitude is 
reached. Thus, a minimum vertical velocity criterion must be added to allow initi- 
ation of Phase B control. A reasonable value is 2 ft/sec. Thus, the Phase B 

transition logic is : 


PB = Phase B Mode Engage 

SB = Switching Boundary (Figure 5-4) has been Penetrated 

HDOTMIN = |h| < 2 ft/sec 
HMIN = |hj > 25 ft 

PB = SB * HMIN + HDOTMIN * HMIN 

Note that for |h e | <25 ft. Phase C, altitude hold is initiated as described 
below. 


It is noted that the seemingly complex computations associated with 
altitude capture are primarily the result of the need to constrain accelerations 
to acceptable values for passenger comfort. An alternative to the Phase A con- 
trol law [ equation (5-18)] is to start immediately with Phase B [equation (5-19)1 , 
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but to constrain 0 , the total pitch command, so that the g limits are not 

C 2B 

exceeded. Thus, 

* ~ 2.25 x 57.3 (5-19a) 

c V 

2B 

max 

where 2.25 ft/sec 2 = 0.07g;s is the g constraint. 

Another alternative which permits the elimination of the Phase A 
control law is to continuously compute the instantaneous altitude reference on 
the basis of the integration of the vertical speed reference. Thus, h g in 
equation (5—19) is defined as 

h e " h - h ref 

and 

h ref- h o + j[\ef dt 

where h is the altitude at which capture is initiated. This approach allows for 
compensation of an off-nominal acceleration maneuver where the h e term will act 
as a vernier on the h ref so that the trajectory is guided toward an asymptotic 
intersection of the h ref = 0 with the selected altitude. 


(5-1 9b) 


(5-19c) 


Phase C - Altitude Hold 

At h less than 25 feet, the altitude hold control law is initiated. 

© 

A block diagram for this phase, which involves a transition to a new control law, 
is shown on Figure 5-5. At transition, the pitch command synchronizer holds the 

last value of 8 as the initial condition for all subsequent 8 ^ computations. 

C 2B 3 

The control law is: 


e 

C 3 




(5-20) 


where 8 is a predictive command for flap deployment compensation. 

P f 
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An explicit value for 0 could be computed if the airspeed were con— 

Pf 

stant. However, flap deployment is associated with a planned speed reduction. 

The character of the speed reduction transient involves too many variables to 
achieve a reasonable computed value of 0 .A compromise approach is the fol- 

P f 

lowing simple compensation 


V 


. a. 


f (r 4 s + 1) (r 5 s + 1) 


(5-21) 


The best values of and r, should be obtained from simulations of reasonable 

f 4 

flap deployment programs. The approximate value of is: 


\ 
* 


f C. 


(5-22) 


a 


The washout time constant should be about 4 seconds and the filter T should be 
about 0.25 second. The gain k^ should be a function of (1/V). Thus, 




(5-23) 


nominal 


for (V in ft/sec) 


where k. 


is approximately 0.05 deg/ft. 


nominal 

4. Altitude Hold 

The altitude hold control laws and block diagram are identical to those 
used for Phase C of the altitude capture (Figure 5-5). 

5. Pitch Compensation in Turns 

To minimize altitude loss in turns, a feedforward or predictive pitch 
command is required to provide the nose-up attitude that yields the necessary 
lift increment. The lift increment Al is: 


* . W „ - cos <A 

cos 0 \ cos <P / 


100 


(5-24) 
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Figure 5-5 

Altitude Hold Control Block Diagram 











Ax and the associated AS in 


This increased lift is provided by the nose-up 
accordance with 

AL = C AxQS + c L as e QS 
“ 5 e 


and 


°M 


Aa + 


a 


% 


AS. 


(5-25) 


(5-26) 


Substituting equation (5-26) and (5-24) in (5-25) and solving for 
A a yields 

6a = I " 1 - p °? *1 ~ 

^ cos <(> J QS 

For constant flight path angle 

Aa m Ad 

where 0 is the required pitch change for lift compensation. 

L 

The predictive pitch command should be filtered so that 

\ 

, ■ , *L 

P L ' (V + ’) 

where r, 88 1.0 second, 
o 

6. Stability Considerations 
a. Vertical Speed Control 



(5-27) 


(5-28) 


(5-29) 


Figure 5-6 shows the block diagram of the various control loops and 
transfer functions involved in the vertical speed mode. The vertical speed 
feedback loop is closed as a pitch command to the pitch stabilization inner loop. 
The pitch loop, Gg ( s) modifies the aircraft closed-loop response to an elevator 
input. Since the h error signal is converted to a pitch command (the output of 
G^(s) is a pitch command), then to make the stability analysis loop consistent. 
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GAIN OF VERTICAL SPEED 

0 LOOP CONTROL LAW 


• THE PHUGOID MODE CAN BE ELIMINATED IF A GOOD 
AUTOTHROTTLE LOOP CLOSURE IS ASSUMED. 


Figure 5-6 

Stability Analysis Block Diagram 
Vertical Speed Control 
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b) VERTICAL SPEED CONTROL 
ROOT LOCUS (DISPLACEMENT 
LOOP ONLY) 


Figure 5-7 

Vertical Speed Control Stability Analysis 


\ 

the pitch gain % must be added to the h loop as shown on Figure 5-6. The pitch 
feedback causes a shift in the aircraft's [0/Sgl poles and zeros to the loca- 
tions illustrated on Figure 5-7a. The closure of the h loop causes the root loci 
illustrated on Figure 5-7b. In general, excessive gains of the h loop excite 
relatively short period pitch oscillations. This is in contrast to excessive 
gains of the altitude control loop which cause longer period flight path 

oscillations. 

b. Altitude Control 

The stability analysis block diagram for the altitude control mode 
is shown on Figure 5-8. In this case, the modified aircraft is shown as a 
[0/0 ] transfer function achieved by the closure of the pitch inner loop. ^The 
effect of this loop closure is approximated as a second order response [ H A (s)] 
defined by two real poles. A nominal displacement, integral, and rate control 
law [G (s)l results in a cubic numerator and quadratic denominator. Typical 
ratios h of integral and rate gains yield the control law zeros and poles shown 
on Figure 5-9. As the loop gain 1^ is increased for this fixed ratio of dis- 
placement, integral and rate feedback, the loci of Figure 5-9 are obtained. 

7. Summary of Control Parameters and Performance Criteria 

a. Control Parameters 

The control parameters identified in the control equations given in 
this section are specified in terms of typical minimum, nominal, and maximum 

values in Table 5-1 • 

b. Performance Criteria 

(1) Vertical Speed Control Modes 
(a) Vertical Speed Hold 

Introduce 10 ft/sec offset or initial condition error. 
Response should he restricted to 15-percent overshoot. Final value should be 
attained within 10 seconds, with static accuracy of 2.0 percent, if airspeed is 
maintained via the throttle loop. 
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ALTITUDE CONTROL LAW 


OhW- [f^TT + ? + ,83 J kh 

1 

H A' (s > “ (T A S+ D <T B s+ 1) 

WHERE — AND — ARE THE MODIFIED SHORT PERIOD POLES OF FIGURE 5-7a 
T A t B 

t 

Figure 5-8 

Stability Analysis Block Diagram - Altitude Control 






Figure 5-9 

Altitude Control Root Loci 
(Displacement, Rate and Integral) 
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TABLE 5-1 
PARAMETER SUMMARY 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 

a 1 

0.20 

0.30 

0.50 

Vertical Speed Control 
Integral Gain Ratio 

K 

0.075 

0.15 

0.25 

Vertical Speed Control Pitch 
Command Gain (degrees 6 per 
ft/sec error) 





[ Note gain reduction 
function of velocity per 
equations (5-4) and (5-5)] 

T 1 

2.0 

4.0 

6.0 

Vertical Speed Complementary 
Filter Time Constant 

•• 

h 

max 

1.5 ft/sec 2 

2.25 ft/sec 2 
(0.07g's) 

2 

4.0 ft/sec 

Acceleration Constraint 

T 2 

0.50 

1.0 

1.5 

Predictive Pitch Command 
Filter (seconds) 

k 1 

5 

10 

12 

Altitude Capture Rate to 
Displacement Ratio 

T 3 

0.30 

0.50 

1.0 

Altitude Error Filter 
(seconds) 

a 2 

0.03 

0.05 

0.10 

Altitude Control Integral to 
Displacement Ratio 

a 3 

0.05 

1.0 

4.0 

Altitude Control Rate to 
Displacement Ratio 

"h 

0.03 

0.05 

0.10 

Altitude Control Gain - 
Degrees 6 ^ per Foot Error 



c_ 


(Function of 1/V) 

k F 

— 

S F 

— 

Flap Compensation Gain 

T 4 

1.0 

2.0 

4.0 

Flap Compensation Washout 
(seconds) 

T 5 

0.25 

0.5 

1.0 

Flap Compensation Filter 
(seconds) 
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(b) Vertical Speed Command 

s 

From a vertical speed reference of 10 ft/sec with aircraft 
stabilized to that value, introduce a step reference change of 30 ft/sec to a new 
reference value 40 ft/sec (2400 f t/min) . The normal acceleration should be re- 
stricted to within 10 percent of the maximum limit value (h ). The final con- 

max 

vergence to the new reference vertical speed should be held to an overshoot below 
15 percent of the A h re ^. 

(2) Altitude Capture 

From an initial h of about 100 ft/sec (6000 f t/min), approach the 
reference altitude. Altitude capture initiate, flare to reference altitude, and 
final acquisition of the Preference altitude should be achieved with maximum nor- 
mal acceleration held to within 15 percent of the specified h • The overshoot 

max 

of the reference altitude should be restricted to a maximum of about 50 feet. 
Convergence to within 10 feet of the reference altitude should be / achieved within 
12 seconds after the first overshoot if the capture trajectory overshoots. If 
the capture trajectory undershoots, convergence to within 10 feet of the refer- 
ence altitude should be achieved within 10 seconds after the altitude error 
reached 50 feet. 

(3) Altitude Control 

(a) Transient Response 

With altitude hold mode engaged and altitude error = 0, 
introduce 50 feet altitude error. The corrective response should have its over- 
shoot restricted to 10 feet maximum. The reference value should be reached 
(within 2 feet) in 15 seconds. 

(b) Turn Compensation 

With altitude hold mode engaged and airspeed maintained 
via throttle control, bank the aircraft to a bank angle of 30 degrees at a roll 
rate of 5 degrees/second. The altitude error transient should not exceed 30 
feet and should be reduced to within 10 feet within 10 seconds after the 30-degree 
bank angle is attained. 
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B. DIGITAL PROGRAM - VERTICAL GUIDANCE (NONLANDING) 


1 

. Control Law 

Conversion 

^^hey 

Four subroutines were written to implement the vertical guidanci 

are: 



MEASUR 

Derives compensated vertical speed information 
from blend of inertial and barometric data 


VERTSC 

Vertical Speed Control Law 


ALTHLD 

Altitude Hold Control Law 


HDTCMP 

Vertical Speed Command Processor - Applies 
acceleration constraints to vertical speed 
reference changes 


The FORTRAN namelists for MEASUR, VERTSC, ALTHLD, AND HDTCMP are given 
in Tables 5-2, 5-3, and 5-4, respectively. The block diagrams for these sub- 
routines (in FORTRAN notation) are given in Figures 5-10, 5-11, 5-12, and 5-13, 
respectively. 


TABLE 5-2 

NAMELIST FOR SUBROUTINE MEASUR 


Name 

FORTRAN Name 

Definition 

£ i 

(-) VDDOT 

Vertical acceleration in local vertical 
coordinate axis 


ALTDTB 

Barometric vertical speed 

• 

h 

c 

HDTC 

Compensated vertical speed 

T 1 

TV1 

Filter time constant 

— 

DT3 

Sample time of computation interval 
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TABLE 5-3 

VARIABLE NAMELIST FOR SUBROUTINE VERTSC 


Variable 

FORTRAN. Name 

. . Definition 

h 

ALT 

Aircraft altitude above sea level 

h * 
ref 

ALTREF 

Altitude reference 

• 

h 

c 

HDTC 

Calibrated vertical speed 

h * 
ref 

HDTRFC 

Vertical speed reference command 

V 

VT 

Aircraft ground speed 

e 

c 

THECOM 

Pitch attitude command 

Q 

QBAR 

Dynamic pressure 

w 

WAIT 

Aircraft weight 

C L 

CLALPH 

Aircraft lift curve slope 3c -J^ a 

L a 


■ 

% 

e 

CLDE 

3C T /35 
L e 

c 

m 

a 

CMALPH 

3C /3 a 
m 

C 
m 6 

e 

CMDE 

3C /3 6 
me 

S 

AREA 

Reference area 

k i 

KHC 

Altitude capture displacement to rate ratio 

H 

KHDT 

Vertical speed control pitch command gain 

a l 

Al 

Vertical speed control integral gain ratio 

L 1 

LH1 

Vertical speed command rate limit 

L 2 

LH2 

Altitude error limit for altitude capture 

h 

max 

HDDMAX 

Vertical acceleration limit 

T 2 

TV 2 

Filter time constant 


DT3 

Subroutine sample time 


111 





TABLE 5-4 

VARIABLE NAMELIST FOR SUBROUTINE ALTHLD 


Variable 

FORTRAN Name 

Definition | 

h 

ALT 

Aircraft altitude above sea level 1 

h 

ALTREF 

Reference altitude j 

ref 

I 

• 

h 

c 

HDTC 

Calibrated vertical speed I 

«£ 

DELTF 

Flap deflection angle I 

0 

PHI 

Euler roll attitude angle 

e 

c 

THECOM 

Pitch attitude command j 

Q 

QBAR 

Dynamic pressure I 

w 

Wait 

Aircraft weight 1 

(L 

CLALPH 

• '■ 1 

*1a 



c. 

CLDE 

1 

L 5 



e 


• • ■ ' .. I 

c 

m 

a 

CMALPH 


c 

CMDE 


m 5 


1 

e 


1 

S 

AREA 

Reference area 1 

a 2 

A2 

Altitude control integral to displacement ratio 

a 3 

A3 

Altitude control rate to displacement ratio J 

k h 

KH 

Altitude control gain I 

k F 

KFLAP 

Flap compensation gain 1 

L 3 

LH3 

Altitude error limit 1 

T 2 

TV 3 

Altitude error filter time constant 

r 4 

TF4 

Flap compensator filter time constant j 

T 5 

TF5 

Flap compensator filter time constant j 

t 6 

TC5 

Roll compensator filter time constant j 


R2D 

Radians to degrees conversion factor 


DT3 

Subroutine sample time | 
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H 
-P- - 


LH2 “ KHC • |HDTC| 



MODE CONTROL 

lALXCP CONDITION AND ACTION 

0 VERTICAL SPEED HOLD OR SELECT, ALTITUDE CAPTURE NOT ARMED 

1 VERTICAL SPEED HOLD OR SELECT. ALTITUDE CAPTURE ARMED: 

IF (ALTERR < DELH), IALTCP - 2 _ 

2 VERTICAL SPEED HOLD, SENSE CAPTURE ALTITUDE: IF [(ALTERR < DELHI 
© (ALTERR >LH2)1 IALTCP=3 

3 PHASE A OF CAPTURE: IF [( |HDTC| 12)© (ALTERR > 25) + (AALTER < LH2) ) ] . 
IALTCP = 4 

4 " PHASE B OF CAPTURE: IF (ALTERR < 25) SWITCH TO ALTITUDE HOLD 


Figure 5-11 

Vertical Speed Hold, Select, 

Altitude Capture Block Diagram (Fortran Notation) 


THECOM 

































2. Program Flow Chart 


The initial condition computations which are performed in the SASIC sub- 
routine for the vertical guidance modes are: 

IC Calculations (MEASUR) 

a) Initialize Filter 

HDTCIO = 0 
HDTC = ALTDOT 

b) Filter Difference Equation Coefficients 

CM1 = EXP (-DT3/TV1) 

DM1 = 1.0 - CM1 
IC Calculations (VERTSC) 

a) Predictive Pitch Constants 

KC1 = R2D/ (CLALPH - CLDE * CMALPH/CMDE) / AREA 
WC = WAIT * KC1 * HDDMAX/32.2 

b) Vertical Speed Reference Rate Limit 

LH1 = HDDMAX * DT3 

c) Altitude Limit Term 

THDDMX = 1.8 * HDDMAX 

d) Control Law Gains 

' • KHDV = KHDNOM * 200 
AIT = A1 * DT3 

e) Control Law Difference Equation Coefficients 

CVERT1 = EXP (-DT3/TV2) 

OVERT 1 = 1 - CVERT1 

f) Set Logic 

ITEST1 = 0 
ITEST2 = 0 
ITEST3 = 0 
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IC Calculations (ALTHLD) 


a) Altitude Error Filter Difference Equation Coefficients 

CAH1 = EXP (-DT3/TV3) 

DAH1 = 1.0 - CAH1 

b) Altitude Error Integrator Gain 

DAH2 = A2 * DT3 

c) Roll Compensation Filter Difference Equation Coefficients 

CPL1 = EXP (-DT3/TC5) 

DPL1 - 1.0 - CPL1 

d) Flap Compensation Filter Difference Equation Coefficients 

/ CTF1 = EXP (-DT3/TF4) + EXP (-DT3/TF5) 

CTF2 = EXP [ -DT3 * (1.0/TF4 + 1.0/TF5)] 

DTF = R2D/TF5 * [ EXP (-DT3/TF4) - EXP (-DT3/TF5)] / (1 .0/TF5 
- 1.0/TF4) 

e) Altitude Error Gain 

KHDVA = KHNOM * 200.0 

f) Initialize Synchronization Logic Variable 

ITEST4 = 0 

The flow charts are shown in Figures 5-14 through 5-17. The vertical speed con- 
trol and altitude capture parameter values for optimum performance are: 

KHC = 10.0 

V 

KHDNOM =0.25 deg/ft/sec 
A1 = 0.3 

HDDMAX =2.0 ft/sec 2 
TV 2 =1.0 sec 
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BLOCK DIAGRAM 


HDTRFC 



HDTRF 


II: FLOWCHART 



III: 1C CALCULATION 

LH1 - HDDMAX DT3 

HDDMAX - VERTICAL ACCELERATION LIMIT 

DT3 “ SAMPLE TIME 


Figure 5-13 

Vertical Speed Reference Processor 






MEASUR 




ALTDTB *> ALTDOT 


[a 

jn 


BAROMETRIC ALTITUDE 
(TO BE SIMULATED LATER) 


'] 


HDTCIN - ALTDTB - VDDOT • TV1 
HDTC = CM1 . HDTC + DM1 * HDTCIO 
HDTCIO ** HDTCIN 


[complementary filter] 


T 

^ RETURN ^ 


Figure 5-14 

Measure Subroutine Flow Chart 





f 


ITEST1 

° 1 

ITEST4 

- 0 

XHOUT 

- THECOM/KHDT 

XHINO 

a 0 

THEP10 

■ 0 

THEP1 

■ 0 

THEPS 

* 0 

THEP2 

= 0 

THPOUT 

= 0 

HDTRF 

■ HDTC 

FSIGN 

-1.0 

FSIGNL 

-1.0 

All 

-AIT 

HDTRFP 

» HDTRF 

HDTRFO 

= HDTRFP 





Figure 5- 15a 

VERTSC (Vertical Speed Control) 
Flow Chart 
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|DVH1|>2 


ITEST2 

>0 


THEP1 = THEP1 - THTSTP 
THSTP = 0 
ITEST2 = 0 



THEP1 = THEP1 + THEP2 


I 


THPOUT = CVERT1 • THPOUT + DVERT1 • THEP10 
THEP10 = THEP1 


F / IALTCP 

~ V > o 


ALTERR ” ALT-ALTREF 
AALTER = ABS (ALTERR) 
LH2 ° KHC» ABS (HDTC) 


IALTCP 


THSTP » FSIGN*WC/QBAR 
THEP1 = THEP1 + THSTP 
THEP2 = LH1*FSIGN/VT*R2D 
ITEST2 = 1 



DELH = ABS (HDTRFC *• 2/THDDMX) | | DELH = ABS (HDTC *• 2/THDDMX) 


f/ AALTER \J_ 
\ > DLH / 


AALTER \F_ 
< LH2 / 


PRINT WARNING 


AALTER NT 
> DELH / ~ 


IALTCP = 2 


HDTRFC = SIGN (4.0, HDTRFC) 
IALTCP =3 

XHIN =0 


Figure 5- 15b 
VERTS C Flow Chart (cont) 










PHASE A 


380 


PHASE B 


390 
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Figure 5- 15c 
VERTS C Flow Chart (cont) 
















T 


[engage] 


ITEST4 = 1 
ITEST1 = 0 
ALT0T1 = 0. 

ALTIN = 0. 

ALTOT2 = THECOM/KH 
ALTIN2 =0. 

THETLO = 0. 

THEPF1 =0. 

THEPL = 0. 

DELTF1 = DELTF 
DELTF2 = DELTF 


ALTITUDE 

ERROR 

CALCULATION 


ALTITUDE 

ERROR 

POSITION 

LIMITER 


ALTITUDE 

ERROR 

FILTER 


ALTITUDE 

ERROR 

INTEGRATOR 


t 



\ 


\ 


Figure 5-1 6a 

Altitude Hold (ALT HLD) Flow Chart 
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ROLL COMPENSATOR 
FOR LOSS OF LIFT 
IN A TURN 


n 


FLAP 

COMPENSATION 

[^FILTER 


[PITCH 

COMMAND 


U 


SUMMER 


Figure 5- 16b 
ALTHLD Flow Chart (cont) 
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Figure 5-17 
Intentionally Omitted 








VERTICAL 

FT/SEC 


RESPONSE TO A 4 SEC WIND PULSE 





c. SIMULATION TEST RESULTS 


- VERTICAL GUIDANCE (NONLANDING) 


1 . Vertical Speed Control 

Figure 5-18 illustrates the vertical speed control command response. 

From an initial velocity of 10 ft/sec (climb) the reference is increased to 40 
ft/sec. The response is almost dead-beat with the maximum error restricted to 
about 1.0 ft/sec. Hie acceleration was held to 2.0 ft/sec . A key factor in 
permitting tight tracking of the vertical speed reference is the inclusion of 
an additional lag filter on h ref . As seen on Figure 5-18, h re£ has a first 
order lag of about 1.0 second after its computation based on the h constraint. 
This lag matches the basic vehicle acceleration response capability. If l> ref 
was a pure ramp, it would require infinite acceleration to maintain zero error 

initially. 

The gust response is illustrated in Figure 5-19. During a 10 ft/sec 

(600 ft/min) rate of climb, a 4-second wind pulse (from above) of magnitude 

10 ft/sec is introduced. The response illustrated in Figure 5-19 demonstrates 

that the vertical speed error is arrested by the time one half of the wind 

velocity has been transferred to the aircraft. The overshoot in the recovery 

is caused by the. integration loop. Although this response is satisfactory, the 

overshoot could be largely eliminated by incorporating integral mode switching 

logic. With such an approach, the integrator gain would be held at zero unless 

an fh + kh] criterion was satisfied. \ ,, 

error 


2. Altitude Capture 

The final phase of altitude capture in terms of the vertical speed and 
-acceleration responses is illustrated in Figure 5-20. An altitude overshoot of 
26.9 feet occurs as the 13,000-foot altitude is penetrated. The phase B control 
sequence starts with an altitude rate of about A0 ft/sec. The acceleration 

exceeds the desired 0.07g limit value by a small and acceptable value 
2 

(peak = 3 ft/sec ). 


A complete capture trajectory starting from a vertical speed of about 
4800 ft/min is shown in Figure 5-21. In this case, the alternate control tech- 
nique described in equations (5-19b) and (5-19c) was used. Altitude capture is 
Initiated at an altitude of about 4500 feet. The reference altitude is 6000 feet 
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Page intentionally left blank 






ALTITUDE 



9990 10 20 30 40 BO 60 


ALTITUDE HOLD PARAMETERS: 
LH3 - 100 FT 

TV3 -as SEC 

A2 - 0.04 

A3 - 1.6 

KHNOM » 0.10 DEG/FT 
TC5 -1.0 SEC 

TF4 - 4.0 SEC 

TF5 -5.0 SEC 


Figure 5-22 

Altitude Hold Step Response Tests 
V = 296 kt 
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ALTITUDE HOLD ENGAGED 


ALTITUDE HOLD ENGAGED - NO 
ROLL COMPENSATION 




Figure 5-23 

Performance of Altitude Hold During 
A Turn with a 30° Bank Angle, 

V = 296 Knots 


ALTITUDE HOLD 
PARAMETER: 

A2 “ 0.04 
A3 ° 1.50 
KHNOM “ 0.1 
TV3 - 0.5 
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SECTION VI 


VERTICAL GUIDANCE - FINAL APPROACH AND LANDING 


A. CONTROL LAW DEFINITION 
1 . General 

The vertical flight path, landing guidance laws covered in this section 
relate to the capture or acquisition of an ILS glide path, the precise tracking 
of that glide path to an altitude somewhat below 100 feet, the maintenance of 
that path below 100 feet until a flareout initiate altitude is reached, and the 
flareout maneuver to achieve desired touchdown velocity and position 
objectives. 

A large variety of guidance schemes for automatic landing have been 
studied and applied in actual operational systems. (See Reference. 3, for example.) 
When viewed mathematically, these different schemes can generally be reduced to 
almost identical systems. The motivations for different guidance laws, however, 
are not usually dictated by performance benefits. System design criteria involve 
availability of state variable information, the quality of the available measure- 
ment and considerations related to the redundancy architecture of the total 
system. In this latter category, for example, the choice of sensor may be 
dictated by the cost of providing a triple or quadruple set to meet fail- 
operative and monitoring requirements. Another consideration in system design 
is the selection of a system configuration that minimizes or avoids the introduc- 
tion of new sensors as the flight progresses toward touchdown. Thus, the 
activation of computation sequences requiring new sensors at flareout is 
undesirable because the validity and integrity of the total system, including all 
sensors, should have been established by actual operation prior to initiation of 
flareout. 

The comparative evaluations of the relative merits of different guidance 
configurations, when viewed from such operational considerations, are beyond the 
scope of this report. We are concerned here only with the problem of static and 
dynamic performance in the presence of reasonable disturbances. Nevertheless, the 
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computation sequences defined herein are realistic in that they invo ve f Iters 
that are employed for practical considerations not apparent if one 
1 problem as a mathematical abstraction. The specification of synthesis procedur 
for developing a vertical velocity signal from inertial, barometric and rad 
measurements is an example of this realistic approach. Mathematically all tha^ 
is required is a vertical velocity term, h; but practically, we are def 
computer's computation load associated with generating a usable h signa . 

Although some consideration is given to the problem of state variable 
synthesis (as mentioned above for the case of h), the scope of this report does 
not, in general, cover the problem of state estimation. For example, any 
law using vertical acceleration, fi, can obtain the required function by measur g 
the variable directly or by employing compensators to synthesise the desired term 
from h data. Choosing the correct approach involves analysis o, the measurement 
processes for noise and bandwidth characteristics. This type of analysis is ^ 
general, beyond the scope of this report. For complete generality e gu 
laws are specified in terms of the various state variables as though 
variables are measured directly. This should not preclude the ultimate consider- 
ation of compensators in place of some of these variables. 

Two generic types of systems are described in this report. The first is 
the pitch command landing system. It is consistent with the concept of pitc 
steering to adjust vertical flight paths. All steering commands nay be viewed 
pitch commands into a pitch attitude inner loop. The other type of system 
eliminates the pitch attitude inner loop entirely and replaces t with 
vertical velocity control system augmented with pitch rate feedback for dampin 
and perhaps vertical acceleration for increased tightness (increased bandwidth)^ 
The vertical velocity system is activated at the time of glide slope capture an 
continues to touchdown. It has some practical advantages in the mechanisation o 
redundant systems that must use all flareout sensors in an active manner pr or : to 
f lareout initiation. It also can have some performance advantages in minimising 
flight path disturbances resulting from gusts and wind shears. These performance 
advantages are not clearly established if compared against an optimally compensate 
pitch system. Both systems are described but only the pitch steering systems were 

simulated for this report. 


135 


2. Glide Slope Control Geometry 


Figure 6-1 summarizes the geometry of an aircraft penetrating the ILS 
glide slope to the point where the capture maneuver is initiated. The parameters 
that are used in the control problem are X, the angular deviation from the center 
of the beam; 7 G/S . the nominal beam angle (2.5 degrees for discussions in this 
report); and R, the range from the intercept of the glide slope with the runway. 
Note that the polarity of X is positive when the aircraft is above the center of 
the beam. We assume knowledge of R in the airborne computer. (This assumes an 
area navigation capability or a DME colocated with the localizer transmitter.) 

The aircraft can penetrate the glide path from above or below, but it is 
obviously desirable to have a standard intercept procedure, usually from a 
constant altitude of about 1500 feet. The guidance laws defined in the subsequent 
paragraphs make provision for penetration at off-nominal altitudes and initial 

flight path angles. 

3. Glide Path Capture Phase — Pitch Command System 

When the aircraft penetrates the outer boundaries of the glide slope beam, 
it should have reached a stable situation along the center of the localizer beam, 
and terminal approach speed and flap settings should have been established. 

If a consistent vertical acceleration maneuver is desired for capturing 
the glide slope regardless of initial altitude, then the value of X for starting 
the capture maneuver should vary with range R. If we make the approximation that 
the distance travelled from the start of the capture maneuver at X q to the inter- 
cept of the beam center (X = 0) is equal to the horizontal component of that 

distance. Ax, then 

RX , 

Ax ^ ° ° — = Vt ' (6-1) 

7 G/s 

where 

Ax = distance to beam center 
V = velocity (with respect to ground) 
t = time to beam center 

y = glide slope angle (from aircraft viewpoint) 

G/ S 
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For a constant vertical acceleration equal to h 


c’ 



(6-2) 


where Ay is the change in flight path angle required to align the aircraft with 
the beam center. 


Ay 



(6-3) 


where 7 , the initial flight path angle, is zero if the intercept proceeds from 
a constant altitude initial condition. 


Solving for X q from equations (6-1) and (6-2) yields 



(6-4) 


2 

A reasonable value of h is 1.0 ft/sec . Thus, a constant altitude beam penetra 

c 

tion at 1500 feet of altitude would define a X. as follows: 

o 

V = 238 ft/sec 


Ay = -2.5/57.3 radians 
7 G y g = -2.5/57.3 radians 


X = -0.145 degrees or 
o 


0.145 

0.70 


x 150 = 31.1 


R = 42,700 

•• / 2 
h = -1 ft/sec 

c y 


microamperes of beam signal 


Some additional logic must be added to the glide slope capture initiate 
computation to cope with unusual initial conditions. If this logic is not added, 
a dangerous condition can exist in which the aircraft does not capture the glide 
path until too low an altitude is reached. For example, if the outer boundary 
of the glide slope is penetrated from below while the aircraft is descending in 
a - 2.2 degree flight path angle, the specified capture procedure would not 
initiate glide path capture until the aircraft is near the ground. The aircraft 
was descending almost parallel to the beam. A corresponding condition exists when 
the aircraft is above the beam at a flight path angle of about -2.8 degrees. In 
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these cases the capture phase should be eliminated and the glide path tracking 
mode should be initiated (without integral control). The logic to cope with such 
situations may be summarized as follows : 


Engage Glide Slope Tracking Mode ** G 
Abort Landing => A 


G = (7 < -2°) . 

(X < 0) . 

(h > 600 ft) . 

(1X| <0.751X1 ) 

max 

too steep a 
descent for 
reasonable 
glide slope 
intercept 

below 

beam 

center 

above 600 
feet 

not too far from 
' beam center for 
reasonable 
acquisition 

+ (y > -2.5°) 

. (X > 0) 

. (h > 600 ft) 

. (I XI <0.751X1 

max 

too shallow 
a descent 
for reason- 
able glide 
slope 
intercept 

above 

beam 

center 

above 600 
ft 

not too far from 
beam center 

A =» (h < 600 ft) 

. (G + C) 


' 

below 600 
ft 

neither the 
tracking or 
capture mode 
have been 
initiated 

’✓ 


Figure 6-2 illustrates the computation sequence to initiate glide slope 

capture [equation (6-4)]. It shows that when the capture sequence is initiated 

(the mode logic is defined by c) a predictive pitch command O' is activated. 

C 1 

This predictive command is equivalent to the required change in flight path angle 
since angle of attack and speed are assumed to remain constant through the 
throttle loops. Thus, 


\ 



• • 

h MAX 


Figure 6-2 

Glide Slope Capture Initiation - 
Block Diagram-Pitch Command System 
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where h is the compensated vertical speed existing at initiation of the capture 
co 

mode. Note that h is a compensated or synthesized vertical velocity signal 

c 

generated by blending inertial and barometric data as defined in Section V on 
vertical guidance (non-landing) . 


The predictive pitch command, 0' , is constrained in rate to correspond 

C 1 

to a vertical acceleration limit, Thus, 


£0 


c 


1 


w 

V 


(57.3) deg/sec 


( 6 - 6 ) 


The constrained pitch command is referred to as 0 

C 1 

At the initiation of the capture mode, c, all previously existing pitch 

commands. 0 , are held at their last value prior to mode transition. In addi- 
c 

tion to the pitch command 0 , a vertical speed mode is activated. This mode is 

Cl 

shown in Figure 6-3 as contributing to the pitch command 0 . The vertical speed 

" C 2 
control law that is activated at initiation of c (and remains active through the 

remainder of glide slope flight) is 

% - [k - \ef] h- [** ' W «/s] "i (6 - 7 > 


where 0 is rate constrained to correspond to the acceleration limit by control- 
c o 

z • • • 

ling the rate of change of h ref and h £ is compensated h. The complete control law 

for capture is therefore 


0 = 0 + 0 

C capture °1 C 2 




(6-8) 


4. Glide Slope Tracking — Pitch Command System 

The glide slope tracking phase begins when the beam capture has been 
completed. The following criteria indicate that the beam has been captured: 


• The reference flight path angle, 7 G / g , has been achieved or 
exceeded and X remains negative (below beam center) — 
(overshoot from above or undershoot from below) . 
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• The reference flight path angle, 7 G/S , has not been fully attained, 
but X reached zero or became positive — (overshoot from below or 
undershoot from above) . 


In analog systems where failures in beam detector circuits may occur and we do 
not wish to abort the approach, a timing circuit is often added as a back-up to 
indicate glide slope tracking. Thus, if a nominal glide slope descent vertical 
speed of 10 ft/sec is assumed and a minimum capture maneuver of 1.0 ft/sec is 
commanded, then the desired descent should have been established in 10 seconds 
so that when t - t Q exceeds 10, the glide slope tracking phase G is initiated - 
even if the other detectors had not operated. This timing logic is included 
in the control logic computations, illustrated in Figure 6-3, although its 
utility is doubtful in the digital system. In summary, the glide slope 
tracking phase, G, is initiated in accordance with the following logic equation: 


where 7 


error 


r 7 ' < ol 

. I 

x < o] 

r error j 

1 

J 

Py >0 

• 

[x >ol 

L error 


L J 

[t - t o > 10 

secondsj 


<y - W * » + 2 - 5 °>- 


(6-9) 

( 6 - 10 ) 


At initiation of the glide slope tracking phase, the h loop remains 
closed and the X (beam) control loop is activated as shown in Figure 6-3. The 
control law pitch command is : 



where the gains 1 ^ and kj are programmed as a function of radio altitude (or 
range to touchdown). 
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Since a radio altimeter is needed for the flareout phase, it is desirable 
that it be used during the glide path tracking phase for the gain programming 
function. Figure 6-4 illustrates the gain function as a gain ratio versus radio 
altitude. Since many terrains are sufficiently irregular to question the use of 
such a gain programmer, some compromises are usually incorporated to minimize 
rough terrain problems. For example, the radio altitude sensing logic incorporates 
a unidirectional blocking function that does not allow the h R input to the gain 
computer to increase, but always holds it at the minimum value of altitude previ- 
ously attained. The function g, shown in Figure 6-4, is applied to both k^ and 
k . Note that the integrator gain reduction must be accomplished at the input 
to the integrator and never at the output of the integrator. It is also noted 
that the integrator may be switched on and off on the basis of an integrator con- 
trol logic scheme designed to improve stability. Such a logic scheme would cut 
out the integrator when the error plus error rate exceeds specified thresholds. 

This type of function is not included in this study. 

Although the guidance equations have been specified with a vertical 
speed loop it is apparent that the vertical speed terms can be replaced directly 
with a flight path angle loop (with the application of the appropriate V gain 
adjustment). Historically, h rather than 7 has been used because of such opera- 
tional considerations as: 

• h is available from simple analog interfaces 

• V data is not generally available to the autopilot 

For the digital autopilot these considerations are not as pertinent and, conse- 
quently, the 7 loop is certainly an acceptable alternate. 

5. Flareout Control — Pitch Command System 

As the aircraft tracks the glide slope below 100 feet, dependency upon 
the glide slope signal diminishes. The gain program reduces glide slope gains 
initially to compensate for the fact that an angular beam results In Infinite 
gain at its origin. Below 100 feet, the validity of glide slope signals become 
doubtful and the gain reduction program drives k^ and k^. to zero at an altitude 
of about 60 feet. The flight path control is primarily maintained by the pitch 


Ikk 


GAIN RATIO 



Figure 6-4 

Glide Slope Gain Programming 
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attitude and vertical speed loops below 100 feet. The radio altimeter will 

activate the flareout mode at about 35 feet so that from 60 feet to 35 feet 

(about 2.5 seconds) no glide slope information can be used. Flareout initiates 

a new set of control laws shown in the block diagram in Figure 6-5. The 

0 summer switches from the previous h control loop, holding the last value of 
c 2 

0 prior to mode transition. 

C 2 


Three types of flareout laws are described. They are illustrated in 
Figures 6-5a, b, and c. The first (Figure 6-5a) is the exponential flare 
controller. 


a. Exponential Flare Controller 
The control equations are: 



Vertical speed at touchdown = h^, 
Flare initiate altitude = h 





(6-12) 


h Q = h 1 - f h (6-13) 

Thus if h = 20, f = 2, and h = -10 ft/sec, the flare initiate altitude 
1 • 

would be 40 feet. If h had been -15 ft/sec, the flare initiate altitude 
would be increased to 50 feet. 


The predictive term or feedforward pitch command ideally creates 
the maneuver that satisfies the closed loop control law. The predictive commands 
are nose-up signals of the form: 


or in the time domain. 


e p ( S ) 



r 2 s + 1 



s 


(6-14) 


V *> -*1 


(■ 



dt, 


Where 6 and $ are constants. 
1 2 


(6-15) 
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TO MODE LOGIC 



FLAREOUT PITCH PREDICTOR 


Figure 6-5a 

Exponential Flareout Controller Block Diagram 


NOTE: refer to 6-5b 



The values of 9 and 9 are determined by the aircraft characteristics and 

I fc 

especially the ground effect aerodynamics. Constraints should be applied to 

the ramp term (<?„) so that an excessive command does not develop if the flareout 
^ • 

produces an extended float characteristic. Thus 9 should be brought to zero 

after a specified time duration. Also, 6 may be changed to compensate for off 

nominal velocity conditions. Note that the h term attempts to oppose the flare- 

«• •• 
out maneuver since the h reference is zero while a finite (and changing) h is 

required for the flareout. The predictive pitch commands contain the 

necessary bias program to compensate for the nominal £ signals. - .Thus, the. useful 

h information will be the result of deviations from the nominal trajectory. A 

tight h loop is essential for minimizing flight path disturbances due to 

turbulence. 

The touchdown reference terminal velocity is h^. A value of about 
-1.5 ft/sec is desirable, but values of about -2.0 ft/sec to -2.5 ft/sec are 
being used to minimize the downrange excursion of the aircraft during the flare. 

b. Vertical Velocity Flareout Controller 


The exponential controller is converted to a vertical velocity 
controller by removing the h input from the control law. Thus, its control 
equation is 

. I k ol 

h (6-16) 


"c , = y t)+k i »c- V [ + r\ +k h 


The rationale for its use over the exponential controller is that it can provide 
tighter control to the touchdown vertical velocity, hp. It may do this, however, 
at the expense of fore-aft dispersion on the runway. 


c. Acceleration Flareout Controller 

A terminal controller that always attempts to satisfy the terminal 
vertical velocity requirement, h^, by computing and controlling to the precise 
acceleration that will allow the aircraft to fly an exponential flare to the 
desired terminal condition can be derived as follows : 


An exponential flareout is achieved if 

b + K p - K] m 0 

Differentiating this equation gives 


h+k^h=0orh 



(6-17) 


(6-18) 
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FLAREOUT PITCH PREDICTOR 


Figure 6-5b 

Vertical Velocity Flareout Controller 
(Same as Figure 6-5a but remove h input) 








For any initial h and h state, if we wish to perform an exponential flareout 
that satisfies both the terminal requirements and the initial state, a specific 
value of in equation (6-17) must be satisfied. From equation (6-17), it is 



Substituting equation (6-19) into (6-18) gives 



(6-19) 


( 6 - 20 ) 


Thus if we always flew the h ref defined by equation (6-20) , we will follow an 
exponential path toward h p regardless of our initial state. The computed value 
of h can be used as a reference acceleration and a tight acceleration loop 

closed to try to attain the value of h computed by equation (6-20) . Such a 
controller would be of the form 


% " V° + ^ <K " K ref ) ( } 

The predictive pitch command sets up the nominal flareouts as in the 
previous types of flareout controllers. In this case however, for best results, 
the predictive term should start the flareout maneuver with ^ = 0 until the 
vertical speed is reduced to about one half its original value. Since the 
h computation diverges at h - 0, the computation is constrained by letting 
h*- 2 be the minimum allowable value of h. Also, equation (6-20) becomes 
erroneous if h should become positive. . To prevent this, h ref is set to zero 
when h reaches about -2.5 ft/sec (for h p - -2 ft/sec). Figure 6-5c shows a 
block diagram of this system with the various logic computations needed to control 
the h loop. The loop is closed on logic state M and the h ref is set to zero on 

logic state N. 
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FLARE MODE LOGIC (FI 


M - [h > — 6ft/sec] + IM) 
N » [A > -2.5 ft/sac] + N 



\ 


Figure 6-5c 

Acceleration Flareout Controller 
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d. Comparison of Flareout Controllers 

The three types of flareout controllers described above can provide 
satisfactory performance when properly optimized. It is difficult, however, to 
make a judgment on which is best. All three systems are actually quite similar. 
They depend upon the predictive term for most of the flareout maneuver while the 
closed loop controls act as a vernier. Ultimately, the best system is the one 
that provides the tightest flight path control in the presence of wind and tur- 
bulence disturbances and perhaps measurement errors. The mode of operation of 
all three systems can be described in terms of the h, h phase plane. Figure 6-6 
shows these phase planes for each system. The exponential flareout system 
(Figure 6-6a) always tries to control to a fixed line on the phase plane. A 
large vehicle has difficulty in achieving an h + 5h line. (It typically can 
achieve an h + 2h line,) Also an h + 2h controller will give higher accelera- 
tions than an h + 5h system. The higher gain in h is desired for control tight- 
ness but it does not give the best trajectory. What is more significant, 
however, is that if the aircraft has deviated from the reference h + k^h line, 
it generally does not have the control bandwidth to reacquire that line in the 
remaining time. This is where the acceleration controller [Figure 6-6(c)] should 
have some advantage. It does not try to recover to the original reference line 
but always computes the minimum acceleration needed to complete an exponential 
flareout. This controller, however. Is also restricted by the large aircraft s 
inability to achieve rapid acceleration changes. 


Figure 6-6(b) shows the vertical speed controller’s) phase plane tra- 
jectory. It only tries to achieve the terminal h reference. It should nominally 
reach h at about 8 feet from touchdown. If it flares too high it will tend to 
land at but with a penalty in fore-aft excursion on the runway. If it flares 
too low it may not reach the touchdown reference h. There are techniques for 
adding additional intelligence to this controller so that it can minimize these 
penalties. Likewise, there are techniques which can improve the performance of 
the other two types of flareout controllers . These techniques are beyond the 
scope of this study. Their application would be for situations having touchdown 
position and h requirements more severe than those now being used for transport 
Category III landings. 


152 


ALTITUDE 


a. EXPONENTIAL 
CONTROLLER 


b. VERTICAL VELOCITY 
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Figure 6-6 

Phase Plane for Three Types of 
Flareout Controllers 




6. Glide Slope Capture and Tracking - Vertical Speed Control System 


All of the guidance laws discussed to this point were defined as part of 

a pitch steering system. The basic autopilot mode is pitch attitude hold and 

all guidance laws are defined as pitch attitude commands. An alternate scheme 

uses vertical speed, h, or flight path angle, 7 , in place of pitch attitude as 

the basic autopilot mode. Pitch rate is used to damp this mode but pitch 

attitude feedback is not used. The elimination of the pitch attitude feedback 

allows tighter flight path control in the presence of turbulence and wind 

shears. The reason for this improved capability will be discussed later under 

"Stability Considerations". Consider now the implementation of such a system 

only for the landing modes. Thus for the cruise modes, the autopilot retains 

the pitch attitude steering. At the start of the glide slope capture phase, 

c, it completely eliminates the pitch attitude loop. As shown on Figure 6-7, 

the pitch loop that is used during the cruise steering modes generates a control 

signal, 0 . The landing computer tracks 9 in a synchronizer and holds the 
£ £ 

value existing at mode transition. This value is retained as in the 

o land 

summing stage. Glide slope capture starts at the time X q is reached where X^ 
is computed as in equation (6-4) for the pitch command system. At the instant 
of the c mode transition, control law changes to 


= G 1 (s) q + 5 e + k^ h 


^land 


+ lo- h ; 
error n error 


( 6 - 22 ) 


where G^ (s) is the pitch rate control law as it existed in the pitch control 
system except that a gain increase may now be needed. 

The notations k^ and kjj are used for the vertical speed and acceleration 
gains as in the previous system, but the gains are not the same as for the pitch 
command system. For the capture mode 


h. = h - h 
error c ref 


(6-23) 


and the acceleration constraints are applied to the change in h ^ 
pitch steering system. 
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as in the 
















The vertical acceleration loop is added to improve control tightness 
(although a lead compensator on h error could have theoretically achieved the 
same purpose) . The vertical acceleration error is 


error 


= h - h 


ref 


(6-24) 


where the computation of ii ^ is as follows: 

h . = f-sign h 1 |h 
ref' L & errorj | i 

if |h I > 2 ft/sec 
| errorj 


and 


K = o 

ref 


if |h ^2 ft/sec 

error 


(6-25) 


(6-26) 


The criteria for terminating the capture phase are identical to those of 

equation (6-9). This starts the glide slope tracking or G phase. The 5 

, land 

control law for the G phase is : 


^land 


qG,<s) +8 e +kf,h e + k £ fi 
o 


+ 


\ 

T^S + 1 


X 



(6-27) 


where h £ is defined as in equation (6-24), and and kj are programmed per radio 

altitude as in Figure 6-4. Note that h reference is zero in the ii loop. 


7. Flareout - Vertical Speed Control System 

The flareout control law is switched into the summing point shown on 

Figure 6-7 while the previous h loop is disengaged. The new computations added 

for flareout are identical to those shown for the pitch command system except 

that here a 5_ feedforward or predictive input is used. The flareout control 
. E 

law is , 


\ . * 1 G 1 (s) + 5 e + 

land o 


h c>5*c-V 

r _ 


i + 


k F + Ve + \ 


(6-28) 
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The three types of flareout controllers discussed previously for the pitch 
steering system can be mechanized with equation (6-28) by defining the h^^ and 
controlling the h and h feedbacks. The S , if required, should be similar dynam- 
ically to the 0^ of equations (6-14) and (6-15). 

8. Throttle Control Considerations 

The flareout control parameters, especially the use of predictive feed- 
forward commands, are sensitive to the throttle control program. A consistent 
throttle procedure is needed to assure consistent performance. Thus, all flare- 
out control parameters must be optimized on the basis of an assumed throttle pro- 
gram. The following throttle control procedure is recommended for the flareout: 

• At h ^ 50 feet when the flight path control is maintaining rate 
of descent, start a ramp throttle retard program. The thrust 
retard should be about 5.0 percent maximum rated thrust per 
second. The retard ramp continues to zero thrust, even beyond 
the normal autothrottle 10 percent lower limit. 

• If thrust started at 50 percent maximum rated thrust, the retard 
program should end with about 10 percent maximum rated thrust re- 
maining at touchdown. Longer flare trajectories can have zero 
thrust remaining at touchdown. 

• Flareout starts at about 35 feet. Thus, the throttle retard 
program preceded flare by about 1.5 seconds. 

9. Stability Considerations 

a. Pitch Command System 

Stability aspects of the glide slope control and flareout control are 
identical to those of the nonlanding vertical flight path guidance laws. The ver- 
tical speed and altitude control stability analysis was discussed in Section V, 
"Vertical Flight Path Guidance Laws (Nonlanding)". The vertical speed control 
problem appears in the glide path capture phase. The glide slope tracking phase 
is mathematically identical to the altitude control problem from the viewpoint of 
stability. The reference altitude line is slanted to the glide path angle. The 
stability block diagrams and associated root loci for these two modes were shown 
previously in Figures 5-6, 5-7, 5-8, and 5-9. 
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b. Vertical Speed Control System 

With the attitude loop removed and the stabilization system feedbacks 
are the various derivatives of h, it can be shown that the basic stability prob- 
lem does not differ significantly from the attitude command system. 

An important point to stress is that the h and h systems discussed 
here do not permit the implementation of tighter flight path control systems from 
the standpoint of higher position gains and higher closed loop position control 
frequencies. Their advantage is in the minimization of wind disturbance tran- 
sients because they sense the disturbance sooner and because the new pitch atti 
tude equilibrium in wind shear can be achieved without opposition from the pitch 
attitude loop. From the standpoint of command response or ability to close ini- 
tial condition errors, they have no advantages over the attitude command systems. 
This can be demonstrated by analyzing and comparing the control laws for both 
types of systems. Consider first the system based on the acceleration inner 
loop. The simplified surface command control law will be of the form 

5 (h) = a 1 f hdt + a 2 h + + ^ + 

>— ■ 

where h. is the deviation from the reference path. 

If we substitute 

7 = 0 - a = | 

where 7, 0, a , and h are incremental values from the equilibrium 

6 (h) = f & 1 hdt + a 2 h + a 3 V0 ” a 3 W + ^ + 

neglecting control surface lift effects. 


( 6 - 29 ) 


( 6 - 30 ) 


C QS0£ QC a 
L ct _ a_ 

W/g " W/gS“ 




( 6 - 32 ) 


a = a' + a 

g 


( 6 - 33 ) 


where o’ is the nondisturbed angle of attack and oc g is the equivalent gust angle 
of attack. 
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The undisturbed angle of attack has a well-defined dynamic relation- 
ship to pitch angle which can be approximated for the frequencies of interest as 


where 


r = MV/C t QS 

y L- 


a 


T S 


Td 


Cl ~ 


T/yS + 1 


e = 


T y S + 1 


(M = aircraft mass) 


(6-34) 


(6-35) 


Thus 


for the no-gust condition when a = a’, the control law reduces to 


(h) 




. ; v a « k i - a 3 v) e 


hdt +'a 2 h + a 3 V0 + a ^0 + T s + -j 


(6-36) 


This control law can be compared with the attitude inner loop control law which is 

* 

of the form 

5 = b f tidt + b h + b 0 + b 9 (neglecting the h feedback) (6-37) 

(0) 1 J 2 3 j 

for glide slope control. 

The two control laws differ only by a lagged 0 term. Typically 
a k » a V so that this term is relatively small. It obviously can be added to 
the pitch system so that both control laws would be completely identical. 

I 

The essential difference between the two control laws is in the dis- 
turbance situation where the « g term in equation (6-33) and horizontal gusts are 
the significant contributor to flight path accelerations. The acceleration con- 
trol system provides feedbacks proportional to the actual flight path accelera- 
tions, while the attitude based system requires integration of these accelerations 
until the position errors of equation (6-37) develop. Also changes in the equi- 
librium wind condition (wind shear) require a change in the equilibrium pitch 
attitude. With pitch feedback this necessitates a flight path offset that must be 
corrected by the integrator. \ 
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Figure 6-8 

Root Locus for h Feedback Loop 
showing Effect of Right-half Plane Zero 
Caused by C T 

6 e 
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* •• 

One peculiar stability problem occurs in the h, h type of stabiliza- 
tion system. The preceding analysis neglected the lift due to elevator, C 


This term, when included, results in a right -half plane zero in the h/5^ transfer 
function. The implication of this zero is illustrated in the root locus for an 
h/5 transfer function (Figure 6-8) which is of the form: 

£j 



(6-38) 


A high gain h loop closure rapidly rims into instability as the short-period roots 
bend back into the right-half plane toward the zero, Increasing pitch rate 

gain bends the locus further toward the regions of higher damping, but the pitch 
rate loop has its own stability limitations when it interacts with actuator dy- 
namics not shown in Figure 6-8. Partial compensation is often achieved by moving 
the h sensor forward to the pilot's station where angular acceleration at the 
aircraft's nose cancels some of the initial acceleration reversal resulting from 
a deflection. High gain acceleration loops, however, are generally more dif- 

hi 

ficult to achieve in real aircraft then they are in simulators. 


10. Summary of Control Parameters and Performance Criteria 
a. Control Parameters 

The control parameters identified in the control equations given in 
the previous section are specified in terms of typical minimum, maximum, and 
nominal values in Table 6-1. 


b. Performance Criteria 

(1) Glide Slope Capture 


• Start aircraft at a constant altitude of 1500 feet and at the 
outer fringes of the glide slope beam (0.7 degree). Capture should be accomplished 
with less than 0.1 degree beam overshoot. If an undershoot occurs, it should occur 
within 0.1 degree of beam. That is, if the descent velocity g is attained be- 
fore the beam center is intercepted, it should not occur any further out than 0.1 
degree of beam deflection. 
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TABLE 6-1 
PARAMETER SUMMARY 


Parameter 

Typical Minimum 
Value 

Typical Nominal 
Value 

Typical Maximum 
Value 

Remarks 

7 G/S 

2.5° 

2.5° 


Glide slope angle 

h 

max 

0.025g 

0.05g 


Maximum acceleration con- 
straint for glide slope 
capture V 

r i 


0.05 sec 

0.10 sec 

0.25 sec 

Glide slope filter 

‘'x 


50 

30 

50 

Glide slope displacement gain 
degrees 0 per degree beam (X) 

4 

- Pitch 

0.025 

• 0.04 

0.08 k^ 

Glide slope integral gain 
degrees per second 6 per 
degree beam c 

Gl^de Slope 
k^'- Pitch 

0.05 

0.1 

■ 0.25 

Vertical speed gain degree 

0 per ft/sec 
c 

g " fOO 




Gain reduction program for 
and kj (applicable to both 










/ 

pitch command and vertical 
speed systems) 

*f 


.0.04 

0.05 

1.0 

Flareout gain 

Flareout 
(kj.k^) - Pitch 

0.10 

0.20 

0.40 

Vertical speed gain degree 

0 per ft/sec foT flareout 
c 



0.1 

0.25 

1.0 

Vertical acceleration gain 
for flareout - degree 6 

• 0 e 

per ft/sec 

k 2 

- Pitch 

0.15 

0.25 

0.4 

Flareout integral gain ratio 

T 2 


0.5 

1.0 

2.0 

. 

Pitch flareout predictor 
time constant 

*1 


1.0 

2.0 

2.5 

Displacement component of 
predictive pitch term 

*2 

1 

0.30 deg/sec for 
1st. 5 sec, 0. 10 
deg/sec for 2nd 
5 sec, or until 
touchdown 

0.40 deg/sec for 
1st 5 sec, 0.10 
deg/sec for 2nd 
5 sec, or until 
touchdown 

0.50' deg/sec for 
1st 5 sec, 0.15 
deg/sec for 2nd 
5 sec, or until 
touchdown 

Pitch rate predictive com- 
mand .... (requires maximum 
constraint on integral 
output) 

“if 

Vert 

1.0 

1.8 

2.4 

Degrees 5 £ per ft/sec h error 

h 

h 

Speed 
. Guide 
System 

100 

0.06 

150 

0.20 

200 

0.40 

Degrees S £ per degree beam 

Degrees per second 5 per 
degree beam E 



0.8 

1.25 

2.5 

Degree per ft/sec^ h error 
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• When the level flight capture is optimized, perform a 
capture with a 2.0-degree initial descent angle intercepting the beam outer 
boundary (0.7 degree) at 1500 feet altitude. The logic procedure defined 
in Appendix A should produce a successful capture with criteria the same as 
above . 

• Perform a steep-angle capture by initializing the aircraft 
with a 4.0-degree descent flight path angle at the outer boundary of the beam 
at an altitude of 1800 feet. Overshoot and undershoot criteria are similar 
to those for level flight, except the limit values of A are increased from 0.1 
degree to 0.2 degree. 

To optimize performance on beam capture, adjust for tightest 

o 

h loop gains consistent with stability and then insert the feedforward compen- 
sation to minimize errors. 

(2) Glide Slope Tracking 

• Steady-st;te tracking errors should be reduced to X = zero 
± 0.02 degree by the time an altitude of 500 feet is reached. 

• With 5 feet per second vertical gust pulses of 2-second 
duration, recovery to zero ± 0.02 degree beam error should occur within 12 
seconds following gust removal. This transient should be inserted at h = 1000 
feet, 600 feet, and 300 feet. Damping of the flight path transient or associ- 
ated inner-loop modes should exceed 0.4. Damping of about 0.7 in these 
responses is desirable. 

j 

• At an altitude of 400 feet, introduce a wind shear of 4 knots 
per 100 feet. The aircraft deviation from beam center at 100 feet shall not ex- 
ceed 0.12 degree. 

To optimize glide slope tracking performance, use the highest 
gains of k^ and k^ consistent with stability. The gain reduction program should 
ensure that instability does not occur below altitudes of about 200 feet. Mini- 
mum damping of flight path oscillations should be about 0.4. 
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(3) Flareout 

• With nominal conditions (zero wind or gusts) and the aircraft 
^.igned with the center of the glide path beam, let the aircraft fly into the 
ground without any flareout. Observe the change in h due to ground effect. Use 
the recommended throttle retard program. 

• Add the predictive pitch program optimized to yield a touch- 
down velocity of about 2 feet per second and a minimum runway excursion beyond 
glide slope runway intercept. Flare should start at about 35 feet to minimize 
runway excursion. 

• Add closed loop flareout law with tightest gains achievable 
without causing instability or oscillatory responses. 

• Test the nominal flareout system under conditions of wind- 
shear (4 knots per 100 feet), nominal turbulence, and combinations of head wind 
and tail wind. Successful landings are those that have touchdown vertical 
speeds of less than 4 feet per second and runway dispersions of -300 feet to 
+1200 feet of the ILS reference point. 


164 



B. DIGITAL PROGRAM - VERTICAL GUIDANCE (LANDING) 


1. Control Lav Conversion 

Figure 6-9 is the block diagram (in FORTRAN notation) oi the glide path 
capture and tracking part of the vertical landing guidance laws. (Table 6-2 
summarizes the FORTRAN designations.) The sequencing of functions in acc °^ a 
with the mode switching logic is summarized in the table on Figure 6- . 
ponential flareout block diagram is illustrated in Figure 6-10. Note that, as 
mechanized here, the acceleration loop has not been included. The accelera on 
loop with the h reference set to zero could tighten the system response to gu 
disturbances. Also note that the predictive pitch command is inserted at are 
initiation which is 1.5 seconds before the closed loop h + h system is activated. 
This is one possible variation of this system although adequate performance can 
also be obtained without this delay. The rationale behind delaying initiation of 
closed loop control is the recognition that satisfaction of the h + h contro aw 
requires a continuous h with the largest h at the start of flare. Since t ere 
a lag in the attainment of the initial h, the closed loop system tends to overcom- 
pensate. This problem can be handled with a properly shaped predictive command 
„ by delaying initiation of closed loop control as in the particular implementa- 

tion shown in Figure 6-10. 

An alternate flare control system based on the computed acceleration 
needed to achieve the specified final value of h is shown in the FORTRAN block 
diagram in Figure 6-11. As in the case of the exponential flare control system, 
closed loop control is delayed for 1.5 seconds. Also, the control law changes to 
a vertical speed control if the vertical speed is arrested too rapidly. 

It is noted that the specific Implementations for both flare controllers 
illustrated in Figures 6-10 and 6-11 have weaknesses which will compromise their 
performance in winds and turbulence. The brief delay before closed loop control 
is initiated represents loss of closed loop control for part of the flare maneu- 
ver. For the nominal, no-wind case, performance can be made perfect. With steady 
winds, shears and gusts, the time for correction is cut short so that the system 
should suffer somewhat from the weaknesses of open loop systems sensitivity to 
disturbances. Another feature not incorporated in the Implementation shown in 
Figures 6-10 and 6-11 is the variable flare initiation altitude based on satis- 
fying an altitude plus altitude rate criterion. The variation in flare initiation 
altitude helps to minimize dispersions due to headwind and tailwind variations. 


VERTICAL LANDING CONTROL, PITCH COMMAND SYSTEM SUBROUTINE Vt.ANDI 


KGS : 


( KGSNOM • (. 
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1 
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3 
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5 
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ALCULATE SAFETY OF CAPTURE (SEE FLOW CHART). IF SAFE, CALCULATE AND 
EST FOR CAPTURE ANGLE. WHEN REACHED, SWITCH TO 2. 

ALCULATE: THEC1P =THECOM2 GSD * HDTC/VT • R2D, VHDTRFC = VT * GS 
THEC1 = HDMAXL/VT * DT3 • R2D, SWITCH TO 3. 

EGIN CAPTURE MANEUVER. TEST FOR CAPTURE COMPLETION. WHEN CAPTURE IS 
OMPLETE, SWITCH TO 4 (SEE FLOW CHART). 

YNCHRONIZE GS TRACKING LAW: SET INTEGRATOR, THCINT = 0 AND BIAS THEC10, 
WITCH TO 5. 

RACK GS, TEST FOR BEGINNING OF THRUST REDUCTION PROGRAM (ALT < 50). WHEN 
,LT 50, SWITCH TO 6. 

nwTiNMP TRACKING START THRUST REDUCTION (SEE FLOWCHART). TEST FOR 
LARE INITIATION. WHEN FLARE INITIATION ALTITUDE IS REACHED, SWITCH TO 7. 

NIITIALIZE FLARE, SWITCH TO 8. 

LARE LAW. 
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1 Figure 6-9 

Vertical Landing Guidance-Glide Slope 
Control (Fortran Notation) 













1. begin FLARE: AFTER 1.5 SECONDS; SWITCH TO © 

2. ENGAGE FEEDBACK CONTROL LAW AFTER SPECIFIED DELAY TIME 

3. HALT FEED FORWARD PITCH RAMP COMMAND 


Figure 6-10 

Standard Flare Control System 
Fortran Notation 








■ THECOM 


1. BEGIN FLARE: START FLARE TIMER (TRAMP). WHEN TRAMP - 1.5 SEC TO 5.0 SEC, SWITCH TO 2 

2. ENGAGE FEEDBACK CONTROL LAW AFTER SPECIFIED DELAY TIME 

3. HALT THE FEED FORWARD PITCH RAMP. WHEN HDTC - -2.5 FT/SEC, SWITCH TO 4 

4. CONSTANT VERTICAL SPEED PORTION OF FLARE, HDTF - -2.0 FT/SEC 


Figure 6-11 

Flare Control System: Acceleration Feedback Control Law 

Fortran Notation 





TABLE 6-2 

VERTICAL LANDING GUIDANCE FORTRAN NAMELIST 


Variable 

FORTRAN Name 

Definition 

h R 

ALTRAD 

Aircraft radar altitude above ground 

V 

VT 

Aircraft ground speed 

A 

EPSGS 

Glide slope displacement error angle 

• 

h 

c 

HDTC 

Compensated vertical speed 

•• 

h 

VDDOT 

Vertical acceleration 

7 GS 

THETGS 

Glide slope beam angle above horizon 

X 

0 

EPSGSO 

Glide slope error angle for capture initiation 

O' 

C 1 

THEC1P 

Predictive step command for glide slope capture 

^ref 

VHDTRC 

Vertical speed reference command for capture 

— 

TIMCP 

Glide slope capture duration time 

— 

DTC 

Thrust reduction increment 

— 

LTHEC1 

Rate limit on O' 
c 

H 

KGHD 

Glide slope capture, tracking vertical speed error 
gain 


KGSNOM 

Glide slope displacement error gain, nominal 

h'K 

DGS (KGSINT) 

Ratio of glide slope integral to displacement gain 

T i 

TAUG1 

Glide slope displacement error filter time constant 

V T 5 

TVL1F4, TVL1F5 

Flap compensation filter time constants 

5 f 

DELTF 

t 

Flap angle 

e , 

THECFS 

Predictive step pitch command for flare 

* 2 

THCFF, THCFR 

Predictive ramp pitch command, fast and slow rates, 
respectively 


HDTF 

Desired vertical speed at touchdown 

K 

KFHD 

Vertical speed error gain for flare 


KFLARE 

Total flare gain 
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TABLE 6-2 (cont) 

VERTICAL LANDING GUIDANCE FORTRAN NAMELIST 


Variable 

FORTRAN Name 

Definition 

k 2 

KF2 

Flare integral gain 

t 3 

TAUF3 

Predictive pitch step command filter time constant 


TRAMP 

Predictive ramp pitch command timer 

— 

DELAY 

Time delay for engaging feedback control law in the 
flare 


KFHDDT 

Vertical acceleration gain 

T 1 

TAUF1 

Vertical acceleration filter time constant 

THEC10 

Pitch command synchronization term 

— 

HDINT 

Pitch command synchronization term 


2. Flow Charts 


The flow charts for the landing guidance program (VLAND1 ) are shown in 
Figures 6-12a through 6-12j. The initial condition calculations are summarized 
on Table 6-3. 

TABLE 6-3 

SUMMARY OF VLANDl 
INITIAL CONDITION CALCULATIONS 


LHDT1 = HDMAXL*DT3 

CTF1 = EXP (-DT3/ TVL 1 F4 ) +EXP ( -DT 3 /TVL2F 5 ) 

CTF2 = EXP(-DT3*(1/TVL1F4+1/TVL1F5)) 

DTF = R2D/TVL1F5*(EXP(-DT3/TVL1F4)-EXP(-DT3/TVL1F5))/ (1/TVL1F5-1/TVL1F4) 
CGI = EXP (-DT3/TAUG 1 ) 

DG2 = 1-CG1 

CF1 = EXP(-DT3/TAUF3) 

DF2 = 1-CF1 

CF2 = EXP(-DT3/TAUF1) 

DF3 = 1-CF2 
DGS = KGS1NT*DT3 
THECFF = THCFF*DT3 
THECFR = THCFR*DT3 
GSD = THETGS 
GS = GSD*D2R 
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ITEST6 ** 1 

I LAND - ISTART 

THCINT “ 0 

DTC “ 0.05 • DT3 

EPSGS1 - 0 

EPSGSF - 0 

HDDTF “ 0 

EPSGS2 '» 0 

THECF1 -0 

VDDOT1 «* 0 

THEPF1 -0 

THEPF2 *0 

DELTF1 » DELTF 

DELTF2 “ DELTF 

VHDTRC ■» VT • GS 

LTHEC1 ■ HDMAXL/VT * DT3 • R2D 

THEC1 - THECOM 

VHDTRF - HDTC 


m 
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Figure 6—1 2d 
VLAND1 Flow Chart 
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PROCESS 
VERTICAL SPEED 
REFERENCE 


'] 


FLAP 

COMPENSATION 

FILTER 


[ PITCH COMMAND"! 
SUMMATION J 


CALL HDTCMP (VHDTRC. VHDTRF. LHD1, DVH1) 


, 

' 

THECF - CTF1 • THEPF1 — CTF2 • THEPF2+ DTF + (DELTF1 — DELTF2) 

THEPF2 - THEPF1 

THEPF1 - THECF 

DELTF2 “ DELTF1 

DELTF1 *> DELTF 


f 

THECOM - THEC1 + KGHD 

. (HDTC - VHDTRF) + THECF 




Figure 6-12e 
VLAND 1 Flow Chart 








GLIDE SLOPE TRACKING 
GAIN REDUCTION 


ALTRAD 


KGS = KGSNOM « 0.000812 • (ALTRAD + 22.21) 


ALTRAD 


KGS = 0.0025 « (ALTRAD -60) • KGSNOM 


ALTRAD > 60 , 



GS ERROR 
INTEGRATOR 


GS ERROR 
PITCH COMMAND 


Figure 6— 1 2f 
VLAND1 Flow Chart 
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Figure 6-12g 
VLAND1 Flow Chart 
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VLAND1 Flow Chart 
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- Figure 6-12i 
VLAND1 Flow Chart 
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c. SIMULATION TEST RESULTS - VERTICAL GUIDANCE (LANDING) 

1. Clide Slope Capture 

Penetration and capture of a 2 . 5 -degree glide slope Iron a constant ^ 
altitude (1500 feet) and from an initial diving trajectory (flight pat * 

-4*, is demonstrated in Figure 6-13. The two cases illustrated cover a capt 

of the glide Slope from helow and from ahove. In hoth Ti “In- 

to a maximum of about 0.015 degree of glide slope beam deviation. This ov 
shoot would be barely discernible on the pilot's display. The acceleration co 
attaint was only 0.025, which indicates that the glide slope capture — 
can be achieved without detectable accelerations. The charac eristics of th 
overshoot are typical of the integral loop. In the trajectories ill. «a cd 
glide Slope tracking is initiated at offsets near 0.05 degree and the integr 
control law starts at that time. Delaying the start of integral c-ro^e r 
on the basis of timing logic or error plus error rate criteria would eliminat 

most of this small overshoot problem. 

2. Glide Slope Tracking 

Response to 5-knot wind pulses applied at 1000 feet, 600 feet and 300 
feet are Illustrated in Figure 6-14. Although the gust vertical velocity^ 
about 8-1/2 feet per second, the disturbances are small re at ve o 
tion angle. The overshoot in the responses is partially the result o the 
integration term in the control law and as discussed in previous sections this 
characteristic can be remedied with additional switching logic on Che integ 
r What appears to be an offset tendency in Figure 6-14 is actually^ con- 
sequence of the converging beam. (A constant beam displacement corresponds 
a convergence toward the center of the beam in distance units.) This is apparen 
InTe le of the wind disturbance at 300 feet of altitude. In this case, the 
final part of the transient occurs as the beam convergence becomes very pro- 
nounced (near .00 feet of altitude). As shown by the locus of points correspond- 
ing to a 1.0-foot offset above the beam centerline, a 1.0-foot displacemen 
begins to look like a rapid divergence. The 100-foot decision height is reached 
at about t = 18 seconds. The glide slope error is about 0.02 degree at that 
point. The Category II and III glide slope window is 35 microamps or about 
0.175 degree of beam. Thus the responses shown on Figure 6-14 are all well 
within Category II/ III window requirements. 
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DEGREES X«e DEGREES Aq S DEGREES 


RESPONSE OF GLIDE SLOPE TRACKING TO A 2 SECOND VERTICAL WIND PULSE OF 5 KNOTS 



t (SEC) 




GLIDE SLOPE TRACKING CONTROL 
LAW PARAMETERS: 

KGSNOM » 50 
KGSINT >0.08 
KGHD - 0.2 DEG/FT/SEC 
TAUG1 -0.1 SEC 


Figure 6-14 

Glide Slope Tracking Transient Responses 






3. Flareout Response 


The nominal no-wind flareout characteristic for the exponential flare 
controller (standard flare) and the acceleration flare controller are illus- 
trated as h vs h phase plane trajectories on Figure 6-15. The effectiveness of 
a flareout controller is determined by performance in a disturbance environment. 
Statistical data is needed to measure this performance. Histograms of the 
standard flare controller performance in turbulence (plus strong headwind) are 
given in Figure 6-16. One hundred runs were used to obtain this data. In gen- 
eral, the performance meets the FAA specified criteria in regard to touchdown 
dispersion (Reference 4). There are no standard criteria for maximum values of 
touchdown vertical velocity since this limit varies with individual aircraft. 

In Figure 6-16, it is seen that 78 percent of the landings had touchdown verti- 
cal velocities below 4 feet per second. No landings exceeded 6 feet per second. 
This would generally be considered satisfactory performance. Performance with 
the acceleration controller was not as successful, with touchdown h in turbu- 
lence tending to run about 30 percent higher than for the standard flare. A 
statistically significant sample of runs for the acceleration controller was not 
obtained. It is noted that a - 6 feet per second vertical gusts represents 
fairly severe turbulence. 
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VERTICAL SPEED (FT/SEC) 


TOUCHDOWN DISTANCE FROM GS BEAM GROUND 
INTERCEPT 

(?) x - 504 FT (ACCELERATION FEEDBACK FLARE) 

( 5 ) x - 594 FT (STANDARD FLARE) (EXPONENTIAL) 


FLARE LAW PARAMETERS: 


© 


ACCELERATION FEEDBACK FLARE 

THECFS -2.18° 

THCFF - 0.36 DEG/SEC 
THCFR - 0 
KFHD - 0.20 
KFHDDT - 0.25 
TAUF3 - 0.5 


(5) STANDARD FLARE 

THECFS - 2.18° 

THCFF “ 0.36 DEG/SEC 
THCFR - 0 
KFLARE - 0.40 
KFHD -0.15 
TAUF3 - 0.5 


Figure 6-15 

Nominal (No Wind) Flare Performance 
(Phase Plane Trajectories) 
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VERTICAL SPEED AT RUNWAY CENTERLINE - FT 

TOUCHDOWN — FT/SEC 



FLARE LAW PARAMETERS: 
THECFS - 2.18° 

THCFF ** 0.38 DEG/SEC 
THCFR - 0 


KFLARE - 0.40 
KFHD - 2.0 
KF2 - 0.15 
TAUF3 - 0.5 


Figure 6-16 

Standard Flare Touchdown Performance in 
Turbulence (Mean Headwind - -30 ft/sec, 
a Gust = 6.0 ft/sec) 
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SECTION VII 


LATERAL GUIDANCE 


A. CONTROL LAW DEFINITION 

1. Heading Control Modes 
a. General 

There are three heading control submodes . From the standpoint of 
control dynamics, they are identical. They differ because of operational pro- 
cedures associated with reference data entry, mode selection, and computation 
requirements. The three submodes are: 

• Heading hold 

• Heading preselect 

• Heading command 

Heading hold is the basic lateral steering mode that is engaged automatically 
when other steering modes are not selected and manually commanded bank angles 
fall below a specified threshold value (about 5 degrees). The heading hold 
mode provides an automatic wing leveling capability. Heading preselect allows 
a reference heading to be entered while some other steering mode (including 
heading hold) is engaged. Heading preselect is not actually a control mode, 
but it is the initializing stage of the heading command mode (also referred to 
as heading select). A desired heading is entered, but control is not initi- 
ated until the heading command mode is engaged. 

Provision is also made for a bank command mode. Various methods 
exist for manual bank angle insertion. Transport aircraft autopilots generally 
have turn knobs or control wheel steering sensing devices for this purpose. 

The computation and logic requirements associated with the manual bank modes 
are not covered in this report except for basic logic provisions that allow 
disengagement and synchronization of heading error signals when a manual bank 
command is received and reengagement of heading hold when bank commands are 
removed. Note that automatic turn coordination is implicit in the lateral/ 
directional stabilization system. That is, all steering commands are executed 
to provide turn coordination by virtue of the lateral stabilization control 
laws defined in Section III. 
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Figure 7-1 

Heading Control Mode Block Diagram 







b. Heading Hold Control 
(1) Control Law 

The basic heading control block diagram is shown in Figure 7-1. 


The heading control law used for all heading submodes is 




r k j_ 

L t a 8 + 1 


* 


ERROR 


(7-1) 


where the roll error is 


<t> E - (0 " $ c > " $ + 


r 

Lv + \ 




ERROR 


The definition of * m£)R Mfe) Spends upon the mode logic. 

mode, is defined as follows. . - - — 

E 

# “ W' " 


(7-2) 

For the heading hold 
(7-3) 


where 


kEF y o 


\p is the heading that exists at 


(7-3a) 


the instant the heading hold mode is engaged. 


The gain k^ should be made a 
turn rate capability with velocity. Thus, 


k = a — 
k * 1 v„ 


function of velocity to compensate for the change in 


(7-4) 


where 


V 88 200 ft/sec 
o 

The limits L 1 and L represent roll and roll rate command limits, respectively 
A typical bank limit (L^ is about 30 degrees, and a roll rate limit is about 

5 degrees per second. 


.See discussion of in following paragraph on heading command mode for method 
of resolving the zero/360° ambiguity. 


189 


1 



(2) Heading Hold Logic 

The following symbols are used to represent mode logic states : 

HH = Heading Hold Engaged 

BC = Manual Bank Command Mode Engaged 

BT = Bank Threshold Exceeded ( I $1 > 5°) 

HC = Heading Command Mode Selected 
( ) = Negation 

HH = [ (BC) + (HC)] • BT + [ (HH) * (BC) • (HC)] (7-5) 


c. Heading Command 

The heading command computation provides the heading reference stor- 
age and synchronization. Heading synchronization is performed in accordance with 
equation (7-3a) . The elimination of the zero or 360-degree ambiguity in the 
computation of heading error may be accomplished as follows: 


<* - '"ref ) -'K 

For -180° < \}j ' <+180°, that is I I < 180‘ 


(7-6) 

(7-7) 




For > +180° 

■E* 


(7-8) 


For < -180* 

Ej 


= 360 ‘ 


\ - r E + 360- 


For = ±180° 

£i 


or 


(7-9) 


<L, - +180°, if <{> < 0 (left bank) 
£ 


^ = -180°, if <P > 0 (right bank) 

£ 


To illustrate the above, consider the case where ^ ^ = 30° and 

the aircraft heading = 340° 

K = * - *REF " 340 ° " 30 ° ” 31 °° 
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since 


ifg > 180° 

^ = 310° - 360° - -50° 

E 

Thus, to reach the heading reference, the aircraft will command a bank angle in 
accordance with equation (7-1). For the case illustrated, <t> c will be 




r .s + 1 
A 


(-50) 


50 k. 

+ rrrr r 

A 


which is a right bank angle command. 

d. Synchronization - Data Hold 

The block identified as the synchronizer on Figure 7-1 provides the 
necessary mode engage and disengage smoothing. At the time of any mode transi- 
tion, the existing value of bank command (0 £ ) resulting from a previously com- 
puted error signal is decayed to zero as follows: 

t r s 

d) = (k T 

C R C oV + 1 

where $ is the residual roll command that is to be reduced to zero, 

C R 

is the value of roll command at the instant of mode transition. 


(7-10) 
and 0_ 


e. Stability Considerations 

Heading control by rolling the aircraft to correct heading errors 
involves a single integration loop (associated with turning kinematics) plus the 
roll stabilization dynamics and filter lag (Figure 7-2a) . Note that the lag 
filter, t^, on the heading error is used primarily to decouple any residual dutch 
roll oscillatory tendencies not fully compensated by the yaw damper and roll 
stabilization loops. This filter also helps prevent excessive rolling as a re- 
sult of turbulence. The penalty for this filter is a significant lag (1 to 2 
seconds) and a compromise in attainable gain. The typical root locus is shown 
in Figure 7-2b. 
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The heading control loop gain should be made proportional to velocity 
to compensate for the velocity term in the aircraft’s turning kinematics (*/*>• 

In some applications, heading error integral is used to improve heading control 
static accuracy. Asymmetric thrust would result in a heading standoff (with a 
resultant wing-down condition). The use of integral control could prevent the 
heading standoff error. It is not recommended because the heading mode does not 
generally have stringent accuracy requirements. A good guidance system readily 
recomputes the heading correction needed to correct a flight path error caused by 
a condition such as asymmetric thrust* 

An important factor in the analysis of heading control stability is 
the definition of the heading angle If * is the euler angle determined by a 
yaw, pitch, roll sequence of rotations from a local vertical coordinate frame, 
then the heading control laws given in equations (7-1) and (7-3) are not ade- 
quately represented by the stability analysis (Figure 7-2) for large angles of 
attack and large bank angles. The problem results from the fact that the turning 
kinematics are only an approximation of the azimuth change experienced by the air- 
craft’s X axis. To illustrate the problem without the required derivation of the 
geometrical relationships, consider the hypothetical case of an aircraft in hori- 
zontal flight with a 90-degree angle of attack (pitch angle = +90 degrees). Now 
perform a zero sideslip bank about the velocity vector. A bank about the velocity 
vector is all body axis yaw rate and zero body axis roll rate. Let the roll angle 
change about the velocity vector be 90 degrees. The initial result is that the 
angle of attack remains 90 degrees (no change in velocity vector) but the azimuth 
angle has changed 90 degrees. This corresponds to a case where the azimuth rate 
is the rate of roll about the velocity vector rather than g/V tan <t> , the relation- 
ship defined by the turning kinematics. 

Thus if the control laws specified in this report are implemented 
using euler angle <//, then serious stability problems can occur for high angle-of- 
attack flight conditions (at high velocities). To overcome this problem, the 
angle can be interpreted as inertial velocity vector heading defined as: 

* ‘ tan '’ V n<.rth /V east 
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£. Control Parameter Summary 

The control parameters Identified In the control equations given in 
the previous paragraphs are specified In terms of typical minimum, nomina , an 

maximum values in Table 7—1 • 


TABLE 7-1 

CONTROL PARAMETER SUMMARY 


Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

Remarks 

a 1 

0.50 

1.0 

1.5 

Heading control gain at 
velocity V q 

V 


200 ft/sec 


Normalizing velocity 

0 

T. 

0.70 

1.0 

2.0 

Heading error filter 

A 

L 1 

20 deg 

30 deg 

40 deg 

Roll displacement com- 
mand limit 


3.0 deg/sec 

5 deg/sec 

8 deg/sec 

Roll rate command limit 

z 

t b 

1.0 sec 

2.0 sec 

4.0 sec 

Mode transition 
smoothing 


g. Performance Criteria 

Heading Control Tr ansient Response 

a. At an aircraft velocity of about 200 to 300 feet per second, 
apply a 45 -degree step heading change command. 

b. The roll rate limit should be held within 110 percent and the 
bank angle Increased to the maximum value (L,) with a maximum 
overshoot of about 3 degrees. 

c. Sideslip should never exceed a value equivalent to a lateral 
acceleration of 0.08g. 

d. Roll out to the desired heading should be achieved with a 
heading overshoot restricted to about 2 degrees. 
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e. When stabilized on the reference heading, apply a step rudder 
command of about 5 to 10 degrees S R for about 5.0 seconds. 
Remove the rudder command and allow the aircraft to settl 
the reference heading. The heading error should be less 
0.5 degrees within 10 seconds after release of the disturb- 
ance. The corrective response should have a maximum of one 

overshoot. 


2. Localizer and Landing Guidance 


a. General 

The ILS steering functions covered in this report relate to the 
capture and tracking of the ILS localizer radio beam from the initial intercept 
path to landing. Two important aspects of the ILS control problem are avoided 
by the simplified statement of control equations given herein. The first in 
volves the terminal area navigation problem associated with establishing the 
proper localizer intercept trajectory. The availability of range and bearing 
to the localizer transmitter (VOR and DME source located near localizer trans 
mitter) could permit the automatic computation of optimum terminal area flight 
paths. For example, the aircraft could be guided toward the proper intercept 
path and then automatically turned to that intercept heading when the downrange 
distance is considered optimum. This terminal area navigation and guidance sys- 
tem would use the aircraft turning radius constraints to compute a fixed bank 
angle maneuver that yields a flight path tangent to the desired localizer inter- 
cept heading. In effect, a circle of radius given by 




g tan 


(7-11) 


(where <t> = the maximum permissible bank angle) could be located so that it is 

tangent to the desired intercept heading. Automatic programs that provide this 
navigation function are not covered in the present report. The problem is 
started on the heading select mode where the selected heading is the desired 
beam intercept angle and the point of intersection is sufficiently far from 
touchdown to permit beam capture without excessive overshoots. 
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The other aspect of the localizer steering problem that is not 
covered in this report involves the considerable complication associated with 
redundancy and reliability considerations for sensors. In a Category III, fail- 
operative control system, the sensors and other reference devices used for flight 
path control must be operating in a properly monitored and redundant configura- 
tion. If some of these sensors fail during earlier phases of the flight, an 
automatically coupled ILS approach should not be precluded. Thus, if the Iner- 
tial Navigation System (INS) provides information used in the control laws, 
should a loss of that INS data eliminate the automatic approach capability? 
Systems in operation today have answered this question by providing back-up modes 
in the event of the loss of such data. The following is a brief summary of typi- 
cal compromises that have been made in operational autopilots in order to cope 
with practical problems of this type : 

• Localizer beams that meet Category II (or higher) standards are 
able to provide adequate lateral velocity information. However, 
the airborne systems must also operate with Category I (or 
poorer) beams, where beam noise precludes the derivation of 
satisfactory lateral velocity data. Hence, autopilots mechanize 
compromise control laws which derive lateral velocity inertially 
from aircraft heading. 

• Systems using heading derived lateral velocity contain large 
crosswind errors. These errors are minimized or eliminated by 
the type of control law used. The crosswind error can be elimi- 
nated if drift angle correction is derived from the INS computer. 
If a localizer steering control law depends upon true velocity 
vector heading, then the loss of the INS will force a breakdown 
of the steering law. Hence, in typical operational systems, 
back-up steering laws are implemented to cope with an invalid 
INS, but these steering laws are only activated upon loss of the 
INS. 

• Localizer gains are programmed downward during the final approach 
phase. The gain programming is often made a function of radio 
altitude. If the radio altimeter is not "valid", a back-up gain 
program is activated. This back-up gain program may be a func- 
tion of time, or it may be a function of marker beacon signals. 
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The activation of an alternate gain reduction program can 
permit penetration to lower altitudes although not to flare- 
out altitudes. 

In this report, the concepts of back-up control laws and alternate 
sensing schemes are not considered. It is assumed that all of the required data 
is available. It is also assumed, however, that the localizer radio signal is 
not of ideal quality. Hence, its use in deriving lateral rates is restricted to 
reasonable applications that minimize noise effects. The dominant lateral, veloc- 
ity source for the specified control laws is drift angle corrected heading. An 
idealized source of gain programming is obtained from range to localizer. Such 
information is not generally available today and a radio altimeter would be a 
more reasonable source of such information. The altimeter, however, is a com- 
promise choice because of its dependence upon a smooth terrain along the final 

approach path. v 

b. Localizer Geometry and Control Phases ..... ~ ■ 

There are five phases of lateral beam steering in regard to control 
laws used. They are: 

1 . Intercept - Heading select mode on a constant beam intercept 
heading. 

2. Capture - Proper penetration of the beam has been sensed and 
the aircraft is turned to align with the beam center line. 

3. On-Course - Alignment has been satisifed with regard to posi- 
tion and rate errors , and tight tracking of the beam is 
initiated . 

4. On-Course/Final Approach - The final phase of the localizer 
tracking is initiated (usually at glide slope penetration) 
and downward gain programming of the beam signal is initiated. 

5. Decrab - The aircraft nose is aligned with the runway center 
line immediately prior to touchdown. 
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Figure 7-3 

Localizer Capture and Tracking Geometry 
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Figure 7-3 illustrates these control phases in terms of the localizer 
beam geometry. Note that for consistently good performance, the intercept range, 

R , should be greater than 12 nautical miles if bank angle limits of about 25 de- 
grees are to be imposed. (For a 10,000-foot runway, an R ± of 12 nautical miles 
represents an intercept about 10 nautical miles from the runway threshold.) The 
intercept angle is It should be about 45 degrees. Intercept angles as large 

as 90 degrees are acceptable if the intercept range is increased to greater than 
14 nautical miles. The capture phase will automatically start when the beam is 
penetrated to a value P q . This capture point is computed as a function of posi- 
tion and rate so that it will vary with distance from the runway and the steep- 
ness of the intercept. For large intercept angles, capture starts near the outer 
boundary of the beam; while for shallow intercept angles, capture is delayed 
until the aircraft is near the beam center (small P Q ) • 

A standard localizer beam for a 10,000-foot runway is assumed for the 
calculation of all nominal parameters. Such a beam is 3.6 degrees wide (7.2 
degrees total). Figure 7^4 illustrates how this converging beam produces an in- 
creasing sensitivity as the runway is approached. From the point of the penetra- 
tion of the glide path center (28,000 feet for 1,500 feet of altitude and a 3.0- 
degree glide path) to touchdown, the localizer sensitivity increases by a factor 

of nearly 5 to 1.0. 
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FT RUNWAY <d - 10,000; e - 1.000) 
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Figure 7-4 

Localizer Position Measuring Sensitivity 
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c. Control Laws 

(1) Localizer Capture 

The control law block diagram for all phases of localizer control 
is illustrated in Figure 7-5. The problem starts on the heading select mode where 
the steering law is 



(7-12) 


where ^ is the velocity vector heading 


ft, = * + t 


and is the selected beam intercept heading. 

s 


(7-13) 


If <// is the runway heading, the 
R 


intercept angle, is 


% - V 




(7-14) 


(Note that the polarity of roll error summing is ^ heading seleCt 

control laws have been described in .equations .(7-2) through (7-9).. 

The aircraft is maintained on the selected heading while the 
capture computer sums weighted beam displacement, beam rate, and heading error in 
accordance with either of the following equations: [Note that (7-1 5a) or (7-1 5b) 

can be optimized for good performance but with different combinations of constants 

c. , c_ and c_.] 


T 2 


c 2 0s 


C 1 sin + (r„s + 1)(r,s + D 


C 3 P 

T . S + 1 

4 


= e 


(7-1 5a) 


or 


c 2 0 s c 3 0 

Cl (1 - cos * E ) + (r 2 s + 1)(r 3 s + iJ + V + 1 


(7-1 5b) 
(7-16) 


where = \ 

and the following sign conventions are observed: * E , fi . fi are positive for those 

conditions that displace the aircraft or cause a rate of change of displacement 
to the right of the beam (from' the viewpoint of the pilot) . The angle P is there- 
fore positive if it is defined as counterclockwise when viewed from the localizer 
transmitter. (In Figure 7-3, \ and P are positive and P is negative at the cap- 
ture angle P . The capture phase starts when e = 0. 
o 
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Figure 7-5 

Localizer Control Block Diagram 

















The heading select mode is disengaged and any residual bank 
angle command decayed to zero in accordance with equation (7-10) . The localizer 
capture control law is now activated. It is 


a^s 


b 1 + (r 2 s + 1)(r 3 s + 1) 


aJ? 

* + TT+T = - *C 


(7-17) 


(2) On-Course (0/C) Tracking 

The capture equation (7-17) aligns the aircraft with the local- 
izer beam. The on-course computer senses beam displacement beam rate and bank 
angle and, when these parameters satisfy certain specified minima, the 0/C phase 
will begin. The 0/C sensing logic equation is 


O/C - 



< e 



< e 


2 



(7-18) 


When the 0/C logic equation is satisfied, the control law 
switches to the 0/C parameters. Note that control variables do not change; only 
the computations performed on these variables change. Hence, there is no need 
for any mode transition synchronization or data smoothing. The 0/C control law is 


T 1 8 

b 1 r 1 s + 1 




(J 2 s + 1 H T 3 s + 1) 


— — - — , + a_ 
r ,s + 1 3 


/ 


0 dt - " 


(7-19) 


As the O/C control law is activated, bank limits are usually reduced from 25 to 
30 degrees to about 10 degrees. Likewise, roll rate command limits are reduced 


(3) On-Course (O/C) Tracking - Final Approach 

This phase starts with penetration of the glide slope center 
line. The control law remains identical to equation (7-19) except that the gains 
are now changed and programmed. Thus, for this final phase, the control law may 

be expressed as: 


b i \ 


r v 


r ^ + ^ + ft dt] 

T-S + 1 
* 

+ 

(r 2 s + 1)(r 3 s + 1) r 4 s + 1 3 J J 


k(R') - - * c (7-20) 
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where the individual constants, a,, a.,, tnay also be modified, but the 
gain reduction program k(R') is activated. The gain reduction is of the form 


k(R') - 1 - 



k 


(7-21) 


where 

R* = distance to runway threshold 

R' = distance from glide slope intercept point to runway threshold 
G/S (see Figure 7-4) 


k = localizer control attenuation factor 
Thus, if k - 0.5, the localizer gain at the runway threshold will be half the 
gain at the glide slope intercept point. 


(4) Decrab (D/C) 


Two techniques 
from work in automatic landing, 
decrab. 


for automatic runway alignment have resulted 
These are the skid decrab and the forward slip 


Skid Decrab 

This technique involves roll control to track the lateral path 
through coordinated turns (aero sideslip down to a "decrab" altitude of approxi- 
mately 8 feet. In the presence of cross winds, a aero sideslip crab angle wall 
develop. At the decrab altitude, the lateral guidance commands are removed and 
aero roll angle is commanded. At the same time, rudder commands are used to 
align the aircraft with the runway heading (decrab) . Predictive commands are 
added to both the rudder and aileron channels to provide surface deflections that 
will compensate for roil and yaw moments resulting from the sideslip developed 
during the maneuver. The system is normally designed so that touchdown occurs 
when approximately 70 to 80 percent of the crab angle is removed. At this time, 
the crab angle is small and the aircraft has a yaw rate established in the di- 
rection of the remaining crab angle. This results in low side forces on the gear 
at touchdown and does not allow time for the aircraft to develop a significant 
cross runway drift velocity. 

Forward Slip 

Hie "forward slip" technique involves aligning the aircraft 
heading with the runway heading by applying roll and yaw commands at an altitude 


20b 


of approximately 200 feet. The roll commands used for lateral guidance combined 
with the rudder commands used for alignment result in a sideslip equal to the 
original crab angle. The vehicle can be landed on one gear truck in the forward 
slip configuration provided maximum bank angle constraints (imposed by wing 
scrape limitations) are observed. In fact, this is the normal manual landing 
technique for transport aircraft. Because of restrictions on roll attitude re- 
sulting from wing and engine pod clearances, techniques have been developed to 
reduce the touchdown roll attitude to an acceptable value with an additional skid 
maneuver. Systems have been developed combining the forward slip and the skid 

decrab maneuver. 

The forward slip maneuver is the preferred runway alignment 
technique for manual control. Pilots have found it easier to minimize lateral' 
drift with this technique than with the skid decrab. For the skid decrab to be 
done properly, a critical and precise sequence of rudder and roll commands must 
occur in the final three seconds prior to touchdown. Automatic systems can, in 
general, perform this maneuver with less difficulty than a pilot because they can 
utilize precise measurements and computations to develop the necessary roll and 

yaw controls. 

For automatic control, the forward slip maneuver has the disad- 
vantage of interacting with the lateral guidance. While rudder control is main- 
tained to keep the vehicle heading aligned with the runway heading, roll commands 
are used to command sideslip for lateral guidance. To avoid large lateral errors 
in the presence of wind shear and gusts, cross feed is required between the rud- 
der channel and the roll channel. Experience has shown that the definition of 
these cross-feed terms is critical and that small errors in these parameters can 
result in lateral guidance errors which are more dangerous than incorrect runway 

alignment at touchdown. 

In this study the skid decrab was selected for implementation 
because of its simplicity. It is acknowledged that pilots tend to prefer the 
other approach for manual control but a skid decrab should be adequate for most 
aircraft when precise automatic control is feasible. The recommended decrab 

control laws are: 

\oM ' \ + + V + ^ <7 ' 22> 

where ip is the runway heading and 8 is the predictive rudder command. 

R 
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(7-23) 



<t - V’ \ 



where 

\p = crab angle at decrab initiate 




a 1 

T yS + l 


P 


+ 6 


AP 


(7-24) 


where 8^ is a predictive or feedforward compensation used to help keep wings 
level in the presence of sideslip build-up or even to drop the wing in the 
direction of the wind to partially compensate for the lateral drift. 

' *AP " k SA ». - V 


where 



The decrab initiate altitude should be adjusted as a function of 
lateral error and vertical velocity. If h^j, is the altitude at which decrab is 
initiated, it should be defined as follows: 


^/C 


h n/ _ + f (A h, tp - tp y, y) 

D/ NOMINAL 


(7-26) 


d. Stability Considerations - Lateral Guidance (Landing) 

The stability of the flight path control loop used to track the loca- 
lizer is described by the dynamics shown in Figure 7-6 . Note that control laws 
using P, P, and are forms of y and y in accordance with the following: 


P 


Z = 1 

R R' + (d + e) 


(7-27) 



(7-28) 





(7-29) 
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Thus, equation 7-19, expressed in terms of y and y, would be 


b 1 

( T ^8 ^ 

+ a ! 1 

+ !iy+!il= -<}> c (7-30) 

V ' 

+ V 

+ R (r 2 s + 1)(r 3 s + 1) 

sy R y R s c 


(Note that equation (7-30) neglects the R part of P .) 

Figure 7-6 simplifies equation (7-30) by assuming that the filters 
do not have a dominant effect on closed-loop stability. The form of the root loci, 
with and without integral control, is shown in Figure 7-7. Since the gains vary 
with R, it is seen that the loop (with integral control) tends toward conditional 
stability when gains are too low. Low gains correspond to large values of R. 

The effects of the low frequency modes associated with the integral part of the 
control law can be minimized by engaging the integrator only when errors are very 
small. The 0/C engage logic is used to ensure that this condition is met. 

The most important factor involving stability of localizer control is 
the increasing gain as R approaches zero [equation (7-30) and Figure 7-4] . The 
gain reduction program used in the 0/C - final approach phase helps compensate 
for this problem. When the damping term (y) is obtained from heading, then the 
gain of that term does not increase as the runway is approached. Hence, the 
ratio of rate to displacement decreases as the displacement gain rises. This is 
one reason why beam rate (a 1 term) is used as well as heading. However, the 
effectiveness of the beam rate damping is limited in bandwidth because of the re- 
quired filtering. Hence, if the total y and y gain is allowed to go too high for 
small R, the closed loop frequency will reach a region where the lags in the beam 
rate filter become destabilizing. The effect of these lags is not shown in 
Figure 7-7. If they were included, the dominant roots would turn into the right- 

half plane as gain is increased. 

e. Control Parameters Summary - Localizer and Landing Guidance 

The control parameters identified in the control law equations given 
in the previous paragraphs are specified in terms of typical minimum, nominal, and 
maximum values in Table 1-7.. 
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Figure 7-7 

Stability Analysis of Lateral Flight Path Control 



Parameter 

Typical 

Minimum 

Value 

Typical 

Nominal 

Value 

Typical 

Maximum 

Value 

C 1 

2.0 

2.5 

3.0 

c 2 

4.0 

7.0 

10.0 


n - R 

1 A TT ^ 

J Q y ^ 

C 3 

u * 3 x 72,000 

72,000 

x 72,000 

r 2 

0.5 

1.0 

1.5 

T 3 

0.1 

0.2 

6.4 



— 

.. -- — 

r 4 

0.10 

0.25 

0.50 

(Capture) 

1.0 

1.5 

2.0 

a 1 (Capture) 

3 

5 

10 

(Capture) 

18 

20 

25 

(O/C) 

1.0 

1.5 

2.0 

a, (O/C) 

3 

5 

10 

a 2 (O/C) 

18 

20 

25 

a 3 (O/C) 

0.20 

0.30 

0.40 

T . 

25 

30 

40 

1 



• 

bj (O/C Final) 

1.8 

2.2 

2.5 

a^ (O/C Final) 

3 

5 

10 

a 2 (O/C Final) 

18 

20 

25 

a 3 (O/C Final) 

0.20 

0.30 

0.40 


k 


0.25 


0.5 


0.70 










TABLE 7-2 (cont) 

LOCALIZER AND LANDING GUIDANCE 
PARAMETER SUMMARY 


. 

Typical 

Typical 

Typical 

Remarks 

Parameter 

Minimum 

Nominal 

Maximum 


Value 

Value 

Value 



1 

2 deg 5 per 

4 

Linear decrab rudder gain 

\ 

deg 


deg 6 r per deg 

c i 

Lj (Capture) 

0.05 

0.10 

0.20 

Decrab integral ratio 

1.0 

1.5 deg 

2.0 

Beam displacement limiter 
(deg P ) 

L 2 (Capture) 

• — 

None 

— 

Beam rate limiter 

L t (O/C) 

1.0 

1.5 deg 

2.0 

Same as above 

L, (O/C) 


None 

■ — 

Same as above 

2 

L 3 (Capture) 

20° 

25° 

30° 

Heading error limiter 

l 3 (0/C) 

20° 

25° 

30° 

Heading error limiter 

Displ (Capture) 

25° 

30° 

36° 

Roll command limit 

L 4 Displ (O/C) 

8° 

10° 

15* 

Roll command limit 

Rate (Capture) 

5 deg/sec 

7 deg/sec 

10 deg/ sec 

Roll command rate limit 

L 4 Rate (O/C) 

3 deg/sec 

4 deg/sec 

8 deg/sec 

Roll command rate limit 


0.15 

0.25 

0.30 

O/C sensor trip logic 

1 



^MAX 

e „ 

0.010 

0.013 

0.015 

O/C sensor trip logic 

2 




^MAX 

e 3 

2 deg 

3 deg 

5 deg 

O/C sensor trip logic <$> 

^/C 

4 

6 

10 

Nominal decrab altitude - 
wheels above ground 


f. Performance Criteria 
(1) Localizer Capture 

(a) Initialize aircraft on a 45-degree beam intercept heading 
12 nautical miles from runway threshold. The localizer should be captured with 
less than 10 percent overshoot. The on-course (0/C) sensor should operate before 
the aircraft is 9.0 nautical miles from touchdown. 
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(b) Initialize aircraft on a 90-degree beam intercept - heading 
14 nautical miles from runway threshold. The same performance as in (a) above 
should be achieved. 

(2) Localizer On Course 

(a) Initialize aircraft so that 0 = 1.0 degree and 0 = 0 at 
6 nautical miles and 2 nautical miles from runway threshold. (Zero crosswind 
should be maintained.) The aircraft should converge to the beam center (0=0) 
with less than 20 percent overshoot. The steady value of 0 should be reduced 
to below 0.36 degree within 20 seconds. The steady-state response should be 
within 0.10 degree by the time the runway threshold has been reached. 

(b) With on-course\ parameters set, initialize the aircraft so 

that 0=0 and 0 = 0 at a distance of 8 miles from the runway threshold. Apply a 
10-knot crosswind step. The peak 0 error should be below 1.0 degree and a steady- 
state value below 0.1 degree should be attained with overshoot held to below 20 
percent of the peak error. The 0 error should be below 0.36 degree within 30 sec- 
onds after insertion of the step crosswind. 

(3) Decrab 

Decrab performance can only be evaluated in conjunction with 

flareout tests. 

(a) Establish the nominal time required for the aircraft to 
descend from the nominal decrab altitude to touchdown. 

(b) Establish a steady-state descent with 5 degrees of crab 
angle. The decrab should be accomplished with the following touchdown criteria. 

• Landing gear drift velocity should be below 
3 feet per second at touchdown. 

• Lateral drift from initiation of decrab should 
be below 20 feet. 
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B. DIGITAL PROGRAM - LATERAL GUIDANCE 
1. Control Law Conversion 

The representation of the lateral guidance block diagram and control 
laws in FORTRAN notation is summarized in Figure 7-8 with a tabulation of the 
FORTRAN namelist given in Table 7-3. 

TABLE 7-3 

VARIABLE LIST FOR SUBROUTINE LGUIDE 


Variable 

FORTRAN Name 

Description 

a 1 

ALG1 

Localizer beam rate gain 

a 2 

ALG2 

Localizer beam displacement gain 

a 3 

ALG3 

Localizer beam displacement integral gain 

b 1 

BL1 

Heading error gain' 

— 

BMPEN 

Glide-slope beam^penetration logic variable 

— 

CB 

Capture localizer beam logic variable 

e 

CBTEST 

Localizer capture trip logic level 

— 

CBTHLD 

Localizer capture trip logic threshold 

R' 


Range to runway threshold 

r g/s 

DMEGS 

Range to runway threshold at glide-slope penetration 

6 1 

EPLB1 

• 

On-course trip logic - P threshold 

€ 2 

EPLB2 

On-course trip logic - P threshold 

€ 3 

EPL3 

On-course trip logic - 0 threshold 

P 

EPSLOC 

Localizer beam displacement in degrees 

— 

HH 

Heading hold logic variable 

C 1 

KCL1 

Localizer capture trip logic: heading error gain 

C 2 

KCL2 

• 

Localizer capture trip logic: P gain 

C 3 

KCL3 

Localizer capture trip logic: p gain 

k(R’) 

KLDCT 

Attenuation factor for localizer gain reduction 
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TABLE 7-3 (cont) 

VARIABLE LIST FOR SUBROUTINE LGUIDE 


Variable 

FORTRAN Name 

Description 

h 

APSI 

Heading hold gain 

L 1 

LGL1 

Localizer displacement signal limiter 

L 3 

LGL3 

Heading error limiter 

V p) 

LGL4P 

Roll command position limit 

l 4 (r) 

LGL4R 

Roll command rate limit 

k 

LOCATN 

Attenuation factor for localizer gain reduction 

p 

max 

LOCMAX 

Localizer beam width in degrees 

1 

— 

ONC 

On course logic variable 

0 

PHI 

Roll angle 

*c 

PHICOM 

Roll angle command , 

— 

PHIMAN 

Manual bank angle command 

0 

PSI 

Heading angle 


PSIDRF 

Aircraft drift angle 

*ref 

PSIREF 

Heading reference angle • - 

*R 

THETRK 

Runway heading 

T 2 

THH 

Synchronizer time constant 

T 1 

TLG1 

Heading error washout filter time constant 

T 2 

TLG2 

• s i 

P filter time constant 

T 3 

TLG3 

• 

P filter time constant 

T A 

TLG4 

P filter time constant 

— 

DT3 

Subroutine sample time 

— 

ANGLE 

Input error angle ) 

> Heading error processor 

— 

ANS 

Output error angle ) 
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PSIDRF 



LGL4P 

LGL4R 

HH 

30° 

50/SEC 

CB 

30° 

70 /SEC 

ONC 

10° 

40 /SEC 



(EPSLOC < 0.2°) • ONC 


ONC-1 IF(CB-I) 
•(|EPLDOT| < EPL2) 

• {|EPLF| < EPL2) 

• ( | PHI 1 < EPL3)] 


Figure 7-8 

Lateral Guidance Block Diagram - 

Fortran Notation 
ro 





















2. Program Flow Chart 


The sequences of operations associated with the LGUIDE subroutine are: 

a. Initial Condition Calculations 

b. Heading Hold Gain Calculation 

KPSI = APSI/200.0 

c. Localizer On-Course Test Thresholds 

EPL1 = EPLB1 * LOCMAX 

EPL2 = EPLB2 * LOCMAX 

d. Difference Equation Coefficients 

Heading Error Washout Filter 
CPSI = EXP (-DT3/TLG1) 

Localizer Beam Displacement Rate Filter 

CEPL1 - EXP (-DT3/TLG2) + EXP (-DT3/TLG3) 

CEPL2 = EXP [ -DT3 * (1.0/TLG2 + 1.0/TLG3)] 

DEPL2 = [EXP (-DT3/TLG2) - EXP (-DT3/TLG3)] 

* (1 .0/TLG2/TLG3) / (1 .0/TLG3 - 1.0/TLG2) 

Localizer Beam Displacement Filter 
CEPF = EXP (-DT3/TLG4) 

DEPF - 1.0 - CEPF 

Roll Angle Command Synchronizer Washout Filter 
CPH1 - EXP (-DT3/THH) 

Localizer Beam Displacement Integrator 
DEPLI = ALG3 * DT3 
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The control law parameters for localizer on-course tracking are summar- 
ized on the figures depicting the transient responses. The localizer capture 
parameters that gave best results are: 

LGL3 =25° 

BL2 = 3.0 
TLG1 = 40 sec 
KCL1 =2.5 
KCL2 =2.8 

KCL3 = 1.05 * (12 NM/R) 

TLG2 =1.0 
TLG3 = 0.2 
TLG4 = 0.25 
ALG1 = 50 
ALG2 =50 

ALG3 = 0.2 ' 

The flow diagrams for lateral* guidance are summarized on Figures 7-9 (a) 
through 7-9 (h). * 
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Figure 7-9b 
LGUIDE Flow Chart 
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Figure 7-9d 
LGUIDE Flow Chart 







Figure 7-9e 
LGUIDE Flow Chart 
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Figure 7-9f 
LGUIDE Flow Chart 
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Figure 7-9g 
LGUIDE Flow Chart 










C. SIMULATION TEST RESULTS 


1. Heading Control 

Figure 7-10 illustrates heading control transient responses which meet 
the performance criteria specified in section 7-A-1f. A 5-degree step rudder 
held for 5 seconds causes a yaw transient of 5.5 degrees. Ninety percent of the 
recovery is achieved in about 6 seconds and the response is well damped. 

Heading command response is almost entirely defined by the process kine- 
matics and the control system constraints as illustrated in Figure 7-10. The 
initial response to a step 45-degree heading change command is dictated by the 
5-degree-per-second roll rate command constraint, the 30-degree bank command 
limit and the coordinated turn rate achieved by the 30-degree bank angle. When 
15 degrees of heading error remains, the control comes out of the command limit 
and the linear system takes over to provide a well-damped acquisition of the 
reference heading. The peak miscoordinaXidnis 0.03g, well within the specified 

criteria. - _ ■ 

2. Localizer Capture 

Excellent capture responses for a 90-degree and a 45-degree intercept are 
illustrated in Figure 7-11. In the 90-degree intercept case the overshoot is an 
imperceptible 0.05 degree of beam while in the 45-degree intercept case, on- 
course control starts when the beam error is about 0.1 degree and the overshoot 
is less than 0.05 degree. In both cases the overshoots would not be perceptible 

on the pilot's display instrument. 

3. Localizer On-Course Transien t Responses 

The fact that accurate localizer control is more a measurement than a 
control problem is illustrated by the responses in Figure 7-12. The aircraft is 
initialized parallel to the beam but offset 1 degree of beam angle. The recovery 
for this position error 6 nautical miles and 2 nautical miles from touchdown is 
shown on this figure. Case © includes the washout filter in the control law 
while case © eliminates the washout filter from the heading feedback term. 
Obviously, when we eliminate the washout we can obtain an excellent, almost dead- 
beat, response. The washout degrades the transient response with an effect 
analogous to increasing the gain of the integral loop. It is used to compensate 
for inaccuracy in the measurement and computation of heading. In a lateral 
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guidance law with lateral velocity weighted 20 times lateral displacement 
(y + 20 y ** K# ), a 1.0 foot per second error in measuring y can cause a 20 foot 
offset in y. Note that 1.0 foot per second lateral velocity error is about 0.2 
degree of heading error at approach speeds. If heading can be measured to 0.05 
degree (or lateral velocity to about 0.25 foot per second) then the washout will 
not be needed. With velocity vector heading accuracies of 0.05 to 0.2 degree, 
the washout may be eliminated if sufficient beam rate compensation can be used. 
The restriction on the use of beam rate (ALG2 in Figure 7-8) is related to beam 
noise and the type of servo system used in the autopilot. 

A s umma ry of transient responses at 8 nautical miles, 6 nautical miles, 
and 2 nautical miles from touchdown is given in Figure 7-13. All of these re- 
sults show adequate responses ‘ that should allow convergence of all errors to 
within a few feet of the runway centerline at touchdown. These results were ob 
tained by eliminating the washout from the heading feedback term. 
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Figure 7-10 

Heading Control Transient Response 
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Localizer Capture Responses 
( 90 ° and 45 ° Intercepts) 





LOCALIZER DEFLECTION 
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Figure 7-12 

Localizer Transient Responses with 
and without Heading Washout Filter 
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Figure 7-13 

Localizer On-Course Transient Responses 
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LOCALIZER DEFLECTION 



Figure 7-13 

Localizer On-Course Transient Responses (cont) 




SECTION VIII 

MODE SELECTION AND DATA ENTRY 


A. INTRODUCTION 

The subject of mode selection and data entry for a digital autopilot involves 
operational procedures and data transmission programming that is beyond the scope 

this report. Nevertheless some comments on these requirements are pertinent. 

The flow chart for a mode select panel is one of the more complicated parts 
of a digital autopilot. A strategy for scanning the status of all input-output 
devices each computation cycle must be established. In a practical mechanization 
the data coding scheme and transmission format to and from the panel is an essen- 
tial part of the system design. The versatility, of the data entry and display 
devices enters into the software design requirements. In this discussion, the 
operational requirements for the altitude select and vertical speed select func- 
tions are presented as typical examples of procedures and techniques that are 
applicable to this problem. In place of the panel and numerical keyboard concept 

described, one could easily envision an alphanumeric keyboard and CRT interactive 
terminal. 

B. OPERATIONAL CONSIDERATIONS AND MODE LOGIC 

To utilize the various guidance modes described in this report, the pilot 
must have a means of entering desired references and a means of selecting a de- 
sired control mode. The elements of such a mode control panel for accomplishing 
the non-landing vertical guidance mode functions are shown in Figure 8-1. The re- 
quirements are a keyboard for entering numerical references, a display for trans- 
ferring a verified keyboard entry to the appropriate control law parameter, and 
mode engage switches. For example, on Figure 8-1, the following procedure is 
used to set an altitude reference and then proceed to the desired altitude: 

• The desired altitude is 12,500 feet. 

• Set 1-2-5-0-0 on the keyboard. 12500 will appear on the keyboard display. 

• Depress SET button under altitude reference window. 

• 12500 transfers to reference altitude window (also is entered as h 1 

• Depress altitude select mode switch. 
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Figure 8-1 

Mode Select Panel Requirements For 
Vertical Guidance Modes 
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The aircraft will be maneuvered toward the desired altitude automatically in 
accordance with the following mode sequencing logic: 

• If altitude hold was previously engaged, that mode is auto 
matically disengaged and the vertical speed command mode is 
automatically engaged. 

• A reference nominal vertical speed is automatically selected. 

A suggested nominal vertical speed is the one that points 
toward the reference altitude with a 2.0-degree flight path 

angle . 

• The vertical speed command mode remains engaged until the Ah 

of the altitude select [ Section V, equation (5-16)] is attained. 

• The altitude select armed light is illuminated following de- 
pression of the altitude selecf engage switch. When capture 
to the reference altitud£_is initiated, the ."armed" light is 
extinguished and the altitude select mode is now considered 

engaged . 

• If the pilot wishes to change the vertical speed reference from 
the value automatically selected, he need only follow the normal 
procedure of entering a new vertical speed via the keyboard and 
the reference vertical speed SET button. 

• If the altitude hold mode were not initially engaged when the 
altitude select switch was depressed, an automatic vertical speed 
steering program will not be commanded if the aircraft is already 
heading toward the desired altitude. However, if the aircraft is 
moving away from the desired altitude, then the vertical speed 
command sequence described above will be activated. 

It is apparent that many computation and logic sequences must be programmed 
in the general-purpose digital autopilot computer to accomplish these functions 
included in these programs are the logic and computations that examine the 
reasonableness of the data entries (selected altitudes above allowable opera- 
ting limits,' for example) and the compatibility of selected modes (altitude 
select or altitude hold cannot be engaged while in an automatic approach and 



landing sequence). The displays and warning sequences following mode select or 
data entry errors are an important part of the autopilot program. 

The other functions which must be handled in a manner similar to that de- 
scribed for altitude select and vertical speed select are: 

• Heading select 

• Airspeed select (for autothrottle or pitch control modes) 

• VOR course select 

• Various mode engage functions. 

Proper programming of these functions are dictated by the cockpit display and 
control concepts which are beyond the scope of this report. 
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APPENDIX A 

DESCRIPTION OF DIGITAL AUTOPILOT PROGRAM 

The stabilization and control subroutines are wzitten in FORTRAN IV. This 
makes the program less efficient, but makes it independent of the computer, and 
also makes it more accessible to the design engineer. 

Farh subroutine is organized in three different sections, initial conditions 
(IC), engage, and operate. The linear filter portions of the control laws are 
calculated in the same manner as described in Reference 5 by using difference 
equations. The IC calculations include the difference equation and rate limiter 
coefficients, which are determined as a function of the filter parameters and 
the sampling time interval. When a control mode is engaged, the engage section 
of the corresponding subroutine is entered once to initialize filters and inte- 
grators. The operate section of the control subroutine is entered once each 
loop time interval and contains the control logic and filters. A list of sub- 
routines is given in Table A-1 . 
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TABLE A-1 
SUBROUTINE LIST 


PITCHS 

Longitudinal stabilization and 
automatic trim 

AUTOTH 

Autothrottle 

LATSTB 

Lateral stabilization 


Roll stabilization 

Yaw Damper 

Turn coordination 

ALTHLD 

Altitude hold 

VERTS C 

Vertical speed hold 
Vertical speed select 
Altitude capture 

LGUIDE 

Heading hold 
Heading Select 
Localizer capture 
Localizer tracking 
Decrab 

VTLAND 

Glide-slope capture 
Glide-slope tracking 
Flare 

MEASURE 

Compensated vertical speed computer 

HDTCMP 

Vertical speed command processor 

HEP 

Heading error processor 

SASIC 

— — ■ ■ — 

, Difference equation coefficient 
calculations for most control modes 




APPENDIX B 


DIGITAL SIMULATION OF REFERENCE JET TRANSPORT (RJT) 

The reference jet transport simulation used as the test vehicle for the auto- 
pilot design and evaluation is that of a current transport aircraft with the sub- 
stitution of hydraulically powered surfaces in place of aerodynamically boosted 
controls used in the actual aircraft. Elevator, aileron, rudder and throttle 
control servos are modeled as 2nd order systems with a damping ratio of 0.7 and 
a natural frequency of 20 rad/sec. The surface power boost actuators are modeled 
as first order lags with time constant of 0.067 second. The engine response is 
modeled as a first order lag with a time constant of 1.25 seconds. Flaps are 
operated at a fixed rate of 2 degrees /second and the horizontal stabilizer is 
driven (by the automatic system only) at a fixed rate of 0.06 degree/second. 

The simulation includes ground effect, landing gear dynamics, and uses poly- 
nominals for non-linear aero dynamic derivatives. 

The equations of motion are programmed in FORTRAN IV to run in nonreal time 
on the IBM 360/67, with line printer output. The simulation is also available on 
the EAI 8400 computer for real time flight evaluation, which includes a full 
scale cab with motion simulation and D/A outputs of all variables of interest. 

The digital program is divided into three basic loops. They consist of: 

1 . Body axis accelerations , Euler angle rates , transformation elements , 
part of landing gear. 

2. Rotational dynamics, and remainder of landing gear. 

3. Translational Dynamics 

For the EAI 8400 an executive timing routine controls the order and point in 
time of execution of these loops. The relationship of the timing of these loops 
can be varied: At the present time the relationship used is 1:1:2. (0.05 sec 

:0.05 sec :0.1 sec). 
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The automatic control system is interfaced with the RJT by FORTRAN calls from 
the second and third loops. The actual location depends upon the particular 
function of the variable being considered. No interface has yet been provided in 
the simulator cab for pilot selection of modes or reference control variables. 

At this time, mode control is accomplished by means of special purpose programs 
for testing individual control modes and mode transitions. 
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APPENDIX C 

RELATIONSHIP BETWEEN ZERO SIDESLIP 
AND ZERO LATERAL ACCELERATION 


From Reference 2, Appendix B or any text on aircraft stability and control 


The aircraft lateral acceleration (a y ) ^- s: 


a ° g sin 0 

y 


+ SS ( c 0 + C 6 \ - 1 T fl 

“ l yp y« R R J m p 


(C-D 


or 


a y = g sin 9 + ^ « R + f- ^ ^ 


R 


(C-2) 


A body-mounted lateral accelerometer provides an output proportional to A y . 


A » a - g sin 0 
y y 

The condition for P - 0 [ from Equation (C-2)] is 


(C-3) 


a - “ 8 sin 0 + ^ C,_ 5 r 


QS 

m y&. 


(C-4) 


The condition for A y - 0 is [ from Equation (C-3)] 

a y = g sin 0 


(C-5) 


Therefore, the condition for simultaneous equality of zero sideslip and zero 
lateral accelerometer output is : 


C 8=0 
R 


(C-6) 


This condition is approached for aircraft with large stable values of C 
(positive) where sideslip is near zero in a turn with zero rudder. 
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